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Effect of Sodium on Mechanical 
Properties and Friction-Wear Behavior 
of LMFBR Materials 


By K. Natesan,* T. F. Kassner,* and Che-Yu Lit 


Abstract: The status of materials—sodium compatibility in- 
formation for the liquid-metal-cooled fast breeder reactor is 
reviewed, with emphasis on mechanical behavior. An attempt 
was made to correlate the creep-rupture properties of un- 
stabilized austenitic stainless steels exposed to sodium with 
available data obtained in air or inert-gas environments based 
on the composition, microstructure, and thermal history of the 
materials. Information on the degradation of the mechanical 
behavior due to fast-neutron irradiation is included for 
comparison with the data in a sodium environment. The 
adhesion, wear, and friction properties of materials in sodium 
are also reviewed, Requirements are summarized for additional 
materials-properties information in a sodium environment to 
reduce the uncertainties in the safety factors applied to present 
design data, 


A number of iron- and nickel-base alloys are being 
considered as construction materials for major compo- 
nents in sodium-cooled fast breeder reactors. Austen- 
itic stainless steels, including the stabilized types, are of 
interest for the fuel cladding, in-core structural compo- 
nents, piping, and intermediate heat exchanger. Sta- 
bilized and unstabilized ferritic steels, with low 
chromium and molybdenum contents, and nickel-base 
alloys have been suggested for the steam generator in 
the secondary sodium coolant system. The selection of 
a particular material and combinations thereof which 
will be used in both the primary and secondary circuits 
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is based upon a number of factors. For in-core 
applications the stability of the material in the fast- 
neutron environment (e.g., swelling behavior) and 
compatibility with the sodium coolant and the fuel 
(including possible fission-product interactions) are of 
primary concern. In the case of the steam generator, 
the resistance of the material to aqueous chloride and 
caustic cracking, as well as minimization of corrosion 
wastage rates due to sodium-water reaction products 
in the event of a tube failure, have an important 
bearing on the material selection. 

In general, mechanical-property requirements for a 
liquid-metal-cooled fast breeder reactor (LMFBR) de- 
sign are the tensile, creep, fatigue, and creep-fatigue 
interaction properties as a function of temperature, 
type of loading, and test environment. The effect of 
strain rate on the ultimate tensile and yield strengths 
and fracture ductility is of interest. Uniaxial and 
biaxial creep-rupture data, where the entire load-strain 
history is determined, yield information on the strain 
and strain rate for the primary, secondary, and tertiary 
stages of creep. Low-cycle fatigue and creep-fatigue 
interaction behavior includes the effects of strain rate, 
strain amplitude, mean strain, and stress-strain hold 
times on the cumulative damage. This information is 
essential to set limits on the service condition of a 
reactor component in design analysis as well as to 
specify failure criteria. 

Most of the available mechanical-property data on 
LMFBR materials have been obtained in air or inert-gas 
environments under basic loading conditions. In com- 
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parison, limited data have been generated in either the 
liquid sodium or fast-neutron environments, and a 
knowledge of the integrated effect of the sodium and 
irradiation environments on the mechanical behavior at 
high temperatures is particularly lacking. As a result, 
various safety factors are incorporated into design 
practice to take into account the various environmental 
effects. For this purpose the materials-data require- 
ments can be conveniently classified into in-core (e.g., 
fuel cladding, subassembly duct, and core structure) 
and out-of-core (e.g., heat exchanger and piping 
components. The combination of complex mechanical 
loading and variations in the total reactor environment 
for different components therefore results in a wide 
spectrum of materials-properties data needed for design 
analysis and the establishment of realistic failure 
criteria. Reactor safety analysis dictates that special 
consideration be given to property data pertaining to 
abnormal operating conditions (e.g., high-temperature 
and high-strain-rate effects). It is also recognized that a 
linear superposition of the fractional damage at ele- 
vated temperatures from basic types of loads is not 
sufficient to account for damage that occurs under 
prototypic operating conditions. Therefore future 
mechanical-property studies must take into account 
the extent and nature of damage that result from load 
interactions and synergetic environmental effects both 
for immediate LMFBR design and for a check on 
proposed empirical formulations. This will enable a 
more critical evaluation of safety factors that are 
applied to existing design data, such as those used in 
ASME codes for out-of-core components. Failure 
criteria for in-core components, such as the fuel 
cladding, have not been well established; however, a 
recent evaluation of possible failure mechanisms has 
been made (e.g., a thin-wall tube can fail by mechani- 
cal instability in the form of a violent rupture or by 
grain-boundary mechanisms in the form of a pinhole 
leak’). 

The purpose of this review is to summarize the 
status of materials-sodium compatibility information 
for the LMFBR with emphasis on the mechanical 
behavior of materials. An attempt was made to 
correlate existing data for austenitic stainless steels 
exposed to sodium with property data obtained in air 
or inert-gas environments based upon the composition, 
microstructure, and thermal history of the material. 
Information on the degradation of the mechanical 
behavior due to fast-neutron irradiation is included for 
comparison with data obtained in a sodium environ- 
ment. The friction and wear properties of materials in 
sodium are also reviewed. 
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REVIEW OF MECHANICAL-PROPERTY 
DATA 


Various degrees of metallic and nonmetallic ele- 
ment mass transfer have been observed when materials 
have been exposed to liquid sodium, depending upon 
the conditions of temperature, temperature gradient, 
sodium velocity and purity, alloy composition, and the 
construction materials of the loop in which the 
experiments were conducted. In many of these experi- 
ments, the mechanical properties of the materials were 
also investigated, and in some instances the properties 
were compared with controlled tests in air or inert-gas 
environments. In this section the compositional and 
microstructural changes and the mechanical behavior 
of materials exposed to sodium will be reviewed and 
correlated with results on similar materials obtained in 
other environments. 


Mass Transfer in Liquid Sodium 


The mass transfer of metallic elements from austen- 
itic stainless steels in liquid sodium is well docu- 
mented.” '® In these studies the effects of temperature 
between 450 and 750°C, oxygen content of the 
sodium from 2 to 30 ppM, sodium velocity between 
0.5 and 40 ft/sec, and specimen position (i.e., down- 
stream effects) on the corrosion rate in the highest 
temperature region of sodium loops, as well as the 
influence of minor alloying variations, have been 
investigated. In general, nominal linear metal-loss rates 
were found after an initial period of less than 
~1000 hr; the corrosion rate increased exponentially 
with temperature and linearly with both the oxygen 
content and sodium velocity to ~10 to 15 ft/sec and 
then became velocity independent at the higher values. 
The downstream or positional effect, metallurgical 
condition (e.g., annealed or cold worked), and minor 
variations in the nickel, chromium, niobium, or tita- 
nium contents of the stainless steels have relatively 
little influence on the corrosion rates in comparison 
with the other variables. 

The compositional changes near the surface of 304 
and 316 stainless steels exposed to sodium at tempera- 
tures between 600 and 677°C in several mass-transfer 
loops are shown? in Fig. 1. The 10 to 20-y-thick 
modified region of the specimens results in the 
formation of a ferrite layer, which, along with the 
corrosion rate, attains a steady-state value after the 
initial corrosion period. At higher temperatures the 
thickness of the modified surface region increases as 
does the corrosion rate.'” Metallic element mass trans- 
fer influences the design of fuel elements and reactor 
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(a) MASS-TRANSFER LOOP NO. 1 
TEST MATERIAL: 304 S.S. 
TEST CONDITION 650°C 

<5 ppM O 

1430 hr 


(c) CULCHETH LOOP 
316 S.S. 


LOOP MATERIAL: 304 S.S. 
TEMP. 600°C (max.) 

LOOP AT: 220°C 
OXYGEN: 5 ppM 
VELOCITY: 1.7 m/sec 
TIME: 1575 br 
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Fig. 1 Metallic element concentration profiles in 304 and 316 
stainless steels upon exposure to flowing sodium in mass- 
transfer experiments.’ 


maintenance operations as a result of the thinning of 
the fuel cladding and deposition of radioactive cor- 
rosion products at various locations in the primary 
sodium system. 

Nonmetallic elements such as carbon and nitrogen 
are known to migrate in sodium loops constructed of 
austenitic stainless steels as a result of the chemical 
activity differences that occur in nonisothermal sys- 
tems. Transfer of these elements also occurs in loops 
where combinations of materials of different composi- 
tion are used (e.g., austenitic and ferritic steels). 
Carbon transfer has been investigated in greater detail 
mainly because the stainless steel used in the core of 
reactors has very low nitrogen contents (<0.01 wt.%) 
to minimize helium generation by means of the n,a 
reaction with nitrogen. 
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Both decarburization and carburization of austen- 
itic stainless steels with initial carbon contents of 0.05 
to 0.07 wt.% have been observed. Carbon concentra- 
tions in these materials ranged from 0.01 (Refs. 18 and 
19) to 3.8 wt.% (Refs. 20 and 21) upon exposure to 
flowing sodium with different carbon activities in the 
temperature range from 600 to 750°C. Stainless steels 
generally lose carbon at temperatures between 700 and 
750°C in sodium loops that operate with a 
AT,*:':19:2223' whereas the material in the lower 
temperature regions tends to carburize. Carbon loss 
from stainless steels at these temperatures also occurs 
in isothermal sodium when a carbon sink such as 
zirconium or vanadium is present in the system.??° 


The microstructural changes that occur in austen- 
itic stainless steel during decarburization in sodium?” ** 
are primarily the dissolution of grain-boundary car- 
bides and enhanced precipitation of the sigma phase, 
which is an Fe-Cr intermetallic compound. 

Carburization of austenitic stainless steels has been 
investigated at temperatures below 750°C for a range 
of carbon concentrations in sodium.?*!%?0?6%? De- 
pending upon the temperature and carbon activity in 
sodium, precipitation of different carbide phases in 
stainless steel occurs during carburization. Bogers and 
Chirer?® and Hiltz”° reported precipitation of 
(Cr,Fe).3C, carbide when the carbon concentrations in 
sodium were low. At high carbon concentrations, 
Hiltz”° and Anderson and Sneesby”® showed that the 
(Cr,Fe),C; carbide formed. Valibus? and Andrews 
et al.?” found that the larger the degree of carburiza- 
tion of stainless steel the richer the carbides become in 
iron. Because of the difference in the temperature 
dependence for grain boundary and volume diffusion 
of carbon, the relative amount of carbide in the grain 
boundaries in relation to that in the bulk material 
increases as the temperature decreases during carburiza- 
tion. 

Equilibrium carbon concentrations in the range of 
0.1 to 1.2 wt.% for several stainless steels exposed to 
flowing sodium with different carbon activities at 
750°C have been reported.*"** Sodium samples ob- 
tained at temperatures of 750 and 400°C during these 
experiments gave carbon concentrations between 0.4 
and 2.5ppM by an oxyacidic flux combustion 
method.*? Although these analyses were independent 
of sampling temperature and were much lower than 
those previously reported?®?’ (i.e., 15 to 30 ppM) for 
experiments of this type, the values could not be 
correlated with the carbon concentrations in. the steel 
because of the limit of sensitivity of the sodium 
analysis technique. However, a good correlation was 
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obtained between the response of a UNC diffusion- 
type carbon meter ** and the equilibrium carbon 
contents of alloys such as 304 stainless steel, as is 
shown in Fig. 2. 
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Fig. 2 Relation between the equilibrium carbon concentration 
in 304 stainless steel and the net response from a diffusion- 
type carbon meter in flame-ionization-detector (FID) units for 
a series of experiments in flowing sodium at 750°C. 


A thermodynamic-kinetic analysis of carbon- 
transfer phenomena in sodium-steel systems which is 
consistent with many of the experimental observations 
in the literature has been reported.*”*’ Equilibrium 
carbon concentrations in austenitic stainless steels were 
computed as a function of temperature, carbon con- 
centration in sodium, and alloy composition. The 
kinetic aspects of the carburization—decarburization 
process were examined in terms of carbon diffusion 
and carbide precipitation in the alloys. The results of 
this analysis provide the technical basis for many of the 
isolated observations of different investigators and lead 
to further predictions of carbon transfer in sodium 
systems. 

Nitrogen transport in sodium loop systems has not 
been investigated so extensively as carbon migration; 
however, the same principles govern both processes. At 


K. Natesan, T. F. Kassner, and Che-Yu Li 247 


present, thermodynamic information on nitrogen in 
stainless steel and sodium is not sufficient to make a 
detailed analysis of nitrogen transfer in loop systems. 

The nitridation of LMFBR materials in the cover- 
gas region of sodium systems has been recognized as a 
potential problem.***? For example, the argon cover 
gas of the EBR-II primary sodium system typically 
contains from 0.25 to 0.50 vol.% nitrogen as an 
impurity. During fuel-handling operations the nitrogen 
concentration of the cover gas may approach ~ 1 vol.%. 
Experiments have been performed to evaluate the 
extent of nitriding of 304 stainless steel at tempera- 
tures of 370, 550, and 650°C in sodium and in sodium 
vapor-Ar-N, gas mixtures with 0.1 and 1.0% nitrogen. 
Chemical analyses of the 5-mil foil specimens exposed 
to 370°C sodium with an Ar-1.0% N, cover gas for up 
to 96 days indicated no increase in the nitrogen 
content of the material; however, a slight amount of 
surface nitriding could be detected metallographically. 
Specimens exposed to the sodium with an Ar-0.1 and 
1.0% N, cover gas at the higher temperatures showed 
negligible weight gains (<0.01 mg/cm?) after 14 days, 
no measurable increase in the nitrogen concentration, 
and only slight nitriding at the grain boundaries. 
Similar foils exposed for 7 days at 550°C in a sodium 
system with a nitrogen cover gas revealed a dense, 
continuous nitride layer to a depth of 0.4 mil for the 
material in the sodium and 0.7 mil for the foils in the 
sodium vapor-nitrogen cover-gas region. 

Grosser*’ investigated the nitridation behavior of 
ferritic and austenitic steels in pure nitrogen, and 
nitrogen that contained limited amounts of oxygen and 
hydrogen both in the presence and absence of sodium 
vapor at temperatures between 300 and 560°C. The 
ferritic Fe-2.25 wt.% Cr-1 wt.% Mo-Nb steel and 304, 
321, and 347 stainless steels showed no increase in 
nitrogen content by metallographic or chemical analy- 
ses after 560°C exposures to the gas mixtures for times 
between 500 and 2000hr. In the sodium vapor- 
nitrogen atmosphere, measurable nitridation of the 
austenitic and ferritic steel started at 370 and 400°C, 
respectively. At 560°C the surface nitrogen content of 
347 stainless steel was 0.8 wt.%, whereas that of the 
ferritic steel was 0.2 wt.%. The results of these studies 
indicate that both the temperature and purity of the 
cover gas are important to establish the nitridation 
behavior of materials in sodium systems. 

Information on the corrosion behavior of nickel- 
base alloys has been obtained primarily from compara- 
tive corrosion tests in investigations of stainless-steel 
mass transfer in liquid sodium.°*”” In general, the 
nickel-base alloys exhibit higher corrosion rates and 
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thicker regions of modified composition near the 
surface,*** but the compositional changes do not 
always result in the formation of a ferrite layer. 
Nickel-base alloys have also been found to lose carbon 
to a lesser extent when exposed to high-temperature 
sodium in AT loops.?”** 

Corrosion rates of ferritic steels in liquid sodium 
are similar to those for austenitic stainless steels.°** As 
a result of the high carbon activity in the low Cr-Mo 
steels, the alloys decarburize at temperatures above 
~475°C in both static and dynamic sodium systems 
constructed of austenitic stainless steel.*°*' The effect 
of stabilizing elements (e.g., titanium, niobium, zirco- 
nium, tantalum, and vanadium) on decarburization of 
the alloys in sodium has also been investigated. *" *? 
Titanium and niobium were found to be particularly 
effective in minimizing decarburization of the alloys in 
sodium. This information has been reviewed in greater 
detail by Bauer and Conner. 


Mechanical-Property Data in 
Liquid-Sodium Environment 


The tensile, creep-rupture, and fatigue properties of 
304 and 316 stainless steels in sodium have been 
investigated to a limited extent where temperature, 
sodium purity and velocity, and the mechanical and 
thermal treatment of the specimens were the major 
variables. A comparatively limited amount of mechani- 
cal-property data for stabilized austenitic stainless 
steels,7°°%° mnickel-base alloys,?”*° and _ferritic 
steels’? 55557 has been obtained in a sodium environ- 
ment. 


Tensile Properties 


The effect of carburization of 304 and 316 
stainless steels in liquid sodium on the tensile proper- 
ties?!:?9:5°58 showed that the ultimate tensile strength 
and the proof stress increased as carbon concentration 
increased in the materials at all test temperatures. For 
carbon concentrations in the alloys greater than 
~0.5 wt.%, the embrittlement was so severe that 
fracture occurred with no measurable ductility at 
temperatures”? between 500 and 800°C. The 0.3% 
proof stress and tensile elongation data”? as a function 
of carbon content in 316L stainless steel are shown in 
Figs. 3 and 4, respectively. The tensile properties of 
stainless steels and nickel-base alloys decarburized to 
carbon levels between 0.005 and 0.05 wt.% have not 
been determined. The tensile properties of ferritic 
steels deteriorate as a result of decarburization in a 
sodium environment.*”*° 
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Creep Properties 


The effects of carburization and decarburization on 
the creep and stress-rupture properties of 304 and 316 
stainless steels have been investigated in both static and 
flowing sodium systems. For conditions of temperature 
and sodium purity (i.e., carbon concentration) under 
which no compositional changes occur in the materials, 
the stress-rupture properties are essentially the same as 
those in a helium or argon environment.*” ® This type 
of behavior is shown in Figs. 5 and 6 (Ref. 59) for 304 
stainless steel with 10 to 15% cold work; similar results 
were reported for 316 stainless steel.°”® A compari- 
son™ of the biaxial stress-rupture properties of cold- 
worked 304 and 316 stainless steels in Fig. 7 shows 
that the former has better strength at temperatures 
below 540°C, whereas the latter is stronger at tempera- 
tures above 650°C. The microstructures of the steels 
exposed to static sodium in the above experiments at 
650°C for 5000 hr showed grain growth and carbide 
precipitation but no evidence of sigma-phase forma- 
tion. 

Limited data are available for the creep-rupture 
properties of austenitic stainless steels under carbu- 
rizing conditions in liquid sodium. The creep-rupture 
strength of 304 and 316 stainless steels in flowing 
sodium at 650°C was reported®® to be the same as that 
in helium; however, the carburization in these experi- 
ments was confined to a 2- to 3-mil surface region of 
the 60-mil-thick specimens for exposure times less than 
~2000 hr. If thinner specimens or longer exposure 
times were used, the alloys would have carburized in 
their entire cross section, and significant property 
changes, when compared with the tests in helium, 
would have been observed. 

The stress-rupture life of several austenitic stainless 
steels has been determined under decarburizing 
conditions in static sodium gettered by zirconium and 
in high-velocity sodium in loops with a large AT. The 
experiments in the zirconium-gettered sodium 
environment resulted in a decrease in the stress-rupture 
life, compared with control tests in helium, as shown in 
Figs. 8 and 9 for 304 and 316 stainless steels at 760 
and 732°C, respectively.4 The decrease in the 
stress-rupture properties was attributed to the change 
in the microstructure that resulted from the loss of 
carbon from the alloys to the zirconium getter. The 
carbon concentration in the alloys decreased from an 
initial value of 0.06 to 0.03 wt.% for 304 stainless steel 
and 0.04 wt.% for 316 stainless steel. There was no 
transfer of metallic elements from the steels in these 
experiments in static sodium. Metallographic 
examination of the specimens indicated that a 
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Fig. 3 Variation of 0.3% proof stress with carbon content for 316L stainless-steel tensile specimens 
carburized in sodium.?°? e, 500°C tests. x, 600°C tests. ©, 650°C tests. 


significant amount of sigma phase was present in the 
material decarburized in the sodium environment. 

The biaxial stress-rupture properties of these steels 
obtained in flowing sodium at 705°C in a AT loop, as 
well as control tests in an argon atmosphere,”?? are 
shown in Figs. 10 and 11. The decrease in the rupture 
time for a given stress in the sodium environment is 
mainly a result of decarburization of the test specimens 
and the associated microstructural changes in the 


materials. The bulk carbon contents of 304 and 316 
stainless steels decreased from initial values of 0.062 
and 0.073 to 0.017 and 0.038 wt.%, respectively, after 
4838 hr at 705°C. A substantial loss of ductility in 
these materials was also observed. The presence of a 
large amount of sigma phase in the materials was 
attributed to the loss in carbon. 

As in the case for the tensile properties, few data 
are available on the creep-rupture behavior of stabilized 
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austenitic stainless steels,” nickel-base alloys,®* and 


bg 64: “ . 
ferritic steels>©™ in a sodium environment. 


Fatigue Properties 


The fatigue properties of LMFBR materials in 
liquid sodium have been investigated to a very limited 
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extent.275*° In Fig. 12, cyclic strain vs. cycles to 
failure for 316 stainless steel tested in air, helium, 
low-oxygen sodium, sodium containing carbon, and for 
specimens preexposed for 272hr to 650°C 
sodium containing carbon is plotted.?’ The low-cycle 
fatigue results indicated that at 2.1% cyclic-strain level 
no significant difference in the number of cycles to 
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Fig. 5 Comparison of the stress-rupture behavior of 304 stainless-steel cladding (10 to 15% cold 


work) in static sodium and helium. *? 


A, sodium tests. @, helium tests. 
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Fig. 6 Stress dependence of the average diametral creep rate of 304 stainless steel (10 to 15% cold 


work) in sodium and helium at 538°C (Ref. 59). 


failure was observed between the high-carbon sodium 
and air environments. At this same strain level, the 
cycles to failure were 34% lower than in the 
low-oxygen sodium and 72% lower than in the helium 
environment. The fatigue life of specimens tested in 
helium, air, and high-carbon sodium at 650°C, after 
exposure to a high-carbon sodium environment for 
272 hr at 650°C, decreased 13, 44, and 65%, 
respectively, compared with the as-received material. 
At all strain levels investigated, the addition of carbon 
to the sodium reduced fatigue life. The high-cycle 
fatigue properties of 304 stainless steel at 593 C were 


, sodium. @, helium. 


found to be the same in air and a sodium environment, 
based on a small number of tests.°° At present these 
data are not sufficient to assess quantitatively the 
effect of sodium on the fatigue properties of stainless 
steels. 

The low-cycle fatigue life of the 2/4 wt.% 
Cr-1 wt.% Mo ferritic steel at 593°C was the same in 
sodium and helium environments and was significantly 
greater than that in air.°° On the basis of averaged 
failures, the fatigue life in helium or sodium was 5, 10, 
and 16 times longer than in air for cyclic strains of 2.1, 
1.0, and 0.56%, respectively. 
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Fig. 8 Comparison of biaxial stress-rupture behavior of 304 
stainless steel in zirconium-gettered sodium and helium at 
760°C (Ref. 24). 
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compositional variations with emphasis on the austen- 
itic stainless steels. We are particularly interested in the 
compositional changes known to result from exposure 
of these materials to flowing sodium, i.e., carbon and 
nitrogen migration as well as selective metallic element 
mass transfer. The information will place in better 
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Fig. 10 Biaxial stress-rupture behavior of 304 stainless steel in 
flowing sodium and argon at 705°C (Ref. 22). 


perspective the property data obtained in liquid so- 
dium and will lead to further understanding of material 
behavior in this environment. 

The microstructural changes that result from ther- 
mal aging of austenitic stainless steels have been 
investigated for various heat treatments and degrees of 
cold work.” *? In these studies the amount, distribu- 
tion, morphology, and composition of various phases 
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Fig. 11 Biaxial stress-rupture behavior of 316 stainless steel in 
flowing sodium and argon at 705°C (Ref. 22). 


have been determined as a function of time and 
temperature using transmission and scanning electron 
microscopy, X-ray diffraction, metallographic tech- 
niques, and chemical analyses of extracted phases. A 
significant portion of the data is frequently reported in 
the form of time-temperature-precipitation (TTP) 
diagrams.***? Representative TTP diagrams for 316 
and 316L stainless steels are shown® in Figs. 13 and 
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Fig. 12 Fatigue behavior of 316 stainless steel in air, helium, and sodium that contains high carbon 
and low oxygen concentrations at 650°C (Ref. 27). 


14. The phases identified in these studies are generally 
denoted by o, x, n, and the various carbides, M,C, 
M,3C,, and M-C3, where M epresents the metallic 
constituent of the carbide. 

The sigma phase is primarily an Fe—Cr compound 
with 40 to 48 wt.% chromium that also contains small 
amounts of molybdenum and nickel.** The chi phase 
has been identified as a complex metallic compound 
(Fe,Ni)3, Cr;sMo, (Refs. 85 and 86) in 316 stainless 
steel. The eta phase is a Laves phase that has a 
composition of Fe,Mo and contains minor amounts of 
chromium and nickel.** In the Fe-C and Cr-C systems, 
the binary carbides Fe3C, Cr,3C,, CrjC3, and Cr,C, have 
been reported. As a result of the intersolubility of the 
binary carbides in the Fe-Cr-C_ system, the 
(Cr,Fe),,C, with up to 30 wt.% iron and (Cr,Fe),C; 
with up to 60 wt.% iron are formed.*’ These carbides 
have been observed in austenitic stainless steels, de- 
pending upon the temperature and carbon concentra- 
tion in the material.°*?! 

Figures 13 and 14 show that a small change in the 
carbon concentration in the alloy has a profound effect 
on the TTP curves and the relative amount of the 
various phases present at any given time and tempera- 
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ture. A schematic representation of the stability of the 
carbides and intermetallic phases in 316 stainless steel 
as a function of the carbon concentration in the matrix 
and temperature is shown in Fig. 15 (Ref. 83). Nitro- 
gen additions to 0.2 wt.% have been shown to have a 
similar effect in decreasing the stability of intermetallic 
phases, such as sigma, in austenitic stainless steels.” 
The state of the material with regard to the degree of 
cold work is also important from a comparison of 
Figs. 14 and 16 (Ref. 83) for 316L stainless steel. 
Similar observations have been made for the 304 
stainless steel.*' Even though these phases have been 
characterized with respect to composition, crystal 
structure, distribution, and orientation in the matrix, 
their relative importance on property changes has not 
yet been clearly established. 

The microstructural changes that result from ther- 
mal aging of nickel-base alloys have been recently 
reviewed by Sabol and Stickler.?* Phase instabilities in 
unstabilized and stabilized ferritic steels with low 
chromium and molybdenum contents have been re- 
ported by a number of investigators.°* 

Extensive data are available for the tensile, creep, 
and creep-rupture properties of 304 and 316 stainless 
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Fig. 13 Time—temperature—precipitation (TTP) diagram for 316 stainless-steel solution treated at 
1260°C for 1.5 hr and water quenched.® * 
steels in air and inert-gas environments”? °° upon nese, vanadium, silicon, and cobalt) concentrations on 
which the design codes are based. Limited data for the the mechanical properties of austenitic stainless steels 
fatigue properties in these environments are also has been investigated.''*'!? These substitutional ele- 
2 Or - 2 ss a ‘ . 4 ~ P 
available’? ''! and are currently being incorporated ments, in general, can be classified into ferrite-forming 
into design practice. Almost no data have been elements (e.g., chromium, molybdenum, vanadium, 
generated on the creep-—fatigue interaction behavior of and silicon) and austenite stabilizing elements (e.g., 
: i : ca ; : a 1 118, 119 
LMFBR materials:''*''? however, this information is manganese, nickel, and cobalt). Nakagawa et al. 
now being obtained in an air environment.'' reported that the tensile properties of Cr-Ni austenitic 
The effect of variations in the substitutional stainless steels improved slightly upon increasing the 
element (e.g., chromium, nickel, molybdenum, manga- concentrations of the ferrite-forming elements, whereas 
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increasing the concentrations of the austenite stabi- steels has been observed because of variations in the 
lizers had virtually no effect. Microstructural examina- substitutional element concentration.!'%'!? For ex- 
tion of the alloys revealed that the sigma phase formed ample, modification of 316 stainless steel by decreasing 
quite rapidly in the former at temperatures between the chromium concentration from 17 to 14% and 
500 and 800°C, whereas in the latter, delta ferrite varying the molybdenum concentration in the range of 
transformed to sigma phase only after prolonged aging. 1 to 3% resulted in slightly better stress-rupture 
The influence of a number of such elements on the behavior at temperatures''’ between 590 and 815°C. 
0.2% proof stress is shown in Fig. 17 (Ref. 115). At present, the influence of minor variations in the 
As in the case of the tensile properties, little change substitutional element concentrations on the fatigue 
in the creep-rupture behavior of austenitic stainless properties of austenitic stainless steels has not been 
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Fig. 15 Schematic representation of stability ranges of carbides and intermetallic phases in 316 


stainless steel, as a function of matrix composition.® * 


determined; however, it may be possible to predict 
whether significant changes can be expected from a 
knowledge of the elevated-temperature tensile and 
stress-rupture data for the alloys using the approach of 
Halford and Manson. '”° 

In contrast to substitutional element effects, rela- 
tively small changes in the concentration of nonmetal- 
lic elements such as carbon, nitrogen, and boron have 
been found to have a profound influence on the 
microstructure and the mechanical behavior of austen- 
itic stainless steels.‘ '°''!!2"'° The tensile and 
creep properties of 304 and 316 stainless steels have 
been investigated over a limited concentration range of 
these elements in air and inert-gas environments. The 
general conclusion from these studies is that increases 
in the concentration of carbon and/or nitrogen over 
the ranges of 0.026 to 0.094 and 0 to 0.2 wt.%, 
respectively, result in increases in the tensile strength 
and the 0.2 or 1.0% proof stress, with corresponding 
decreases in the tensile elongation and reduction in 


area at any temperature up to 700°C (Ref. 121). 
Investigations of the stress-rupture properties of 
austenitic stainless steels have indicated that nitrogen 
was essentially equivalent to carbon in its contribution 
to the high-temperature strength.'*'** Although bo- 
ron is not very soluble in austenite, it produces a 
solid-solution-hardening effect and was found to con- 
tribute to an improvement in the stress-rupture 
strength.'”' The strengthening effect of boron and 
nitrogen cannot be realized in stainless steels used in 
the core of nuclear reactors because these elements are 
set at their lowest possible concentration limits to 
minimize helium generation by means of the 1,a 
reaction and subsequent embrittlement of the material. 
The combined effect of carbon and nitrogen on the 
100,000-hr rupture strength of 304 stainless steel is 
shown in Fig. 18 (Refs. 122 and 123). The stress- 
rupture data of a number of investigators”” '°"'°*'?! 
can be used to make a similar plot for 316 stainless 
steel for a narrower composition range. 
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Fig. 16 Time—temperature—precipitation (TTP) diagram for 316L stainless-steel solution treated at 
1260°C for 1.5 hr, water quenched, and cold worked 20% prior to aging.®* 





From these studies it is evident that an increase in 
the concentration of carbon and nitrogen in austenitic 
stainless steels results in an increase in the stress- 
rupture strength while the elongation at rupture 
decreases. It should be noted that no property mea- 
surements on stainless steels with <0.03 wt.% carbon 
and nitrogen have been made in air or inert-gas 
environments; however, such data would be useful 
because stainless steels in the high-temperature sections 
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of sodium heat-transport systems attain levels in this 
range. Similarly, information on the low-cycle fatigue 
and creep-fatigue interaction behavior of these mate- 
tials is also required. 

In addition to compositional effects, thermal- 
mechanical treatments influence the microstructure 
and mechanical properties of austenitic. stainless 
steels, 1° 103 105,127,128 The stress-rupture strength of 
304 and 316 stainless steels increases as the solution- 
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Fig. 18 Influence of carbon and nitrogen on the 100,000-hr rupture strength of 304 stainless steel in 


air at 650°C (Ref. 122). 


2 : 10¢ Of: 
treatment temperature increases ” to ~1060°C and 


é : ae a + 103 
remains constant — thereafter.'°°'°> Garofalo et ai.’ 
observed that matrix precipitation of carbides was 
enhanced by cold work in 316 stainless steel. The creep 


rate decreased by a factor of 250, and the rupture life 


, laboratory heats. ¢, commercial material. 


increased by a factor of 10 at 705°C. Furthermore, the 
rupture life and ductility of this steel at elevated 
temperatures were improved at high strain rates when 
the grain-boundary carbide precipitates were allowed 
to coarsen by preheat treatment.'?” For comparison 
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with data presented in a later section, representative 
creep-rupture data” are shown in Figs. 19 and 20 for 
304 stainless steel and in Figs. 21 and 22 for 316 
stainless steel. 

Driver'*® has shown that M,3C, carbide at the grain 
boundaries in 304 stainless steel appreciably improved 
the fatigue strength at 700 and 800°C and changed the 
failure mode from intergranular to transgranular at 
700°C. This result is related to the inhibition of 
grain-boundary sliding by precipitates, according to 
Driver’s interpretation. 

Although mechanical-property data for stabilized 
austenitic stainless steels (e.g., 321 and 347) in a 
sodium environment are meager, a large amount of 
information has been published on the effect of 
thermal treatment,?”'°° '°>!?%!3° compositional varia- 
tions in carben, nitrogen, and the stabilizing elements 
such as titanium and niobium,'!?' and micro- 
structure *!29!3'!33 on the mechanical behavior in air 
and inert-gas environments. Mechanical-property data 
for stabilized®*'** and unstabilizedS* 9°19" ferritic 
steels with low chromium and molybdenum contents 
and nickel-base alloys''*'°*"'*" are also available. 
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Fig. 19 Sherby—Dorn master curve of stress-rupture data for 

304 stainless steel constructed from information in Ref. 104. 
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Fig. 20 Long-term 1.0% creep and rupture strengths of 304 
stainless steel, as a function of temperature.’ °* 


The elevated-temperature mechanical-property data 
cited in previous references are primarily for well- 
characterized materials in terms of thermal-mechanical 
history, composition, microstructure, and thermal 
aging conditions. Welding of these materials introduces 
compositional and microstructural changes and post- 
weld stress conditions that can alter the mechanical 
behavior of the material in the heat-affected zone. The 
variations in welding procedures, filler metal, and 
postweld heat treatments are additional variables that 
must be considered. Information is available on the 
mechanical properties of weldments for the austenitic 
stainless'°%'47"148 and ferritic'***7'5' steels and 
nickel-base alloys'®*'** for comparison with the 
properties of the base materials; however, a detailed 
review of this area is not within the scope of this 
article. The effect of environment on the properties of 
weldments, particularly in aqueous media where stress- 
corrosion cracking occurs, would also have to be 
considered. This phenomenon does not occur in 
materials exposed to a liquid-sodium environment. 
However, stress-corrosion cracking in a caustic environ- 
ment is of fundamental importance when considering 
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Fig. 21 Sherby—Dorn master curve of stress-rupture data for 316 stainless steel.° ° 


the removal of hardware from a reactor, removing the 
sodium by a water dissolution process, repairing the 
component, and placing it back in service. 


Discussion of Mechanical-Property Data in Sodium, 
Air, and Inert-Gas Environments 


The compositional and microstructural changes 
that result from the exposure of LMFBR materials to 
high-temperature sodium were described in a previous 
subsection, as were the mechanical-property data for 
these materials. The effects of compositional variations 


on the microstructure and properties for a number of 


these alloys have been determined in alloy develop- 


ment studies that were not related to the sodium 
environment. The premise in this subsection is that the 
mechanical behavior of materials (cui exceptione, 
fatigue) tested in sodium and gas environments is the 
same, to good approximation, provided the comparison 
is based upon composition, thermal-mechanical treat- 
ment, microstructure, and test parameters (e.g., tem- 
perature and stress conditions). At present, a semi- 
quantitative comparison of this type can be made for 
the creep-rupture properties of austenitic stainless 
steels. 

As shown in the preceding subsection, the influ- 
ence of small variations in the substitutional element 


concentrations (e.g., chromium, nickel, manganese, 
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Fig. 22 Long-term 1.0% creep and rupture strengths of 316 
stainless steel, as a function of temperature.° ° 


and molybdenum) on the tensile and creep-rupture 
properties of austenitic stainless steels was minimal. 
For this reason, and also because the steady-state 
thickness of the modified surface region is small in 
comparison with the section thickness, the effect of 
selective mass transfer of these elements in sodium on 
the mechanical properties is likely to be small in many 
instances. In the case of fuel cladding near the 
end-of-life, corrosion and possible wear and fretting by 
the helical wire wrap or other types of fuel-element 
spacer and support systems must be considered in the 
design allowance. 


Variations in the nonmetallic element concentra- 
tions in austenitic stainless steels, particularly carbon 
and nitrogen, were shown to have a profound effect on 
microstructure and mechanical behavior in air and 
inert-gas environments. We now show the correlation 
between these results and the data obtained in a 
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sodium environment from both a microstructural and 
property viewpoint. However, before proceeding, it is 
important to note that when no compositional and 
related microstructural changes occurred, the tensile 
and creep-rupture properties determined in the three 
environments (i.e., air, inert gas, and sodium) were 
essentially the same. This indicates that sodium per se 
does not alter the mode or rate of fracture under 
tensile and creep-rupture testing conditions. The micro- 
structural changes that result when 304 and 316 
austenitic stainless steels undergo decarburization in 
sodium are very similar to those observed when low- 
carbon heats of these steels are exposed to the same 
thermal treatment. The most notable characteristic is 
the enhanced formation of sigma phase; i.e., the 
amount of sigma phase at any time is significantly 
higher compared with that in materials with initial 
carbon concentrations of 0.05 to 0.07 wt.%. A quanti- 
tative comparison of the amount of sigma phase in the 
alloys decarburized in sodium with that in the low- 
carbon stainless steels for similar aging conditions 
cannot be made because the volume fractions of the 
various phases were not determined in the materials 
exposed to sodium. Because other intermetallic phases 
are involved in the thermal aging process, property 
changes cannot be interpreted solely in terms of sigma 
phase or carbide formation. 

The 650°C creep strength (10,000-hr rupture 
stress) of annealed 304 and 316 stainless steels is 
plotted as a function of the final carbon and nitrogen 
contents of the alloys in Figs. 23 and 24, respectively. 
The rupture data for 304 stainless steel, represented by 
the circles in Fig. 24, were obtained by interpolating 
the 1000- and 100,000-hr values reported by Goodell 
et al.'*? The dashed lines represent the uncertainty in 
the stress-rupture values. The data in sodium and the 
corresponding control tests in inert-gas environments 
were obtained by extrapolating the results of Spalaris 
et al.2”?3 and Lee”* to 10,000 hr. Inasmuch as the 
latter results were determined at temperatures between 
538 and 760°C, the values were normalized to 650°C 
based upon the temperature dependence of the 
10,000-hr creep strength of 304 stainless steel given 
earlier in Fig. 20. Information on the carbon and 
nitrogen content of the steels from the control tests in 
inert-gas environments and after exposure to sodium is 
required for comparison with Goodell’s results in 
Fig. 18, in which the composition parameter is 
(C + 1.25 N) wt.%. However, the initial and final nitro- 
gen contents of the alloys were not reported in the 
work of Spalaris.”” Therefore, for the purpose of this 
comparison, a typical nitrogen content for the alloys 
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(the same as that reported by Lee)”* was assumed, and 
the percentage decrease in the nitrogen concentration 
in the alloys exposed to sodium was the same as that 
observed for carbon. The observation that these mate- 
rials lose carbon and nitrogen at similar rates in flowing 
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sodium at high temperatures’** forms the basis for the 
above assumption. 

The same procedure was used to construct Fig. 24 
for 316 stainless steel, where the out-of-sodium data 
were obtained from the work of Hopkin and Taylor,’ 
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Mazza and Willoughby, '®? and Garofalo et al.'°? The 
uncertainty bands are the same magnitude as those for 
304 stainless steel shown in Figs. 18 and 23. 

In the case of 304 stainless steel, the stress-rupture 
data obtained in sodium and the control tests in 
inert-gas environments (Fig. 23) fall below the un- 
certainty band given inGoodell’s work. This difference 
may partly be due to variations in the solution 
treatment of the specimens, i.e., half-hour at 1035 and 
1120°C in Lee’s”* and the work of Goodell et al.,!* 
respectively, and a “mill anneal” in the Spalaris and 
Zebroski investigation.””?* The effect of the solution- 
treatment temperature on the stress-rupture strength of 
304 and 316 stainless steels has been investigated by 
White and Freeman.’ An increase of 40 to 65°C in 
the solution-treatment temperature can result in a 
1000- to 1500-psi increase in the rupture strength of 
these alloys. 

The results in Figs. 23 and 24 indicate that the 
10,000-hr stress-rupture strengths are decreased by 
approximately 25 to 35% as a result of carbon and 
nitrogen depletion upon exposure to high-temperature 
sodium. Lower carbon and nitrogen concentrations in 
these materials have been observed after exposure to 
flowing sodium.'® The stress-rupture behavior of alloys 
with nonmetallic element concentrations <0.04 wt.% 
has not been evaluated in either sodium or gas 
environments. 

A further observation can be made concerning the 
difference in the slope of the stress vs. time-to-failure 
data in sodium relative to the control tests in argon 
(see Figs. 10 and 11). Because the test specimens were 
undergoing decarburization during the period of the 
creep tests, the amount of carbon (and nitrogen) lost 
from the specimens increased with the test duration. 
The time required to reduce the carbon concentration 
in the 15-mil- wall 364 stainless steel fuel-cladding 
tubing from the initial value of 0.062 wt.% to the final 
value of 0.017 wt.% at 705°C is approximately 2000 to 
3000 hr. Therefore the specimens tested at the higher 
stress levels had undergone a relatively small amount of 
decarburization, and the rupture times are only slightly 
less than those for the control specimens. At the lower 
stress levels, and correspondingly longer rupture times, 
the specimens were at the lower (or equilibrium) 
carbon level for a large fraction of the test duration. 
Consequently, only the data obtained at times consid- 
erably greater than the equilibration period should be 
considered when extrapolating the results to longer 
times. 

The decrease in the diametral strain during de- 
carburization of the stainless steels in sodium? is in 
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contrast to the increase in ductility with a decrease in 
carbon concentration observed out of sodium.'?? Other 
results showed almost no change in the rupture 
ductility when decarburization in sodium occurred.***! 
The ductility at rupture in all cases was on the order of 
5%. Although we cannot resolve the differences ob- 
served in the various investigations, Kramer'®® has 
reported that different test environments (e.g., Ar, N>, 
O,, H,, air, and vacuum) can exert an influence on the 
creep rate and stress-rupture life of materials such as 
304 stainless steel, an 80 wt.% Ni-20 wt.% Cr alloy, 
nickel, and aluminum. 

The tensile and low-cycle fatigue properties of 304 
and 316 stainless steels have not been evaluated for low 
nonmetallic element concentrations in the materials. In 
addition to the compositional and microstructural 
effects, surface condition and environment can influ- 
ence the fatigue behavior. The relative contributions of 
these variables cannot be determined from the limited 
information obtained in a sodium environment. Several 
investigators have shown that the fatigue life is 
significantly better in vacuum than in air for 0.5% 
molybdenum steel,'®’ several nickel-base alloys,'** '® 
316 stainless steel,’ and a 26 wt.% Ni-15 wt.% 
Cr-2 wt.% Ti steel.'®' It was also found that fatigue 
cracks produced by slow cycling at elevated tempera- 
tures were nearly always filled with oxide prod- 
ucts.'°7!© In the latter material a highly localized 
surface oxidation was identified as the nucleation site 
for fatigue cracks.'©" 

Inasmuch as oxidation products of this type are 
not thermodynamically stable at high temperatures in 
sodium that contains low concentrations of oxygen, 
which is typical of LMFBR systems, fatigue data 
obtained in air may be conservative for establishing 
design criteria. The small amount of fatigue informa- 
tion for the 2% wt.% Cr-1 wt.% Mo ferritic steel in 
sodium, helium, and air environments®° strongly sug- 
gests that the data in air are conservative. Therefore 
parallel investigations of the fatigue and creep-fatigue 
interaction behavior of all these materials in well- 
characterized sodium are necessary to minimize un- 
certainties in design data. 


IRRADIATION EFFECTS ON 
MECHANICAL PROPERTIES 


The effect of fast-neutron irradiation on the 
mechanical behavior of reactor materials is briefly 
reviewed for comparison with results from exposure to 
the sodium environment. Possible interaction between 
the sodium and fast-neutron environments leading to 
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synergistic effects on the mechanical-property behavior 
is explored. 

The tensile properties of unstabilized austenitic 
stainless steels have been studied extensively by post- 
irradiation tests. The properties deteriorate signifi- 
cantly after exposure to fast-neutron irradiation: how- 
ever, indications are that the property changes may 
saturate after some intermediate level of exposure. A 
number of investigators'®* '”' have shown that both the 
yield and tensile strengths increase and the ductility 
decreases as fluence increases. Typical data'® which 
illustrate these changes are shown in Fig. 25. The 
general conclusions from these studies are: 

1. At low fluence levels [about 10°° neutrons/ em? 
(E >1MeV)], displacement damage results in radia- 
tion hardening if the irradiation and postirradiation 
test temperatures are below one-half the absolute 
melting temperature (0.5 7,,,) of the steel. 

2. Above 0.5 7);,,, hardening no longer exists, but a 
ductility loss due to helium embrittlement is observed. 

3. The irradiation temperature above which radia- 
tion hardening is no longer significant increases as 
fluence increases. 


Postirradiation creep-rupture properties of un- 
stabilized austenitic stainless steels have been studied 
by a number of investigators.'°%'°*!7'* 
conclusions from these studies are: 


General 


1. Strain at failure depends on irradiation and 
postirradiation test temperatures. 

2. The greatest ductility losses occur at the lowest 
irradiation temperature where displacement damage is 
the major contributor. 

3. For temperatures below 580°C, the displace- 
ment damage controls creep ductility; between 580 
and 650°C, ductility losses are due to the combined 
effects of helium production and displacement damage; 
above 650°C, only helium embrittlement is controlling. 

4. The stress-to-rupture (1000- or 10,000-hr rup- 
ture stress) in postirradiation tests decreases with an 
increase in fluence and attains about one-half the 
unirradiated value near 10** neutrons/cm?, but little 
change occurs with higher fluence. 


Typical tube creep-to-rupture data'®' for 316 stain- 
less steel are given in Fig. 26. The dependence of 
rupture life and elongation at rupture on fast-neutron 
fluence is shown in Figs. 27 and 28, respectively,'”* for 
304 stainless steel. It is interesting to note in Fig. 29 
that no strong difference is evident in the stress-rupture 
behavior of 304 stainless-steel base or weld metal in 
either the irradiated (7 x 107! neutrons/ cm?) or un- 
irradiated conditions'” at 593°C. 
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Fig. 25 Effect of fast-neutron irradiation on the tensile 
properties on 304 stainless steel irradiated at 371 to 460°C and 
tested at 550°C (Ref. 162). 


Limited data are available for the effect of irradia- 
tion on the low-cycle fatigue behavior of austenitic 
stainless steels. At low fluences the available data' 
(Fig. 30) show that fatigue life is reduced by irradia- 
tion. 

The effects of fast-neutron irradiation on the 
swelling behavior and mechanical properties of stabi- 
lized austenitic stainless steels'8*'8* and _nickel-base 
alloys'**'87'®? have not been evaluated to the same 
extent as the unstabilized austenitic stainless steels. In 
the case of several stabilized austenitic stainless steels, 
the influence of alloying additions such as niobium and 


84 


: P . ee 1 ° 
titanium on the swelling characteristics“ and tensile 


5 


properties’ and their dependence on the irradiation 


temperature have been investigated. Collins et al.'”° 
have reviewed the comparative swelling and property 
behavior of a number of these materials. Interactions 
between the stabilizing additions and nonmetallic 
elements, e.g., carbon, nitrogen, and boron, in the 
alloys that result in precipitation of carbides, nitrides, 
and borides, of different morphology and their relation 
to helium accumulation at grain boundaries and precip- 
itate particles are considered to be important in this 
regard. 

From the above and discussion in the previous 
subsections, both irradiation and sodium environments 
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Fig. 26 Biaxial stress-rupture behavior of 316 stainless steel after fast-neutron irradiation at 480°C 
(Ref, 181). 
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Fig. 29 Effect of fast-neutron irradiation on mechanical properties of austenitic stainless-steel 
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Fig. 30 Fatigue behavior of 304 and 304L stainless steels at 500°C after fast-neutron irradiation at 
450 and 750°C (Ref. 183). 
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influence the microstructure and mechanical behavior 
of LMFBR materials. With the accumulation of ap- 
preciable fluence, the effect of void, dislocation loop, 
or helium bubble formation on the mechanical be- 
havior usually overwhelms the property changes that 
result from sodium exposure. However, the develop- 
ment of irradiated microstructures can be influenced 
significantly by mass transport of nonmetallic elements 
in sodium, i.e., loss of carbon and nitrogen from 
austenitic stainless steels in high-temperature sodium 
results in dissolution of carbides and nitrides and 
promotes sigma-phase formation. Second-phase parti- 
cles can influence the size and distribution of helium 
bubbles at high temperatures, which have an important 
role in the mechanical behavior and the fracture mode. 

It has not been established whether the micro- 
structural changes that result from mass-transport 
processes in sodium environment can be influenced by 
irradiation. 


ADHESION, WEAR, AND FRICTION 
BEHAVIOR OF MATERIALS 
IN A SODIUM ENVIRONMENT 


Materials used for a variety of mechanisms that are 
required to operate in a high-temperature sodium 
environment must possess adequate corrosion and wear 
resistance as well as antigalling and low-friction char- 
acteristics. For components such as fuel-element- 
transfer devices, values, shaft seals, subassembly wrap- 
per cans, or control-rod mechanisms, the designers are 
faced with a common problem of finding pairs of 
materials, which, under the desired operating con- 
ditions, are not corroded by sodium, do not react with 
each other, do not adhere to each other, and can slide 
over each other with low friction when lubricated only 
by liquid sodium or its vapor. The three phenomena of 
importance in a sodium environment are (1) adhesion 
or self-welding, (2) wear, and (3) friction. 


Adhesion 


Adhesion or self-welding is a phenomenon that 
generally occurs at the contact point of two metals or 
alloys. The extent to which adhesion can occur 
depends primarily on the material couple, the surface 
conditions, and the environment; and to a lesser extent 
on time and temperature.'?' The absence of protective 
oxides and the presence of exceedingly clean surfaces 
generally promote adhesion. On the basis of adhesion 
theory advanced by Bowden and Tabor,'®!' the magni- 
tude of adhesion is proportional to the load normai to 
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the contact surfaces and is independent of the contact 
area. 

In a sodium environment, potential adhesion prob- 
lems can occur in valve seats, various mechanisms in 
control-rod drives, between cladding and its spacer-wire 
wrap, and at interfaces between subassemblies. Very 
localized adhesion or self-welding of 316 stainless steel 
fuel-cladding tubing to spacer wires as well as adhesion 
of loose wires to neighboring pins has been reported by 
Shively.'°? An evaluation of numerous material couples 
as to their adhesive tendencies in high-temperature (to 
815°C) sodium has been made by a number of 
investigators,'°* '> and results of this type have been 
summarized '*'%’ in the literature. 


Wear 


Sliding contact in liquid sodium usually induces 
adhesive wear, as described by Rabinowicz.'*® The 
volume of material removed by the wear process is 
given by V=CL x, where C is a coefficient that 
depends on the materials, L is the contact load, and x 
is the sliding distance. It is frequently observed that 
C=k/p, where p is the hardness of the softer material 
of the sliding configuration and k is another coefficient 
not related to hardness. The effect of a protective film 
on the material or the presence of a sodium environ- 
ment, as compared with helium or argon, is to decrease 
the coefficient k. Increasing the hardness of a material 
couple would also lessen the wear. With an accurate 
knowledge of k, it should be possible to predict wear 
rates for a given sliding system. 


Wear coefficients have been determined experi- 
mentally for a number of different materials. Values of 
k for the softer member of a sliding couple under 
continuous unidirectional sliding are being reported in 
the literature. The k values given by Rabinowicz'” for 
unlubricated surfaces vary from 10’ to 10‘, where the 
median for a majority of the couples falls between 10° 
and !07. It is anticipated that, for a sliding system 
immersed in sodium with low concentrations of non- 
metallic elements (e.g., oxygen, carbon, and nitrogen), 
surface chemical effects on the material couple would 
have a strong influence on the k values. 

A number of investigators have measured wear 
rates in sodium and gas phases for a variety of material 
couples such as tool steels,'°”? '9 molybde- 


copper-, 
Stirred ‘tai  — 
alloys, austenitic steels, 


base 199200 202-204 199, 200, 202-204 
cemented carbides,'"°* °° precipitation hardened 
199 199, 206-208 a : 
steels, “* and pure metals. Representative wear 
data for a number of material couples are included 








EFFECT OF SODIUM ON LMFBR MATERIALS 


with friction results in the next subsection. A con- 
siderable amount of information of this type has been 
developed by component manufacturers and is propri- 
etary and cannot be referenced. 


Friction 


Friction is the resistance to motion that exists 
when a solid object is moved tangentially with respect 
to the surface of another that it touches, or when an 
attempt is made to produce such motion. The coeffi- 
cient of friction, defined in terms of the bulk proper- 
ties of the materials, is the ratio of the shear strength 
(or cleavage strength of the bond of weaker material in 
the couple) to the yield pressure of the softer material 
(indentation hardness).”° Therefore, by this defini- 
tion, the friction coefficient should be independent of 
load, apparent area of contact, and sliding velocity, but 
is usually influenced strongly by the environment and 
surface finish of the materials. This has been found to 
be true only within limited ranges for sliding velocities, 
loads, and surface areas. Since high-purity sodium 
tends to clean metal surfaces of native films, very high 
friction coefficients are expected. Values as high as 1.0 
have been obtained from friction measurements for 
alloys such as Inconel 718 on 304 stainless steel, and 
3.8 for stainless steel on stainless steel.?'° 

The friction coefficient is not closely related to 
wear. For example, Barrau etal.” reported that 
nitrided steel against nitrided steel gave friction coeffi- 
cients less than 0.40 and a large total wear; however, 
WC-8 wt.% Co against WC-6 wt.% Co gave no detect- 
able wear, but the friction coefficients were as large as 
0.8. 


A number of investigators 
evaluated friction coefficients for a variety of material 
combinations in a sodium environment, and these 
results have been summarized by Baumgartner,'®’ 
Balkwill,?’? and Hoffman et al.'*? In general, stellites, 
tungsten carbides, and titanium carbides vs. each other 
and in certain combinations yield the lowest coeffi- 
cients of friction (0.3 to 0.6) in sodium. Typical 
information for the mean friction coefficient and 
specific wear rates for several of these alloys at 
temperatures between 200 and 500°C in liquid sodium 
and gaseous environments with and without sodium 
vapor are shown 7” in Figs. 31 to 34. 


199-203,207,210,211 rave 
.< 


Depending on the material combinations, tempera- 
tures, and sodium purity, namely, oxygen concentra- 
tion, the formation of binary or ternary oxide films on 
the materials may or may not occur. Free energies of 
formation of binary and ternary oxides of a number of 
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Fig. 31 Wear rates and friction coefficients for Nimonic 75 
against spray-fused Stellite 6, in sodium, carbon dioxide, and 


200 


helium, as a function of temperature. * 


metals of primary interest are shown in Fig. 35 along 
with free-energy data for sodium oxide and sodium 
that contain 1, 5, 25, and 50 ppM oxygen by weight. 

In EBR-II, which may or may not be prototypic of 
future LMFBRs from the standpoint of sodium chem- 
istry, the oxygen concentration in sodium is main- 
tained at a level of ~0.8 ppM,7'? and, of the materials 
in Fig.35, only the oxidation of chromium and 
titanium to form a complex (NaCrO;) compound and 
TiO, respectively, can be expected at oxygen levels of 
~1 ppM. At temperatures above 475°C in low-oxygen 
sodium, even NaCrO, is unstable, and thus various 
mechanisms would operate with essentially clean sur- 
faces. The most extensive data on the friction and wear 
behavior of material combinations in low-oxygen so- 
dium of controlled purity at high temperatures have 
been obtained by Hoffman et al.'” 

Glaeser?'* has recently made a detailed survey and 
analysis of material wear and friction in sodium and 
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also’: summarized the estimated operating conditions 
and failure criteria for a number of Fast Test Reactor 
(FTR) core components. It was concluded that data 
from conventional wear tests could not be applied to 
the detailed design requirements of Fast Flux Test 
Facility (FFTF) core and in-vessel fuel-handling ma- 
chine components. The complexity of wear and fric- 
tion problems requires that simulation tests be per- 
formed under prototypic conditions for a specific 
component (e.g., contact geometry, type of motion, 
area ratios, loading conditions, sodium purity, and 
temperature) prior to finalizing the design of a reactor 
component. 
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Fig.32 Wear rates and friction coefficients for various 
Stellites against Stellite 1, in sodium, argon, and in argon that 
contains sodium vapor, as a function of temperature.’ °° 


CONCLUSIONS 


An extensive amount of mechanical-property infor- 
mation is required for the design of LMFBR compo- 
nents. Most of the available data have been generated 
in air and inert-gas environments, and the relatively 
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Fig. 33 Wear rates and friction coefficients for a nitrided low 
Cr—Mo steel against itself, in sodium, carbon dioxide, and 
helium, as a function of temperature.” °° 


limited data in either the sodium or irradiation 
environment are used to establish the necessary safety 
factors. In many instances the property data obtained 
in air do not encompass the entire range of composi- 
tional variations (namely, carbon and nitrogen) in the 
base materials that are observed upon exposure to a 
high-temperature sodium environment. Test programs 
are in progress in which materials from special heats of 
controlled composition are being used along with 
materials preexposed to sodium to develop a range of 
compositional and microstructural variations in the 
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materials prior to testing. The result when correlated 
with in-sodium property data for the same tempera- 
ture, sodium purity, and test parameters will aid 
considerably in reducing the uncertainties in the safety 
factors for property data used in the design of 
out-of-core components. 

The extent to which these compositional variations 
affect the microstructures, swelling behavior, and 
mechanical properties of in-core materials has not been 
established mainly because of the low coolant tempera- 
ture presently available in EBR-II. Future experiments 
with bypass flow subassemblies in EBR-II will result in 
higher coolant and fuel-cladding temperatures and 
therefore will yield some data for the combined effect 
of high-temperature sodium and irradiation on the 
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Fig. 35 Thermodynamic stability of oxides for a number of 
metals that have application in LMFBR sodium systems 
relative to sodium which contains a range of oxygen concentra- 
tions. 


properties of in-core materials. Experiments in the 
FFTF will provide most of the information that results 
from synergistic effects on the swelling and mechanical 
behavior of these materials. Because the magnitude and 
rate of carbon and nitrogen migration in sodium 
heat-transport system materials in EBR-II and FFTF 
will be different, mechanical-property data obtained in 
these facilities may require interpolation or extrapola- 
tion when applied to other reactor systems. A mod- 
eling approach based on both in-core and out-of-core 
experimental results, where a wider range of variables 
can be explored, will be necessary to establish reliable 
design and performance limits for materials in com- 
mercial LMFBRs. 

Much of the data on the wear and friction behavior 
of materials in a sodium environment has been ob- 
tained under inadequately defined test conditions, due 
to the lack of capability to measure and control the 
concentration of nonmetallic impurities in sodium. 
Recent developments in sodium technology have re- 
sulted in standardized sampling and analytical 
methods, tighter specifications for sodium and cover- 
gas purity, and a capability to monitor continuously 
the nonmetallic impurity levels in sodium. Incorpora- 
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tion of this stronger analytical base into wear and 
friction studies will reduce the uncertainties in predict- 
ing the behavior of material couples in the sodium and 
cover-gas environments and will aid in the selection of 
materials capable of better performance under specific 
service conditions. Better characterization and control 
of the environment will enable a more quantitative 
evaluation of the effect of contact geometries, types of 
motion, and loading conditions on the friction and 
wear behavior and further enhance the usefulness of 
these data for the design of reactor components. 
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Nucleation of Liquid Sodium 
in Fast Reactors 


By Hans K. Fauske* 


Abstract: Analytical models and experimental studies pertain- 
ing to superheating of liquid sodium including both homoge- 
neous and heterogeneous nucleation are reviewed. Their 
applicability to analysis of hypothetical fast reactor accidents 
including channel voiding in relatively intact core geometry, 
as well as violent boiling due to molten fuel—coolant inter- 
action, is discussed. 


In the current design of sodium-cooled fast breeder 
reactors, a coolant temperature at least 200 to 300°C 
below the saturation temperature is anticipated during 
normal operation. Despite this large safety margin and 
the fact that current and future designs of liquid- 
metal-cooled fast breeder reactors (LMFBRs) will 
assure that a loss-of-cooling incident is an event of 
negligible probability, reactor safety analysis must 
consider coolant-temperature excursions to the extent 
that boiling may occur. Two rather distinct boiling 
processes are of major interest in the analysis of 
abnormal reactor behavior. The first, termed channel 
voiding, concerns sodium boiling and motion within 
relatively intact reactor-coolant channels. The second, 
termed violent boiling, concerns vapor production and 
pressure generation as collapsing fuel elements or 
molten fuel-element material contact sodium. The 
primary interest in channel voiding results from the 
need to assess rates of sodium-voiding reactivity feed- 
back. The primary interest in violent boiling results 
from the need to assess the magnitude of potentially 
destructive pressure pulses. A basic understanding of 
nucleation and superheating of liquid sodium is needed 
to develop an analytic, predictive capability of these 
processes. 





*Reactor Analysis and Safety Division, Argonne National 
Laboratory, Argonne, IIl. 60439. 
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In the case of channel voiding, the nucleation may 
be initiated from microcavities present on the heating 
surface and/or gas bubbles in the liquid. The bubble 
growth in either case is accomplished against the liquid 
pressure and the surface-tension forces by an eleva- 
tion of the pressure of the vapor. Since the vapor 
and liquid across the bubble surface must be in thermal 
equilibrium, the temperature must be higher than the 
saturation temperature of the liquid. The liquid super- 
heat is defined as Ty4)—T saturation at the initiation 
of nucleation. On the other hand, if the boiling process 
results from the contact of molten fuel and subcooled 
liquid sodium, superheating as high as that correspond- 
ing to spontaneous nucleation (explosive boiling) may 
be possible as a result of the excellent wetting 
properties of these fluids, assuming a preexisting 
gaseous phase is not present. However, in reactor 
systems using commercial-grade sodium, it may be 
argued that significantly lower superheats will be 
required for nucleation. This behavior may result from 
the large variety of surface conditions present in the 
engineering system, from gas bubbles in the sodium 
which serve as nucleation sites, or from radiation 
effects. 

Theoretically, under identical surface conditions, 
superheats for sodium systems are about an order of 
magnitude greater than for nonmetallic systems. Figure 
1 illustrates this order of magnitude difference in 
superheats. The higher results for sodium are primarily 
due to the higher boiling point, higher surface tension, 
and lower vapor density of sodium in the low-pressure 
region. However, as will be discussed, the governing 
mechanisms for sodium vs. water boiling require 
careful consideration of the effective size of the active 
wall cavity and/or bubble size in the liquid. Further- 
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more, the consequences of boiling initiation at a given 
degree of superheat in sodium are considerably differ- 
ent than for nonmetallic systems. As a result of the 
low conductivity of water, near the wall a very steep 
temperature profile characteristic of high-Prandtl- 
number (20.7) fluids is observed. Hence only a small 
portion of the nonmetallic fluid near the wall is 
superheated, whereas the rest remains subcooled. For 
liquid metals, with Prandtl number 0.004 to 0.005, the 
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Fig. 1 Superheat values for sodium and water for a given 
cavity size. 


temperature profile is much less steep at the wall so 
that a much larger fraction of the liquid can generally 
be superheated. The more shallow temperature profile 
and the small film temperature drop for sodium 
(illustrated in Fig. 2) allow the superheat energy that 
can be stored in a sodium system to exceed that of a 
water system. The events following nucleation in the 
two systems (flow regime and void growth) can, 
therefore, be expected to be quite different. In the case 
of sodium, the first bubble that nucleates will contain 
vapor at a higher pressure than the static pressure in 
the surrounding liquid because of the superheat, high 
thermal conductivity, and shallow temperature gra- 
dient. This pressure increase will tend to deactivate any 
other potential nucleation sites in the neighborhood of 
the first bubble and therefore gives rise to a single large 
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Fig. 2. Comparative energy storage in water and sodium 
reactor channels. 


bubble rather than very many small bubbles obtained 
in reactor water systems. 

The subject of liquid-metal superheating has been 
investigated in a number of laboratories, in view of its 
close connection to problems of fast reactor safety. In 
these experiments, superheats of the order of 0 to 
900°F have been observed for a variety of liquid 
metals, surface finishes, and other variables. The large 
spread in data is indicative of the fact that many 
variables may affect superheat. Figure 3, presented by 
Fauske,’ demonstrates the state of the art as of the 
International Conference on the Safety of Fast Reac- 
tors at Aix-en-Provence, France, September 1967.* 

The purpose of this article is to critically review 
existing nucleation models and experimental data 
which have appeared since 1968 as well as to offer 
some new interpretations concerning the governing 
mechanisms. 


NUCLEATION THEORIES 


The potential sources of nucleation include the 
following: 

1. Homogeneous nucleation due to the spontane- 
ous formation of critical-size nuclei containing no 
foreign matter and disallowing the possibility of 
nucleation at the container walls. 

2. Vapor-free heterogeneous nucleation due to wet- 
ted suspended microscopic particles or wetted surfaces. 

3. Heterogeneous nucleation within the liquid due 
to suspended microscopic particles that contain perma- 





“lor detailed discussion of these early data, consult Refs. | 
to 16. 
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Fig. 32 Superheat data presented at the Aix Conference. 


nent gas pockets or, more particularly, due to sus- 
pended small bubbles of inert gas. 

4. Heating at the walls with surface crevices or 
cavities containing gas and/or vapor. 


Homogeneous Nucleation 


The work of formation of a vapor nucleus of 
critical size is given by 


W = oS — (pg — Po)V (1) 


The term oS represents the work required to form the 
surface, and the term (Pg — Pg)V is the work gained in 
forming the inertia mass. The second term on the 
right-hand side of Eq. 1 is always two-thirds of the 
first. W is therefore positive, and the original phase can 
therefore be said to be stable, where W is a measure of 
its stability. Using the Laplace equation, 


20 = (Dg — Pe We (2) 


lets the value of W be expressed without involving any 
geometrical magnitudes: 
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Fig. 3b Superheat data published previous to the Aix Con- 
ference. 
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For a pure fluid we see that the stability, thus 
measured, is infinite for an infinitesimal difference of 
pressures; i.e., a pure liquid can exist in a superheated 
state. It is also clear that the stability decreases very 
rapidly as the difference of pressures increases and the 
surface tension decreases. 

The rate of nucleation, J, i.e., number of critical- 
size nuclei per unit volume per unit time which grow to 
macroscopic size, is approximately proportional to the 
product of the number of vapor embryos per unit 
volume and the frequency of collision between the 
embryo and individual liquid molecules w. Thus the 
rate of formation of embryos of a critical size is given 
by the general rate expression 


_W 
J=wN,e KT (4) 


where N, is the number of molecules per unit volume 
of liquid and K is the Boltzmann constant. Depending 


NUCLEATION OF LIQUID SODIUM 


on the limiting process in the dynamics of a vapor 
bubble, the coefficient w has different forms, retain- 
ing, however, a numerical value close to 10'° sec! 

The interesting aspect of Eq. 4 is that the superheat (or 
fracture pressure) is remarkably insensitive to the 
waiting time of the first bubble. This implies that very 
small temperature changes occur for large changes in 
the rate of nuclei formed per unit volume per unit 
time. This is illustrated in Fig. 4 for sodium. It shows 
the calculated temperature T at which bubble forma- 
tion due to spontaneous nucleation is to be expected. 
It is clear that the calculated superheat depends very 
slightly on the rate of nuclei formation. Excellent 
experimental verification for the above theory has been 
obtained in some special cases. However, in most 
engineering systems experiencing moderate heating 
transients, incipient nucleation is generally observed to 
occur considerably before the threshold corresponding 
to spontaneous nucleation is reached. Figure 4 includes 
the region of measured superheat values for liquid 
metals and that region is seen to fall considerably 
below the theoretical values; this fact suggests that 
some other mechanism is controlling the nucleation 
process. It should be noted that the data shown in 
Fig.4 include only experiments dealing with slow 
heating transients from a solid heating surface. With 
very rapid heating transients and/or where the heating 
surface is a liquid, the threshold of homogeneous 
nucleation may well be reached. For further detailed 
discussion of homogeneous-nucleation theories, see 
Refs. 17 and 18. 
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Vapor-Free Heterogeneous Nucleation 


Nucleation that is initiated by foreign nuclei 
(heterogeneous nucleation) arises from the catalytic 
effect of their surfaces, as well as from the walls of a 
vessel, grain boundaries, pores, etc. The essential 
condition is that the surface be wet. 

Volmer'? considered the case of a flat surface, S, 
upon which a phase transformation from liquid to gas 
takes place with a contact angle, 8, measured through 
the liquid phase (Fig. 5a). On the basis of a force 
balance at the triple interface, the following interfacial 
energy equation can be written as 


iis can i 5 

915 = 9G 5 — %G_1 6088 (5) 

Then the reversible work required to form an 
embryo having the shape of a spherical sector is given 


as 


16707? 





Ww =————- f(0) (6) 
3D, ag De) 
where 
anes _ wn 9 
fig) === F882 cos” 6 (7) 


For 0 = 0° (perfect wetting), when the vapor embryo 
just touches the surface, f(@)= 1, and the limiting 
equation for homogeneous nucleation is obtained. 
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Fig. 4 Homogeneous nucleation for liquid sodium based on Eq. 4. 
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Fig. 5a Embryo on flat surface. 


Conversely, if @ = 180° (no wetting), the new phase 
spreads perfectly upon the solid surface, and no 
activation energy is necessary. Since liquid metals 
generally have contact angles very much less than 90° 
at temperature levels of interest in current LMFBR 
design, heterogeneous nucleation on a flat surface in 
tne absence of a preexisting gaseous phase would 
appear not to be able to explain existing data. 

Bankoff?° extended Volmer’s nucleation theory to 
curved surfaces and reached a similar conclusion. More 
recently, Kottowski?' proposed that necked cavities 
might act as nucleation sites. For a necked cavity, such 
as in Fig. 5b, the function f(@) can be expressed in 
terms of the apex angle 6 and an equivalent contact 
angle @ on a flat surface = 7 — 8 + 6’. The geometrical 
factor in Eq. 6 is thus 


_2+3 cos (m— B+ 6’) —cos® (m7 — B+ 0’) (8) 


f0') 7 





As a result of the chosen geometry, the activation 
energy required to generate a new phase is considerably 
reduced compared to the flat surface case of Volmer, 
even if the contact angle remains small. In fact, as the 
apex angle 6 > 0, the work W- 0. Therefore, at first 
glance, it would appear so that cavities of the shape 
illustrated in Fig.5b5 would require a rather small 
superheat in order to initiate boiling even for a liquid 
having a small contact angle. However, a closer look at 
the cavity illustrated in Fig. 5b reveals that, although 
the superheat required to generate vapor may be small, 
the internal pressure must be extremely high in order 
for the bubble to escape the neck of the cavity; in fact, 
that pressure must exceed the vapor pressure corre- 
sponding to the superheat temperature required for 
homogeneous nucleation, since the neck radius is 
indeed less than r,, as calculated from Eqs. 2 to 4. This 
is obviously not possible, since bulk nucleation would 
take place first. Therefore well-wetted cavities can be 
rejected as nucleation sites. 
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Fig. 5b Embryo within a necked cavity. 


Heterogeneous Nucleation; Bubbles in Liquid 


Because of the solubility characteristics of inert 
gases in liquid metals, gas bubbles are likely to be 
present in the coolant in most convective systems, 
including reactor systems. The presence of gas bubbles 
could be due to either gas precipitating from solution 
in the heat exchanger of the system (decreasing 
temperature) and/or from mechanical entrainment at 
the free surface in the system. The strong influence of 
inert-gas concentration in the circulating liquid can be 
illustrated by first considering pure vapor bubbles. For 
curved interfaces, mechanical equilibrium requires that 
surface-tension forces balance the pressure difference 
between vapor and liquid: 


| | 
P,~Pp=o(++2) (9) 
where r, and rz are the principal radii of curvature. 
For a spherical surface, 7; =r, and 
2 
P, - Pyp=— (10) 


For a small bubble surrounded by a liquid at 
uniform temperature, thermodynamic equilibrium then 
requires that 


T,=T, (11) 


If the saturation temperature, 7,, is defined with 
respect to the ambient or reference pressure, P.., the 


liquid will be superheated by an increment 7), — 7, in 
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order to maintain mechanical equilibrium with a 
bubble of radius r. An approximate expression for the 
superheat using the Clapeyron equation is 


20T sat. 
T, a T sat. rz dg 





(12) 


Equation 12 and a nuclei radius r = 0.003 in. were used 
to obtain the superheat values already noted for 
sodium and water in Fig. 1. 

For ~ fixed system pressure, the relation of the 
vapor and liquid states for various equilibrium-bubble 
radii is illustrated in Fig.6. Note that as 7), = Tp 
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Fig. 6 Relation between vapor and surrounding liquid states 
for various stable bubble radii. 


increases the difference between the vapor pressure and 
the system pressure increases, and thereby smaller and 
smaller bubbles are permitted to exist in equilibrium 
with the surroundings. The equilibrium state is un- 
stable, since a small increase in bubble radius or 
ambient temperature will cause the bubble to grow 
indefinitely, whereas small corresponding decreases will 
cause it to collapse. 

Particularly with respect to incipient boiling, inert 
gas will usually be present in the bubbles. This gas may 
diffuse between the liquid and vapor regions during 
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thermal transients, affecting the pressure within the 
bubble; 


P, +P, —P, == (13) 


and the amount of liquid superheat required to be in 
equilibrium with a bubble of radius r will be reduced 
accordingly. 

The well-known Hsu criterion”? for the necessary 
condition for nucleation in a flowing system has been 
extended by Bankoff and Fauske?® to include the case 
when entrained and/or precipitated gas bubbles diffuse 
from the bulk stream toward the heated wall. The 
equilibrium condition for a gas-containing bubble, 
given by Eq. 13 together with the Clausius—Clapeyron 
equation, can be used as the starting point to obtain 
the equilibrium relation for a bubble containing a 
constant number of moles of inert gas: 


mA 
AT =~ 5 (14) 


where AT = T;, — T, = bubble superheat 
r = equilibrium bubble radius 
¥ = 20T,/p, 
B =nRT),/4nXp, 
ng = (constant) moles of inert gas 


Ap, =latent heat content per unit volume of 
vapor at the saturation temperature 


This equation, shown as the solid line in Fig. 7, 
indicates that there is a maximum bubble superheat, 
AT,, = Tm — T;, at which a bubble containing a given 
amount of gas can exist in equilibrium with the 
surrounding liquid. However, below this temperature, 


AT 
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Fig. 7 Criterion for bubble nucleation. 
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the bubble may coexist stably with the liquid at more 
than one bubble radius for each liquid temperature. 
Thus the bubble is stable to small temperature pertur- 
bations below 7),,,, with respect to both growth and 
collapse, but is unstable to positive temperature fluctu- 
ations at 7; =7};,,. The broken line represents the 
temperature profile in the liquid next to the wall, 
assuming pure conduction heat transfer: 


aT=AT, -& (15) 


where A7,,, = wall superheat 
g = wall heat flux 
k = liquid thermal conductivity 


y = distance from the wall 


The intersection of this line with the equilibrium 
radius line gives the stable bubble radii (assuming no 
inert-gas mass transfer) which a bubble diffusing from 
the circulating liquid with a given initial volume (and 
hence a given gas content) will assume in the thermal 
boundary layer next to the wall. Note that two bubble 
radii are theoretically possible over a range of wall 
superheats and heat fluxes, unlike the pure vapor 
bubble, and that both of these intersection points are 
stable points of equilibrium. This means that gas 
bubbles diffusing into the wall region will either grow 
or collapse, tending to follow a path toward the nearest 
equilibrium point. The domains of attraction for each 
equilibrium point are separated by the bubble radius, 
rm, corresponding to T,,,. As the wall superheat 
increases at constant heat flux, a tangency condition 
gives the threshold condition beyond which bubble 
equilibrium is no longer possible, and the result is 
nucleation of a growing vapor-gas bubble. It is readily 
shown that this tangency condition is given by 


(16) 


6=0 corresponds to the Hsu criterion. For large k, 
Eqs. 15 and 16 reduce to AT, = AT,,,, implying that 
nucleation occurs whenever the wall temperature 
reaches T;,,,.. If one assumes a supply of entrained gas 
bubbles all with the same radius, incipient boiling will 
occur reproducibly close to the point where 
AT,, = AT,,. From Eq. 14, one can obtain the relation 


2 
AT. = 4no°T; 


m ~ IP ngR (17) 
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indicating that the wall superheat at which bubble 
nucleation can occur is inversely proportional to the 
initial mass of gas in the bubble. The mass of gas in the 
bubble is likely to depend on the system design, 
magnitude of the velocity, temperature ramp, heat 
flux, etc. In many cases, this dependence might 
therefore be largely responsible for the observed effects 
of these variables upon the incipient superheat (dis- 
cussed later). 


Heterogeneous Nucleation; Surface Cavities 


In the absence of gas bubbles in the liquid, the 
above discussion suggests that nucleation at those 
surface cavities which have not previously been fully 
wetted and flooded by the liquid metal may be favored 
as the dominant mechanism for explaining existing 
data. It is therefore necessary to examine the condi- 
tions by which vapor or noncondensable gas can be 
trapped and forced to remain inside the cavities on a 
heating surface. It is to be expected that the stability is 
strongly determined by the contact angle and the 
cavity geometry. 

Bankoff** considered the necessary conditions for 
entrapment of gas in a wedge-shaped surface groove by 
a sheet of liquid advancing across it [Fig. 8(@)] and 
showed that the gas will not be completely displaced if 
6 > 26. Even smaller contact angles can result in 
entrapment if a conical or cylindrical cavity is involved. 
Once gas is entrapped by a liquid phase, the curvature 
of the meniscus will determine whether or not the gas 
will remain trapped over a long period of time. As 
pointed out by Fabic,?° the gas will eventually diffuse 
out completely if the meniscus is convex (viewed from 
the liquid phase), if it is assumed that the ambient 
liquid is saturated with respect to gas content. Thus, 
for cylindrical cavities, the meniscus shape will be 
concave, viewed from the liquid, if the contact angle is 
greater than 90° +6 [Fig. 8(c)]. The following force 
balance can be written for the meniscus under the 
condition of equilibrium: 


Cylindrical: P, + P, + cos(m—6)=P, (18) 
; 20 
Conical: P, + Py += cos (m7 — 8 + B) = Po (19) 


Examining first the gas-free case (P,~ 0), Eqs. 18 
and 19 state that, if the terms cos (7 —@) and cos 
(7 —6 +8) are positive, vapor can be trapped inside 
cavities even though the system temperature is lower 
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Fig. 8 Illustration of trapping of gas and stability of surface 
cavities. 


than the saturation temperature corresponding to the 
liquid pressure, Py. The necessary conditions are thus 
Cylindrical: 6 > 1/2 (20) 
Conical: 6 > 2/2 + 8 (21) 
These equations also apply when the liquid is saturated 


with gas, which can be demonstrated through Henry’s 
law 


C*=K(T)Pe (22) 
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In a heated system the meniscus will be at a higher 
temperature than the bulk liquid. For alkali metals and 
most organic fluids, this implies that the constant K in 
Eq. 22 will be higher at the meniscus than in the bulk. 
If Eqs. 20 and 21 are not satisfied, this implies that 
C* > C#. where here m refers to the meniscus and 5 to 
bulk fluid. There will thus be a concentration gradient 
at the meniscus which will eventually cause all the gas 
trapped in the cavity to be removed, and the cavities 
will flood. Only if Eqs. 20 and 21 are satisfied is it 
possible that C#, = C#*. so that a stable equilibrium is 
achieved. This is because the meniscus [in the position 
shown in Figs. 8(b) and (c)] is at a higher tempera- 
ture, but lower pressure, than the bulk fluid, so that 
the two effects may cancel (see Eq. 22). Consequently, 
whether the above-type cavities remain active or not 
depends entirely upon whether Eqs. 20 and 21 are 
satisfied. However, as pointed out by Bankoff.?® it is 
not necessary to consider the effect of contact angle, if 
the geometry of the cavity is of the reentrant type. As 
shown in Fig. 9a, the minimum radius of curvature, 
upon penetration of liquid under prepressurization, 
will exist when the meniscus forms a hemisphere at the 
sharp corner at position 1. The deactivation or 
flooding of such cavities is seen to depend only upon 
subcooling and not on contact angle. For a gas-free 
system, the condition for stability can be stated as 
follows: 


fe) , 
ee an th (23) 


and, in the case that inert gas is present, as 


(24) 





’ , ' 2a( T') 
Pi — @y + Pe) << 


The above equations illustrate the importance of 
knowing the pressure—temperature history of the 
surface prior to the transient. For a well-wetting liquid, 
and if reentrant cavities are available on the surface, 
Egs. 23 and 24 will estimate the maximum size cavity 
available for nucleation. The incipient superheat can 
then be calculated from 








- 
P, -P, =) (25) 
if no gas is present and 
“ 
P, —(P;, + Py) = 2a) (26) 
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if noncondensable gas is present. Equations 23 to 26 
can be written in terms of advancing and receding 
contact angles for advancing angles greater than 90° 
(see Eqs. 18 and 19). These equations represent the 
so-called pressure—temperature history model of nucle- 
ation advanced by Fabic?® and Holtz’? and later 
modified by Chen.?” 

The model proposed by Holtz was interpreted by 
Collier?® as illustrated in Fig. 9b. The metal surface is 
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Fig. 92 Penetration of liquid into reentrant cavity. 
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Fig. 9b Formation of an active site. 


initially oxidized, and the surface is then flooded by 
the liquid metal that is forced into the cavity by the 
imposed pressure. Beyond a certain temperature, the 
liquid metal may interact with the heat-conducting 
wall according to 


M,.Oy + 2yNa + yNaO + xM (27) 


where M, O, and Na denote metal surface, oxygen, and 
sodium, respectively. This chemical reaction consti- 
tutes wetting. By further heating, nucleation will be 
initiated from the size of vapor embryo given by the 
maximum radius, r*, to which the liquid metal has 
penetrated the cavity. Very little information is 
available on the rate of reduction of metal oxide films 
by alkali metals. Furthermore, the contact angles 
required for this model are generally not known. 
Therefore, in lack of detailed information, the simpli- 
fied reentrant model proposed by Bankoff is preferred. 
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The above discussion indicates that the stability of 
a cavity is strongly dependent on the contact angle and 
history. Furthermore, for a given history, smaller 
cavities are more stable than larger cavities, and when 
the contact angle approaches 0°, only reentrant cavities 
can be stable. However, systems may exist where some 
of these observations may not be applicable; in fact, 
larger cavities may be more stable than smaller cavities 
as illustrated below. In a pool of isothermal liquid (see 
Fig. 10) where the pressure at the free surface is 
essentially the vapor pressure (as a result of continually 
bleeding the system), the following force balance can 
be written for the indicated surface cavity if perfect 
wetting is assumed (contact angle equal to 0°): 


P, +P, =P, +=2=P, + pH + (28) 


In this case the pressure inside the cavity is greater 
than the vapor pressure as a result of inert gas dissolved 
in the liquid phase. In fact, if the following inequality 
holds 


20 


Py > pg + — (29) 


& 


the cavity will be stable despite a contact angle of 0° 
Also in fact, Eq. 29 shows that, for a given liquid 
height, H, the stability increases with increasing cavity 
radius, so long as the cavity radius does not exceed the 
capillary instability radius. (Here the liquid would 
drain down the wall of the cavity and displace the gas 
bubble.) As shown later, Eqs. 28 and 29 may be useful 
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Fig. 10 Isothermal pool of sodium with a specified cylindrical 
cavity. 
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in interpreting certain available superheat data for 
liquid sodium. 

As discussed earlier in the subsection “Heteroge- 
neous Nucleation; Bubbles in Liquid,” an additional 
requirement for a given cavity of radius r, to produce 
vapor bubbles is that the temperature in the fluid 
thermal layer at a distance r, from the heater surface 
be greater than or equal to the temperature required by 
the equilibrium equation (see Eq. 12).7? Curve A in 
Fig. 11 plots this equation, and, with net heat flux, 




















A, Bubble-equilibrium eq., T vs.r, 
B, Boundary-layer temp. profile which 
satisfies Hsu's criterion 
C, Boundary-layer temp. profile which 
does not satisfy Hsu's criterion 
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Fig. 1f Hsu’s criterion for bubble nucleation in ordinary 
fluids, 


there would be a temperature profile in the fluid layer 
as illustrated by curves B or C. Hsu’s criterion then 
requires that this temperature curve be tangent to or 
intersect with curve A in order to initiate nucleation. 
The temperature profile in the liquid next to the wall 
can be estimated if pure conduction is assumed (Eq. 
15), and the minimum wall temperature necessary for 
evaporation at a given heat flux is shown by the 
tangent curve B. An excellent discussion of Hsu’s 
criterion and its applicability to liquid metals, given by 
Chen,” ” is quoted below. 
In the case of saturated water at 1 atm pressure, with a 
q of 100,000 Btu/(hr)(ft?), the model predicts that the wall 
temperature, 7,, ,,, must exceed 226°F (14°F superheat). 
The corresponding Tom would be 3 x 10% in., sufficiently 
small to have a good probability of entrapping vapor. On a 
heating surface of normal “‘as-received”’ finish, cavities of 
such dimensions would be plentiful so that incipient boiling 
would be likely to occur as predicted. 
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In the case of liquid metals, however, the situation 
becomes more complicated. Because of the higher thermal 
conductivities, the.temperature profile at a given heat flux 
tends to be much less steep for liquid metals. This leads to 
intersections of curves A and B at much greater values of 
om: Thus, for potassium at 1 atm and q = 100,000 
Btu/(hr)(ft?), Tom Would be > 4 x 10° in., with a corre- 
sponding wall superheat of 4°F. These larger cavities are 
much more likely to be flooded, especially with well- 
wetting fluids such as the alkali metals. It then becomes 
questionable whether there would exist any nonflooded 


cavities of radius r, ,,, to initiate boiling at the predicted 


superheat. This possible lack of potential sites seems to be 
born out by experimental evidence [see Fig. 3] which 
indicates incipient superheats much higher than predicted 
by the foregoing criterion. Apparently, larger cavities 
become flooded to some extent and only smaller cavities 
and/or bubbles [see the subsection “Heterogeneous Nucle- 
ation; Bubbles in Liquid’”’] remain as potential sites—thus 
requiring a superheat higher than the minimum superheat 
predicted by Hsu’s criterion. Therefore, in the case of 
well-wetting liquids (especially the alkali metals and organic 

liquids), it may be said that Hsu’s criterion imposes a 

necessary condition but does not define a sufficient or 

limiting condition for establishing incipient-boiling super- 
heat. 

However, Hsu’s criterion may be useful in deter- 
mining the cessation of boiling, in view of hysteresis 
effects which generally result in reactivation of the 
heating surface. 

The above discussion of surface cavities and the 
role they play in determining the incipient superheat is 
particularly important in pool-boiling systems. How- 
ever, in forced-circulation systems, gas bubbles in the 
liquid as well as surface cavities might be available for 
nucleation. It is important to note that the cavities 
which remain unflooded over a long period of time in 
the presence of suspended gas bubbles in the circu- 
lating stream depend upon several variables, such as 
bubble diameter, cavity shape, cavity radius, and 
contact angle. If the meniscus radius of curvature lies 
in the range of the entrained-gas-bubble radius distribu- 
tion, the smaller bubbles will tend to transport gas to 
the cavity while the larger bubbles will receive gas from 
it. The net effect will depend upon a complex rate 
process. The end results of these considerations, 
however, are that natural cone-shaped or cylindrical 
cavities or cracks may fill with liquid (assuming small 
contact angles) after long periods of operation even if 
entrained gas bubbles are present. These bubbles would 
themselves serve as nucleation centers. On the other 
hand, reentrant or sharply necked cavities, whether 
occurring naturally or artificially, would always be 
stable in the presence of entrained gas bubbles. If the 
system is in equilibrium, the active wall cavity and 
bubble radius should therefore be of the same size, 
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making it very difficult to determine the actual 
nucleation source in a convective system. In either case 
the value of the incipient-superheat temperature would 
be roughly the same. 


EXPERIMENTAL INVESTIGATIONS 


For discussion purposes, experimental investiga- 
tions concerning nucleation in liquids will be divided 
into the following categories: 

1. Homogeneous nucleation as a result of either 
slow (perfect wetting) or rapid heating (preexisting 
gaseous phase may be present). 

2. Nucleation from artificial cavities in pool sys- 
tems; no gas bubbles present in the liquid. 

3. Nucleation from naturally occurring cavities in 
pool systems; no gas bubbles present in the liquid. 

4. Nucleation in convective flows; surface cavities 
as well as gas bubbles in the flowing liquid may be 
present. 


Homogeneous or Spontaneous Nucleation 


Several sets of reliable experimental data exist for 
one-component systems involving slow heating which 
are in excellent agreement with theoretical 
predictions?®-*! (see Eq. 4). The experimental tech- 
nique used in these investigations is rather unique and 
is reviewed briefly below. 

The method consists of the generation of an 
emulsion of small drops of volatile liquids in an 
immiscible fluid medium. The individual drops rise into 
a column of the same immiscible medium, which is 
warmer at the top than at the bottom. The drops 
become warmer as they rise, and eventually the liquid 
superheats enough to nucleate and boil explosively. 
When the homogeneous nucleation temperature is 
reached, vaporization is rapid enough to produce shock 
waves. 

Figure 12 is a schematic diagram of the Pyrex 
apparatus used in Ref. 30. For hydrocarbons, sulfuric 
acid, and water, a silicone oil was used as the 
immiscible fluid medium through which the droplets 
rose. So long as good wetting (absence of an absorbed 
gas layer) existed, essentially all of the droplets of 
n-pentane and n-hexane exploded in a narrow range of 
temperatures (total range of about +0.3° for n-pentane 
and +0.5° for n-hexane). The measured temperatures at 
the point of nucleation are in excellent agreement with 
temperatures calculated from Eq. 4 based upon nucle- 
ation rates which vary by several orders of magnitude 
(see Table 1). 
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Fig. 12 Homogeneous bubble-nucleation apparatus (taken 
from Ref. 30). 


No large differences existed between small and 
large drops concerning the threshold of nucleation 
temperature, which is consistent with homogeneous- 
nucleation theory where J for these one-component 
fluids is a very steep function of temperature. This 
was illustrated (Fig. 4) for liquid sodium. 


Table 1 Measured and Calculated Temperatures 
for Different Nucleation Rates 











T(calculated) 
J = 10* J= 10° 
T(measured) cm ® sec’ cm ® sec! 
n-Pentane 147.8 + 0.3 147.7 148.3 
n-Hexane 183.8 + 0.5 183.6 184.3 
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In the cases that had a preexisting gaseous phase 
present in the form of bubbles, the drops nucleated at 
much lower temperatures and obviously did not 
nucleate homogeneously. However, if extremely rapid 
heating is used, the threshold of homogeneous nucle- 
ation might be reached even in the presence of 
preexisting gas pockets. Reference 32 reports an 
experiment using a thin platinum wire, heated rapidly 
by a pulsed current. With a heating rate of 10°°C/sec, 
the liquid in the layer near the wall was heated 
strongly, and the boiling was reported to be of an 
explosive character. Experiments were carried out with 
a number of organic liquids and with water. If the 
heating time of the wire required to reach the tempera- 
ture for spontaneous nucleation T* <79, the start of 
explosive boiling was clearly discernible. The character- 
istic time tT, depends on the fluid used and the 
preparation of the heating surface which determines 
the size and number of active nucleation sites. With 
large active cavities of high density, the degree of 
superheat required for a high rate of nucleus formation 
by liquid-density fluctuations is attained only with 
sufficiently rapid heating. If the heating rate is insuffi- 
cient, the bubbles growing on the already existing 
centers have time to coalesce in a vapor film before the 
liquid is heated to the temperature of homogeneous or 
violent nucleation. In experiments with water, Skripov 
and Pavlov found 79 © 40 usec. For the organic liquids 
they found to © 1000 yusec, which is an indication of 
the much better wettability characteristics (small con- 
tact angles) of these fluids as compared to water. If 
T*>T 9, the boiling has a comparatively quiet charac- 
ter, and at slight degrees of superheat, the wire is 
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blanketed by a vapor film. When 7* <9, the tempera- 
ture corresponding to explosive boiling, 7*, depends 
only weakly on the heating rate. With a change of 7* 
by 10 times, 7* varies by only 3 to 5°C. This is in 
excellent agreement with the homogeneous nucleation 
theory, where large variation in J has only minor effect 
on the superheat requirement for spontaneous nucle- 
ation. Table 2 (taken from Ref. 32) illustrates mea- 
sured and calculated 7* from Eq. 4 based on J* which 
was evaluated experimentally. Excellent agreement is 
noted between theory and experiment, thereby illus- 
trating that homogeneous nucleation (explosive boiling 
conditions) for a liquid can be created by heating the 
liquid from a heat-evolving solid surface. 

The characteristic time 7) for liquid sodium, as 
well as for the organic fluids, in contact with stainless- 
steel cladding and/or fuel is likely to be large in view of 
the excellent wetting properties of these materials. 
However, for hypothetical loss-of-flow accidents, 
T* > 7,9, and nucleation is likely to occur at tempera- 
tures considerably below the homogeneous nucleation 
threshold. On the other hand, in cases where molten 
fuel comes in contact with liquid sodium, the possibil- 
ity of homogeneous nucleation must be considered 
(either as a result of rapid heating or perfect wetting 
between the two liquids). 


Nucleation from Artificial Cavities 


Experiments with artificial cavities located on the 
heated surface have the advantage that their size is 
known. From measurements of the superheating the 
system can sustain, it can be determined whether a 


Table 2 Measured and Calculated 7* from Eq. 4 























T* = 35 usec T* = 100 usec 7T* = 850 usec 
Log Jt, = Log Jf, i i Log Jf, Ls se 

Liquid cm ® sec! Expt. Theory cm? sec! Expt. Theory cm * sec! Expt. Theory 
Diethyl ether 19.5 152.5 152.0 17.5 150.0 150.5 13.5 147.5 148.5 
Dipheny] ether 19.0 435.5 17.0 430.0 
Acetone 18.5 189.5 17.0 187.5 13.5 185.5 
Benzene 18.5 237.0 2353.5 17.0 232.0 234.0 13.5 232.0 232.0 
Water 15.5 302.0 SYE5 17.0 
Methyl alcohol 18.0 193.0 197.0 17.0 191.5 196.0 14.0 190.0 194.5 
Ethyl alcohol 18.5 199.0 201.5 17.0 197.0 200.5 14.5 192.0 192.0 
Propyl alcohol 18.0 2225 16.5 220.5 13.0 219.0 
Butyl alcohol 18.0 245.0 16.5 243.0 13.0 240.0 
n-Pentane 18.0 153.0 153.0 16.0 4$25 152.0 325 150.0 150.0 
n-Hexane 18.5 189.0 188.5 16.5 186.5 187.0 £25 185.5 185.5 
n-Heptane 18.5 220.5 222.0 16.5 220.0 221.0 13.0 218.0 219.0 
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given size cavity of a given geometry remains stable or 
not. Of several studies performed with alkali metals, 
the one by Schultheiss** is particularly noteworthy. 
This study exhibits excellent reproducibility and is at 
present by no means fully explained. 

The experimental equipment (Fig. 13) centers 
upon a pool of sodium above a highly polished 
horizontal test surface, 80 mm in diameter, which 
contained a single cylindrical artificial cavity. (Cavity 
sizes of 0.2 and 0.4 mm in radius were employed in 
these experiments.) The sodium was purified and its 
oxide content regulated by means of a hot trap. The 
test plate had a separately heated “hot finger” from 
the outside which was used to initiate boiling at the 
artificial cavity. Uniform temperature profiles in the 
sodium were obtained by placing the entire arrange- 
ment, with the exception of the condenser, inside an 
electrically heated furnace. 

As shown in Fig. 14, boiling was reprodugibly 
initiated from the cylimdrical, cavity at the superheat 
given by the Laplace equation (see Eq. 10), where the 
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(a) Flow diagram of superheat test rig. 
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radius of the cavity is used to determine the superheat 
vapor-pressure difference. What is remarkable about 
the Schultheiss work is that, on the basis of measure- 
ments of the wetting properties of sodium at elevated 
temperatures for various metals,?* one might expect 
the artificial cavity to be fully wetted prior to the 
initiation of boiling. In this case, in order for the vapor 
meniscus to fill the bottom of the cavity, it is necessary 
for a vapor embryo to start at some location of the flat 
bottom and then grow until it intersects the cavity 
walls. For a well-wetted cavity, this corresponds to 
forming a bubble by statistical density fluctuations, 
precisely as in homogeneous nucleation and hence 
requires a large work of formation (see the subsection 
“Homogeneous Nucleation’’). Indeed, even if the con- 
tact angle is 90°, corresponding to a very poorly 
wetted surface, these theories predict extremely high 
superheats in the absence of a preexisting gaseous 
phase. 

Rather than discard these theories, for which a 
considerable body of experimental evidence has been 
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Fig. 13. Test arrangement used by Schultheiss.* * 


REACTOR TECHNOLOGY, Vol. 15, No. 4, Winter 1972—1973 





NUCLEATION OF LIQUID SODIUM IN FAST REACTORS 





100 








Yr Trey 


bb eae 


ro) 


| ame | PTTTTy 


SUPERHEAT (Ty-Teot,), °C 








ITT] 





500 550 600 650 700 750 800 850 900 
Tsot. °C 


Fig. 14 Incipient-boiling superheat of sodium at artificial 
cylindrical cavities of 0.4 and 0.2 mm in diameter. 


built up (see the subsection “Homogeneous or Sponta- 
neous Nucleation’’), we are thus led to consider two 
possibilities in connection with this series of experi- 
ments: (1) a small amount of inert gas was stably 
trapped within the cavity for a long period of time, and 
(2) the wetting properties of sodium within the cavity 
may be considerably different than on the bulk 
surface. Turning to the first possibility and recalling 
Fig. 10 and Eqs. 28 and 29, we can make the following 
interesting observation. For the conditions of the 
Schultheiss experiments, which corresponded to a 
liquid head of about 10 cm, the required inert-gas 
partial pressure for a stable meniscus was about 0.015 
atm, even though perfect wetting was assumed. This 
would mean, if we assume that the sodium became 
saturated with argon cover gas during the loading 
process, that a small fraction remaining as residual 
dissolved gas would be sufficient to stabilize the 
presence of a gas phase in the cavity. It is well known 
that it is difficult to remove all traces of dissolved gas 
by boiling in a closed system, even after repeated 
purging. Note that this behavior would not be exhib- 
ited by natural cavities, which are several orders of 
magnitude smaller than these artificial cavities, since 
the required gas partial pressure, in order to achieve a 
stable trapped vapor gas phase, is correspondingly 
several magnitudes greater. The implications of this 
hypothesis are that artificial-cavity incipient-boiling 
superheat measurements may be misleading concerning 
the results to be expected from natural surfaces. 
Another piece of evidence in favor of the above 
hypothesis is that some of the superheats measured by 
Schultheiss were actually smaller than would be pre- 
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dicted by Eq. 10 in the absence of inert gas. Dis- 
counting the possibility of contact angles in excess of 
100° on the flat metal surfaces, there seems to be no 
other explanation than the presence of small amounts 
of dissolved inert gas in the liquid. 

Schultheiss also observed a rather remarkable effect 
of the nature of the heating surface on the incipient- 
boiling superheat in the presence of trace amounts of 
sodium oxide. Chromium and chromium-containing 
stainless steel exhibited a distinct maximum in the 
superheat with Na2O concentration at about 2 ppM. 
This behavior was not at all exhibited by nickel or iron 
surfaces, which, in fact, showed superheats below that 
given by Eq. 10 in the absence of inert gas. He noted 
that, based upon the free energies of formation, 
sodium chromite (NaCrO,) is in equilibrium with 
sodium oxide concentrations in the range of 5 to 10 
ppM O, at 700°C, whereas the oxides of iron and 
nickel were less stable and would be dissolved by the 
sodium. The implication is that in the confined space 
of an artificial cavity (or, for that matter, a natural 
cavity) the oxides of sodium would equilibrate with 
the chromium of the metal wall to form a nonwetting 
oxide surface. This would further reduce the dissolved 
inert-gas concentration required to maintain stable gas 
trapping in the artificial cavity. 


Nucleation from Naturally Occurring Cavities 


The advantage of pool-boiling experiments with 
naturally occurring cavities is that this type of cavity 
can be studied without being influenced by the 
presence of inert-gas bubbles in the liquid (as compared 
to forced-circulation systems). However, as noted by 
Singer,** the effect of inert gas dissolved in the liquid 
is still a major parameter if an inert cover gas is used to 
pressurize the system. The data discussed in Ref. 35 
and shown in Fig. 15 are particularly noteworthy since 
they represent true incipient-superheat data.* Further- 
more, a series of tests was always started from frozen 
sodium that had been melted and heated to a specific 
initial temperature and brought to boiling repeatedly, 
always starting from the same initial conditions in each 
sequence. Figure 15 shows two such test series which 


*Temperatures were measured; in addition, a low-power 
and constant direct current was passed through the test-tube 
wall and sodium, and measurements of the electrical potential 
at various locations along the length of the tube permitted the 
actual point of bubble initiation to be detected, If the location 
of the boiling initiation is not detected, one can only refer to 
maximum superheat values in the system, This may lead to 
difficulties in interpretation of data, as is illustrated later. 


REACTOR TECHNOLOGY, Vol. 15, No. 4, Winter 1972—1973 








292 NUCLEATION OF LIQUID SODIUM IN FAST REACTORS 























Hans K. Fauske 


continuous loss of inert gas from the cavities due to 


repeated boiling. As pointed out in Ref. 35, ignoring 
this effect can lead to erroneous conclusions regarding 
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Fig. 15 Variation of the incipient-boiling wall superheat with 
time and heat flux. 


employing an a-c, three-phase annular—linear induction 
pump. This apparatus is illustrated in Fig. 17. In this 
way the effect of pressure—temperature history upon 
the incipient superheat could be investigated without 
being influenced by inert-gas effects. The data are 
summarized in Fig. 18, indicating a strong effect of the 
applied deactivation pressure. The trend of these data 
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Fig. 16 Early data showing heat-flux effect upon incipient superheat. 


are identical except that the first series used a 
successively larger heat flux from test to test, whereas 
the second series involved the successive reduction of 
the heat flux. It is clear that the incipient-boiling 
superheat observed for the first test immediately after 
melting of the sodium was always much smaller than 
that required for subsequent boiling tests in the same 
series. With repeated tests the superheat is seen to 
increase, reaching a nominally unvarying value after 
three or four tests. This behavior can be attributed to a 
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is in agreement with Eqs. 23 and 25. The correct 
magnitude of the data, however, is not predicted, 
which might be attributed to the fact that the data 
correspond to maximum superheat rather than incipi- 
ent values, since the location of incipience was not 
detected. Because of the large axial temperature 
gradient present in the test tube, this uncertainty could 
be at least as large as the difference between maximum 
superheat values and predictions noted in Fig. 18. The 
possibility of superheat values larger than those corre- 
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Fig. 17 Schematic diagram illustrating experimental apparatus used by Holtz.°° 
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Fig. 18 Plot of incipient-boiling superheat vs. deactivation 
pressure. 


sponding to the noted deactivation pressures given in 
Fig. 18 can be explained by generation of high local 
pressures, following an incipient-superheat run, as a 
result of bubble collapse. Such pressure pulses were 
indeed measured during the experiments and would 
tend to deactivate a certain fraction of the test-tube 
surface beyond that corresponding to the known 
applied deactivation pressures listed in Fig. 18. This 
hypothesis illustrates the need to measure location of 
incipient boiling as well as the temperature. Further 
reference to this point is made in the discussion of 
convective superheat experiments. 

Holtz also demonstrated that boiling from a surface 
with a previous pressure—temperature history can erase 
the memory of the surface of its previous treatment, as 
indicated by subsequent boiling initiation tests. This is 
equivalent to the well-known fact that, once nucleation 
has been initiated at a heating surface, the surface 
temperature falls and the superheat required to main- 
tain boiling at a surface is considerably less than that 
needed to initiate boiling. The reason for this phenom- 
enon is that the formation of vapor patches close to 
the surface enables previously flooded cavities to dry 
out and hence become reactivated. Similar observations 
were made in connection with cavitation studies by 
Harvey et al.*” and in connection with the boiling of 
organic liquids by Bankoff et al.?® 

Harvey et al. found that prepressurization of 
normally unwetted surfaces, such as paraffin, in water 
at pressures of over 1000 atm, would result in very 
high thresholds, in terms of tensile stress exerted on 
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the liquid, for formation of vapor bubbles, and they 
attributed this to gas solution effects. Bankoff et al. 
noted a reproducible hysteresis effect in the wall 
superheat in boiling from a platinum wire and attri- 
buted this to a drift in contact angle over a period of 
several hours due to gradual loss of sorbed gas at the 
wire surface. Holtz'? postulated that the prepressuriza- 
tion effect had some unique features with alkali liquid 
metals due to their cleaning action on the cavity 
surfaces (see Fig. 9b). This leads to large differences 
between the advancing and receding contact angles 
(previously postulated by Fabic?*), with a relatively 
sharp dividing line between the “cleaned” and “un- 
cleaned” zones. Whether or not this large difference 
between advancing and receding contact angles is 
realistic remains open for discussion. However, the 
necked-cavity model proposed by Bankoff (Eqs. 23 to 
26) has the advantage of eliminating this contact-angle 
dependence. 


Nucleation in Convective Flows 


Forced-circulation sodium systems differ from pool 
systems in that entrained inert-gas bubbles are likely to 
be present. These bubbles can originate either by 
mechanical entrainment at liquid—gas interfaces in the 
system or by nucleation in heat exchangers where the 
liquid metal may become supersaturated with respect 
to the dissolved gas. Furthermore, the validity of 
applying pool-boiling data to reactor systems may be 
questioned. In several studies*®-*? of the incipient- 
boiling superheat of liquid metals during forced con- 
vection, it has been observed that an increase in the 
liquid veiocity results in a decrease in the wall 
superheat, whereas an increase in the wall heat flux 
results in an increase in the wall superheat. Several 
mechanisms have been proposed which attempt to 
explain the decrease of the incipient superheat with 
increasing velocity, among which are the entrainment 
of gas bubbles in the flowing liquid and turbulent 
fluctuations at the heated wall. However, the available 
experimental data do not conclusively prove that the 
liquid velocity per se is the cause of the apparently 
decreasing incipient wall superheat. An alternate inter- 
pretation of the reported data in which the velocity of 
the liquid and the wall heat flux are unimportant in 
determining the actual nucleation conditions has been 
suggested by Fauske.** 


In the forced-convection liquid-metal boiling exper- 
iments that have been interpreted as showing a velocity 
effect upon incipient superheat, it has been assumed 
that the maximum observed liquid superheat is, in fact, 
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that superheat required for incipient nucleation. This 
usually is associated with conditions at or just down- 
stream of the outlet of the heated section. This is a 
necessary assumption since only liquid temperature 
measurements were made. No measurements of the 
location of incipient nucleation were made. In reality, 
boiling inception may occur anywhere in the heated 
section where the liquid is superheated. If inception of 
boiling ‘s assumed to occur at some fixed location 
(“favorite cavity”) other than the end of the heated 
zone at a small superheat independent of liquid 
velocity, then the observed flow-rate effect upon the 
maximum superheat can be explained without consid- 
ering turbulence effects on nucleation phenomena. 
This is illustrated by noting that the data discussed in 
Refs. 39 and 40 can be represented by a model of the 
form, Re+ AT=constant or U* AT=constant, as 
shown in Figs. 19 and 20. The above interpretation of 
the velocity effect is in agreement with a recent 
analytical study of the turbulence effect upon super- 
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Fig. 19 Data of Chen at P, = 1.04 atm compared to the 
correlation of Pezzilli and the favorite-cavity hypothesis. 


are almost an order of magnitude greater than the 
experimental ones for a given percentage decrease in 
wall superheat at incipient boiling. 

In Refs. 41 and 42 a marked increase in superheat 
is also reported with increasing heat flux. It is apparent 
from the authors’ interpretation of measured liquid 
temperatures that it was also assumed that the location 
of incipient boiling was at the end of the heated zone. 
Again, for this case, if it is assumed that incipient 
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Fig. 20 Data of Fornari compared to the correlations of 
Pezzilli and the favorite-cavity hypothesis. 


boiling occurs at a small superheat at some fixed 
location other than the end of the heated zone, the 
data can also be represented by a model of the form 
Q/AT = constant, as shown in Fig. 21. In view of these 
observations, it can be concluded that experiments 
reported in Refs. 39 to 42 have not proven or 
unambiguously demonstrated a velocity or heat-flux 
effect on incipient superheat. 

More recent work with convective sodium systems 
has included detailed measurements of the location of 
boiling inception by measuring the local variation in 
the electrical potential along the test section. When 
their data were correlated in terms of maximum exit 
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Fig. 21 Interpretation of the observed heat-flux effect upon 
the incipient superheat. 
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superheats, Henry and Singer** noted the usual veloc- 
ity effect, which disappeared when the variation in 
location of the inception was accounted for (illustrated 
in Fig. 22). This observation supports the hypothesis 
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Fig. 22 Incipient-boiling superheats in a forced-convection 
sodium system, 


that incipient nucleation does not necessarily have to 
occur at the exit where the temperature is the highest. 
Furthermore, Henry*® analyzed the data of Ref. 45 by 
assuming that a loss of history was experienced after 
each boiling event (in agreement with Ref. 36) and that 
the deactivation conditions governing the next event 
were established by the operating procedures used 
between the runs. The system pressure was held 
constant at the boiling pressure; thus the deactivations 
for these tests were accomplished by decreasing the 
temperature. Figure 23(a) shows the comparison be- 
tween predictions from Eqs. 23 and 25 and data for 
three separate histories, and they are in good agree- 
ment over the velocity range investigated. Henry 
applied the same analysis to the “K”’ runs given in 
Ref. 39. For this series of tests, a pressure history was 
applied prior to the next boiling event. These compari- 
sons are shown in Fig. 23(b), and the agreement is 
certainly encouraging. 

As discussed earlier in this article, noncondensable 
gases can be entrained in the liquid sodium by 
dissolving gas at a high-temperature sodium—gas inter- 
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face and then cooling the sodium which drives the gas 
out of solution. This effect must be considered in 
interpreting any forced-convective superheat data. In 
the experiment of Henry and Singer, a large unheated 
bypass flow resulted in a temperature at the sodium— 
gas interface which was only 100°F above the test- 
section inlet temperature and considerably less than 
the saturation temperatures within the test section. 
Also, from the outlet of the pump to the test-section 
exit, the system pressure was more than 10 times the 
partial pressure of the dissolved gas. Therefore any 
microbubbles coming out of solution in the heat 
exchanger should collapse prior to reaching the exit of 
the test section. Nucleation due to entrained bubbles 
was therefore avoided, at least in the lower range of 
tested velocities. 

Dwyer et al.4” also measured the location of 
incipient boiling in their convective sodium-superheat 
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experiments. A temperature ramp, as well as a velocity 
effect, is reported. The data exhibiting a dependence 
on the rate of temperature rise are shown in Fig. 24, 
with the inlet temperature increasing in the order of 
1°F/min. The likely explanation for this behavior is 
probably due to the presence of gas bubbles in the 
liquid sodium. Unlike the experiment reported in 
Ref. 45, no bypass was used, and consequently the 
temperature at the free surface was essentially equal t 
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this behavior would thereby explain the trend noted in 
Fig. 24. Likewise, for an increase in the velocity, the 
inert-mass concentration present in a bubble at the exit 
would be increased since less mass is lost from the 
bubble as it travels from the inlet to the outlet of the 
test section, a behavior that suggests another explana- 
tion for the noted velocity effect. It would therefore 
appear important that, in future forced-convection 
sodium-superheat tests, provisions be made for detect- 
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Fig. 24 Illustration of temperature ramp on incipient superhea 


the test-section exit temperature, which could give rise 
to sizable gas-bubble entrainment. A brief consider- 
ation of the thermal time constants of the system 
indicates that a quasi steady state prevails at all times 
with respect to temperatures of the system, but the 
transfer of gas into or out of bubbles is diffusion 
controlled and exhibits much longer time constants. 
Equation 17 indicates that the wall superheat at which 
bubble nucleation can occur is inversely proportional 
to the mass of gas in the bubble. It follows that, for a 
constant velocity, the inert mass of the bubble would 
tend to decrease with increasing temperature ramp and 


t.4° 


ing the presence of any gas bubbles in the coolant in 
addition to measuring the location of incipient boiling. 
The method proposed by Grolmes and Fauske**® 
would appear to be sufficiently sensitive to small gas 
concentrations that it could be used for this purpose. 


APPLICATION TO REACTOR SYSTEMS 


From the previous discussion it can be concluded 
that inert-gas effects play a very powerful role in 
determining incipient-boiling superheat in alkali-metal 
loops. Although the physical laws involved are rela- 
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tively simple, the interactive effects may be complex 
and may mask other variables such as velocity, heat 
flux, and temperature ramp. From a stability point of 
view, only reentrant or sharply necked cavities appear 
to be stable because of the low contact angles 
experienced with the sodium-—stainless-steel system. 
These cavities may be in equilibrium with gas bubbles 
in the coolant stream, and either one might act as 
nucleation sites. 

After long-term operation in reactor systems, the 
effects of surface corrosion on surface sites might 
become an important consideration. Corrosion was 
experimentally evidenced by Hopenfeld*? of sodium 
on 304 stainless-steel surfaces under LMFBR condi- 
tions [sodium velocity, 20 ft/sec; heat flux, 1 x10° 
Btu/(hr)(ft?); 1200°F]. By means of electron- 
microscopy techniques, it was shown that the chro- 
mium and nickel are removed from the surface, leaving 
cavities of 40 x 10° to 80x 10 ° in. in diameter in 
the grain boundaries. In the absence of entrained gas 
bubbles, cavities of this size range would result in 
considerable superheating in a typical LMFBR system. 
A value of approximately 100 to 200°C is calculated 
from Eq. 26 if the partial pressure of inert gas in the 
cavity is set equal to the cover-gas pressure in the 
reactor. However, as discussed in Refs. 50 and 51, 
entrained gas bubbles are likely to be present in the 
reactor core. The gas bubbles may be due to mechani- 
cal entrainment at the free surface and/or gas precipi- 
tating from the liquid as it is cooled in the heat 
exchanger. The latter mechanism occurs because the 
solubility decreases with decreasing temperature, and 
the coolant is likely to be saturated with gas when it 
enters the heat exchanger. Likewise, dissolution of the 
gas bubbles will take place as they are flushed through 
the core because of the heating of sodium. However, 
Thormeier®' demonstrated that the rate of dissolution 
is relatively low by accounting for transient mass 
diffusion of gas in the liquid. For a completely blocked 
subassembly and reactor conditions typical of the 
primary circuit of the German 300-MW(e) liquid-metal 
fast breeder reactor (SNR),°? he calculated the life- 
times for both argon and helium bubbles. The results 
are shown in Fig. 25. All the helium bubbles were 
considered to be completely dissolved prior to reaching 
the saturation temperature. By contrast, argon bubbles 
above a certain size showed a different behavior. 
Because of the lower solubility and diffusion coeffi- 
cient compared with helium, all argon bubbles of initial 
size greater than 10° cm (10 y) still existed after the 
saturation temperature was reached. In the analysis of 
loss-of-flow accidents, gas bubbles appear to be avail- 
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Fig. 25 Behavior of bubbles initially filled with gas during 
heating of sodium from 380 to 945°C in 1.5 sec (Ref. 51). 


able as nucleation sites, at least if argon is used as cover 
gas, and therefore will prevent any significant super- 
heating of the liquid sodium. 

If molten fuel and liquid sodium come into 
contact, however, the possibility of high superheating 
of the sodium should be considered. Experimental 
laboratory studies have shown that vapor explosions 
can occur upon mixing of subcooled liquid sodium 
with molten UO,,°? as well as with many other hot 
and cold liquid combinations.°**’ However, there 
appears to be no consensus of opinion concerning the 
mechanism that causes the explosive vapor growth 
observed in the experiments. Indeed, there is some 
difference of opinion concerning the role of nucleation 
in explaining vapor explosions. From the previous 
discussion the following observations can be made: 
vapor explosions can occur upon mixing of a cold 
liquid with a hot liquid if the temperature of the cold 
liquid is raised to the temperature limit for sponta- 
neous nucleation. When this temperature is reached, 
vaporization is rapid enough to produce shock 
waves.”° >! If a preexisting gaseous phase is present, 
rapid heating is required to reach the homogeneous- 
nucleation temperature.*? If perfect wetting exists 
between the cold and hot liquids, slow heating is 
sufficient.?°->! The sudden contact between a hot and 
cold liquid is an example of rapid heating, where the 
contact temperature 7; will be established instanta- 
neously according to 


Tj—T:_ (Keun) ts 


Ty — T;  \ KepeCe 





(30) 
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where K is the thermal conductivity, p is the density, C 
is the specific heat, and subscripts H and c represent 
hot and cold liquids, respectively. If 7; is greater than 
the spontaneous-nucleation temperat.re, a vapor ex- 
plosion can occur. The excess temperature above 
saturation temperature (7; — T,) (calculated from Eq. 
30) as a result of sudden contact is illustrated in 
Fig. 26 for two systems. For the Al—H,O system, it is 
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occurrence of a vapor explosion between UO, fuel and 
Na are considerably more restrictive than systems like 
molten Al—H,O and that a large-scale coherent vapor 
explosion between UO,—Na is not possible. Further- 
more, fission fragments would be present which are 
likely to promote boiling prior to reaching the thresh- 
old for homogeneous nucleation in liquid sodium. 
This nucleation mechanism, however, is of little conse- 
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Fig. 26 Comparison between spontaneous-nucleation temperatures and sudden contact 


temperatures for two different systems. 


observed that the interface temperature 7; established 
upon contact exceeds the homogeneous-nucleation 
limit for H,0O, and a vapor explosion can occur 
without delay. This is in agreement with the shock- 
tube experiments reported in Ref. 58. On the other 
hand, for the UO,—Na system, the contact tempera- 
ture is considerably less than that required for sponta- 
neous nucleation of Na, and hence a vapor explosion 
due to homogeneous nucleation cannot occur immedi- 
ately following contact. However, if wetting occurs 
between liquid UO, and liquid Na (which is quite 
possible), bulk heating may result in reaching the limit 
for spontaneous nucleation. This would indicate a 
delay time after bringing two liquids like UO, and Na 
together before a vapor explosion can occur, which is 
in agreement with limited experimental data on this 
system.° > 

If homogeneous nucleation can be demonstrated to 
be a necessary requirement for producing vapor explo- 
sions, it can be said that the conditions for the 


quence in describing channel voiding with intact core 
geometry, where gas bubbles and/or wall cavities result 
in much smaller incipient superheats.! 59 »6° 


CONCLUSIONS 


e Excellent experimental verification of the 
homogeneous-nucleation theory is available. If a pre- 
existing gaseous phase is present, rapid heating is 
required to reach the homogeneous-nucleation temper- 
ature. If perfect wetting exists between the cold and 
hot liquids, slow heating is sufficient. 

e The role of nucleation in explaining observed 
vapor explosions deserves further studying. If homoge- 
neous nucleation can be demonstrated to be a neces- 
sary requirement for producing such explosions, it can 
be stated that the occurrence or a large-scale coherent 
vapor explosion between UO, —Na is not possible. 

e Inert-gas effects play a very powerful role in 
determining incipient-boiling superheat in reactor 
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systems. Tiny entrained gas bubbles are likely to be 
flushed through the core and will thereby prevent any 
significant superheating during a loss-of-flow accident. 

eIn the absence of gas bubbles, the pressure— 
temperature history is of great importance in evalu- 
ating the incipient superheat. However, further infor- 
mation is required regarding surface aging, wetting, and 
contact-angle behavior, including the rate of reduction 
of metal oxide films by sodium. 

e Examination of available data claiming strong 
velocity and heat-flux effects indicates that these 
effects are likely to be caused by inert gas present in 
wall cavities and/or as gas bubbles in the liquid. 
Future experiments should include means for de- 
tecting the presence of tiny gas bubbles in the 
coolant. 


NOMENCLATURE 


C = specific heat 
C* = concentration of inert gas in the liquid 
J=number of nuclei formed per unit volume, per 
unit time 
k = thermal conductivity 
K = Boltzmann constant or Henry’s constant 
ng = moles of inert gas 
N, = number of molecules per unit volume of liquid 
P = pressure 
q = heat flux 
r= radius 
R = gas constant 
S = surface area 
T = temperature 
V = volume of a spherical bubble 
W = work of formation 
Y = distance from the wall 


Greek Letters 


p = density 

B = apex angle 

6 = contact angle 

o = surface tension 

w = constant in Eq. 4 

A = latent heat of vaporization 


Subscripts 
b = bulk 
c =cold 
c = critical 
2 = liquid 
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&§ = gas 

H =hot 
i = interface 

m = maximum or meniscus 
s = saturation 


= vapor 
w = wall 
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Effects of Oxygen Concentration 
on Properties 
of Fast Reactor Mixed-Oxide Fuel 


By C. E. Johnson, |. Johnson, P. E. Blackburn, and C. E. Crouthamel* 


Abstract: This review assesses the current understanding of 


the effects of oxygen concentration on properties of fast 
reactor mixed-oxide fuels, The differences in fission-product 
yields from either ?7°U or ??°Pu fission greatly affect the 
oxygen potential, For fast fission, 7 ?* U produces higher yields 


of oxide formers, whereas ?*°Pu produces higher yields of 


noble metals, The effects of oxygen and plutonium concentra- 
tions on thermal conductivity of ceramic fuels are examined, 
Increases in oxygen and plutonium concentrations appear to 
produce increases in the thermal conductivity. A mathematical 
model for predicting the oxygen potential is described and 
compared with mass-spectrometric measurements, Excellent 
comparisons are found for the urania and mixed urania— 
plutonia systems. The influence of oxygen potential on 
actinide, oxygen, and fission-product transport in irradiated 
fuels is discussed. A method utilizing the Mo—MoO, indicator 
couple is given for evaluating the oxygen potential in the 
columnar and equiaxed regions of irradiated mixed-oxide fuels, 
The contributions of oxygen and fission products, especially 
cesium, to cladding corrosion are evaluated. The cladding acts 
as an oxygen sink, thereby increasing corrosion during the 
irradiation of the fuel. The effect of oxygen potential on the 
sodium—fuel reaction, in the event of a cladding breach, is 
examined, 


The purpose of this article is to assess the current 
status of our understanding of the effects of oxygen 
concentration on the chemical reactions taking place in 
irradiated fast reactor mixed-oxide fuels. In subsequent 
discussion these effects are noted as stoichiometric 
effects since they arise from the nonstoichiometry in 
the urania—plutonia system. The discussion of this 
article is limited to the mixed-oxide fuel attaining a 
targeted burnup of 10 at.% of the heavy metal content. 





*Argonne National Laboratory, Argonne, Ill. 60439. 


Early recognition of the special properties of UO, 
led to its use in the earliest reactor assemblies by 
Fermi. By 1956, bulk UO, clad in Zircaloy-2 had been 
used in its first power-reactor-fuel application. A 
review of the early research on uranium dioxide, its 
properties, and its nuclear applications, edited by 
Belle,’ was published in 1961. Enriched urania has 
been used in a variety cf thermal reactors, and it has 
always been realized that, when a significantly large 
block of the world’s power is generated by nuclear 
reactors, the supply of 7°5U will become limiting. This 
will then require the use of either *°°U or ?°°Pu bred 
from the relatively abundant *’Th or 7°°U. The 
decision was made to utilize firs: the mixed urania— 
plutonia fuel, and serious work on this mixed-oxide 
fuel began about the middle 1960s. 

In comparison with urania the mixed urania— 
plutonia fuel is significantly more oxidizing. The full 
effect of this fact has been only recently recognized. 
The steep thermal gradient and consequent oxygen- 
potential gradient which are set up in the fuel during 
irradiation have a profound influence on the fuel 
chemistry. The oxygen potential and fuel temperature 
are more strongly coupled with the fuel performance 
and life than any other parameters, with the possible 
exception of radiation damage to the cladding. The 
oxygen potential is a measure of the driving force for 
the numerous reactions that take place in the fuel 
element during irradiation. This potential does the 
following: 

e Provides the chemical driving force for attack of 

the stainless-steel cladding. 
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e Controls the chemical state of many fission 
products, their interaction with the fuel, their 
contribution to fuel swelling, their volatility, and 
their redistribution. 


e Controls the vapor pressure of many fuel com- 
ponents, especially that of the uranium oxides, 
and thus the redistribution of uranium in the 
mixed-oxide fuel matrix. 


e Affects the potential for the interaction of the 
liquid-sodium coolant with the fuel matrix to 
form a low-density sodium uranate in the event 
of minor breaches in the cladding. 


e Affects the intrinsic thermal conductivity of the 
fuel, the rate and temperature for fuel re- 
structuring that takes place early in fuel life, and 
the form of the radial oxygen-potential gradient 
which is established in the operating fuel 
element. 


Thus this article places considerable emphasis on 
understanding the generation of oxygen potentials in 
various regions of the fuel. 

With regard to the fuel—cladding interface region, 
particular attention has been given to the total oxygen 
available to this region during early restructuring and 
during the burnup period. We have attempted to 
illustrate the manner in which oxygen reacts with the 
cladding and to delineate those parameters which can 
accelerate cladding corrosion. 


There has been considerable confusion and varia- 
tion in the role assigned to the cladding in controlling 
fuel stoichiometry during irradiation. Rand and 
Markin? and Markin and Mclver? recognized that 
slightly hypostoichiometric mixed UO,—PuOQ, can 
oxidize stainless steel, but, nevertheless, the oxide 
layers observed in irradiated fuel on the inner steel 
surface did not indicate any appreciable reaction. The 
conclusions drawn?** from these observations were that 
the cladding did not react at the prevailing tempera- 
tures and that no control of the fuel stoichiometry was 
effected by the cladding. Those conclusions are at 
variance with observations at Argonne National Labo- 
ratory® on prototypical fuel for a fast test reactor. The 
Argonne data indicate that the cladding is functioning 
as an oxygen sink by uniform oxidation or by 
intergranular penetration of the cladding. The surface 
chromium oxide also reacts to stabilize cesium oxide, 
forming stable mixed oxides that function as an 
additional oxygen sink. Even the stainless steel in the 
plenum region of the fuel is observed to function to 
some extent as an oxygen sink. These observations may 
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have been made at higher cladding temperatures than 
those of earlier fuel, but it is important to include all 
the oxygen reacting in the bulk cladding which may 
not be visible as an oxide layer of appreciable thickness 
on the cladding surface. 


FISSION-PRODUCT YIELDS 


Fissile-lsotope Effects 


The injection rate of fission products into the fuel 
pin is double the fission rate. As a result the com- 
position of the fuel system changes with burnup. 
Further, the ability of the fission products to react 
with the oxygen released in fission can have a 
significant effect on the fuel stoichiometry as burnup 
progresses. 

The mass-yield distribution of the fission products 
is important insofar as it determines the amounts of 
oxygen getters relative to the more noble fission 
products. The distribution does not change appreciably 
with neutron energy within the energy range of major 
concern for fast reactor irradiations. This can be seen 
by comparison of the data for specific fissile isotopes 
in Figs. 1 and 2. In these figures the subscripts F and T 
indicate data for fast fission and thermal fission, 
respectively. 
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Fig. 1 Mass-yield curves for fast-neutron fission of 7°°U, 
737 and >>? Ba. 
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Fig. 2 Mass-yield curves for thermal-neutron fission of ?°*°U, 
22°C ane >>" Pu, 


From Figs. 1 and 2 a major shift in the mass yields 
is evident for **°Pu as compared with 7°°U or 7°°U for 
both fast- and thermal-neutron energies. In both cases 
a significant drop occurs in the zirconium isotopes, 
from a total yield of 30 to 20%, with corresponding 
increases in the noble-metal fission products, i.e., 
rhodium and palladium. As a result, when *°Pu is the 
fissioning fuel, significant increases in the quantity of 
available oxygen are to be expected with increasing 
burnup. The present test programs on fuels irradiated 
in fast reactors have utilized fully enriched 7*°UO,—20 
at.% *°°PuO, mixed-oxide fuel to achieve prototypical 
power densities. However, future large breeder reactors 
will utilize mixtures of plutonia and natural or 
depleted urania, which will increase the relative 
amounts of 7*Pu fission. Also, in future recycled 
plutonium fuel, a significant amount of fission will 
come from **'Pu. The data available for 7*’Pu fission 
yields indicate a fission-product distribution very 
similar to that for 73Pu. Table 1 gives the *°°U and 
29Pu total-fast-fission yields for the isotopes of each 
element. These yields are needed to calculate the 
available oxygen as a function of burnup. 

In thermal reactors with large-diameter fuel pins 
(~%, in.), there is a very significant gradient in the 
burnup between the center and periphery of a fuel pin. 
This behavior contrasts significantly with fast reactor 
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fuels in that, with a smaller diameter pin and a harder 
neutron spectrum, burnup is uniform across the fuel 
radius. 

Utilizing currently available fast flux irradiation 
facilities, it is very difficult to design an experiment 
that will reproduce the effects of plutonium fission in a 
fuel pin operating with prototypical internal thermal 
gradients. It is important in such an experiment that 
the internal fuel-pin temperatures be reproduced, since 
the thermal gradient and oxygen poiential are the 
driving forces for redistribution of material to the 
cooler periphery of the fuel pin. The initial restruc- 
turing period also sets up an oxygen-potential gradient 


Table 1 Elemental Fast-Fission Yields 
(In Pile 365 days, Out of Pile 180 days) 














235 239py 
Fission Fission 
Element yield, % Ref. yield, % Ref. 
Se 0.6 6 0.4 6 
Br 0.4 6 0.3 6 
Kr 37 7 1.9 7 
Rb 3.9 7 1.8 7 
Sr 9.1 7 3.7 7 
Y by ? 1.8 * 
Zr 30.8 + 19.1 7 
Nb 0.0 6 0 6 
Mo 27.0 7 22.9 7 
Tc 5.8 8 6.0 6 
Ru 11.8 ! 23.1 7 
Rh Bo * 7 * 
Pd 4:5 6 14.8 6 
Ag 0.03 6 2 6 
Cd 0.2 6 l 6 
In 0.02 6 0.1 6 
Sn 0.1 6 0.3 6 
Sb 0.1 6 0.2 6 
Te 1.0 6 1.0 6 
I 2.0 6 1.5 6 
Xe 20.6 7 9 By 7 
Cs 19.4 7 21.0 7 
Ba 6.6 7 5.0 t 
La 6.1 8 5.3 - 
Ce 12.7 8 11.9 7 
Pr 6.0 8 5.3 7 
Nd 18.7 7 15.0 7 
Pm 2.4 3 ya * 
Sm 1.1 7 3.2 7 
Eu 0.2 > 0.8 6 
Gd 0.8 * 0.4 6 





*Values estimated from adjacent-chain mass yields. 
+Thermal values from Ref. 7. 
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in the fuel with the cooler outer region increasing in 
oxygen potential. This increase is, however, limited by 
reaction of oxygen with the cladding to produce 
slightly hypostoichiometric fuel compositions. There- 
fore the introduction of extra oxygen into fully 
enriched fuel would not be expected to produce the 
same result as the equivalent additional oxygen gen- 
erated in plutonium fission. Probably the closest 
reproduction of the actual results of plutonium fission 
in future liquid-metal-cooled fast breeder reactors 
(LMFBRs) would be made in experiments employing 
higher plutonium fuel content and slightly larger 
diameter pins in conjunction with flux peaking in an 
Experimental Breeder ReactorII (EBR-II)  sub- 
assembly. 


Reactor Environment Effects 


The reactor environment in large LMFBRs will 
encompass a broad neutron-energy spectrum which will 
shift in different regions of the reactor core. The future 
fuel will be natural or depleted U—Pu mixtures in 
which the competition between 7°8U neutron capture 
to produce 7°°Pu and the depletion of **Pu by fission 
will be counterbalancing. The differential equations 
expressing the time rate of change of the atom content 
are as follows: 


238 











ill iaccill a 
aN = oN — o3°6N™ @) 
IN = o4N® — o39N™ (3) 
_- = 92% gN™ _ g241gy24t _ 241241 (4) 


dN242 
= o2419N*4 ~ of” N?42 =. o22oN2” (5) 


dt 
xNo — UN 
Burnup = —*—— x 100 (6) 
. 
where the superscript identifies the atoms of 


uranium (238) or plutonium (239, 240), ¢=the 
neutron flux, 0, = capture cross section, of = fission 
cross section, and 0g = d¢ + of. 

A number of experimental reaction rates for 
various fast reactor spectra have been measured 
directly by Dudey,? and, using these experimental (0¢) 
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reaction rates, he calculated *°Pu/?°8U ratios as a 
function of burnup for three test cases of 
239Py/ 258 = 0.25, 0.33, and 0.50. This amounts to 
solving Eqs. 1 and 2 using the experimentally deter- 
mined reaction rates that are for the spectral locations 
described in columns | and 2 of Table 2. 


The median neutron energy of the core center of 
EBR-II is 500 kV, whereas out at row 7 the median 
neutron energy is <100 kV. Because the core center of 
EBR-II is harder (higher energy) than the core center 
of the Fast Flux Test Facility (FFTF), the plutonium 
content will decrease less for FFTF than for EBR-II at 
a constant burnup. However, in the FFTF reflector 
edge, the spectral conditions have reversed the trend to 
produce a continual increase in 79?Pu/?9*U ratios 
through 10 at.% burnup. In future large LMFBRs as 
the spectral conditions produce more moderated neu- 
tron spectra, this result would be larger increases in the 
239Py/?38U ratios as burnup progresses. In turn, this 
behavior is expected to begin to play a more important 
role in affecting oxygen potentials of the fuel as the 
fluxes are moderated in future large LMFBRs. 


THERMAL GRADIENT OF OXIDE FUEL 


The use of U—Pu oxide as a fast ceramic reactor 
fuel requires a thorough understanding of the thermal 
behavior of this material. The large temperature 
gradients that occur in a fuel during irradiation 
contribute directly to the variety of structures and 
elemental redistributions that have been observed in 
postirradiation studies. Transport of material down the 
temperature gradient and migration of pores up the 
gradient alter the initially uniform distribution of 
porosity, uranium, plutonium, and oxygen in the fuel. 
The occurrence of equiaxed and columnar grain growth 
also results in variation of density and grain-size 
regions. All these conditions produce a complex fuel 
matrix that can appreciably alter the heat transfer 
across the fuel pin. Postirradiation examination and 
analysis require knowledge of the radial temperature 
distribution in mixed-oxide fuels. The effects of 
stoichiometry, composition, porosity, etc., on the 
thermal conductivity need to be understood before 
reliable temperature calculations can be made. The 
thermal-conductivity data generally fall short of 
dealing with the various conditions and temperature 
ranges found in irradiated fuel. Many of the current 
data are in the form of thermal-diffusivity measure- 
ments that require heat-capacity data for conversion to 
thermal conductivity. 
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Table 2 Calculated 77° Pu/?**U Ratios and 7°? Pu Relative Atom Content for 
Various LMFBR Spectra Based on Measured Reaction Rates 





Initial test-case conditions 





239 Py atoms = 1 
238) atoms= 1 


239 Pu atoms = 1 
238) atoms = 3 


239 Py atoms = 1 
238 atoms = 4 











Fluence, . 
Reflector arbitrary Burnup, 7°°Pu/ 7°°Pu, Burnup, 7°°Pu/ 73°Pu, Burnup, ?°°Pu/ 23°Pu, 
conditions Location units % 238 = atoms % 238U atoms % 238 atoms 
EBR-II Row 2 0.1 2.27 0.2352 0.9304 2.72 0.3114 0.9240 4.99 0.9214 0.9112 
full-power Axial = 
test, Run —1.6 cm 0.3 6.41 0.2092 0.8097 7.63 0.2731 0.7925 13.73 0.7836 0.7580 
50H Radial = 
7.2 cm 0.5 10.10 0.1875 0.7098 11.93 0.2409 0.6840 21.07 0.6682 0.6324 
EBR-II Row 3 0.1 2.35 0.2347 0.9282 2.80 0.3110 0.9224 5.04 0.9212 0.9107 
full-power Axial = 
test, Run -1.9 cm 0.3 6.63 0.2080 0.8036 7.84 0.2718 0.7879 13.87 0.7830 0.7565 
50H Radial = 
10.2 cm 0.5 10.45 0.1855 0.7007 12.25 0.2390 0.6771 21.29 0.6673 0.6301 
EBR-II Row 7 0.1 2.10 0.2388 0.9437 2.57 0.3152 0.9340 4.88 0.9259 0.9145 
full-power Axial = 
test, Run —1.6 cm 0.3 5.99 0.2192 0.8452 7.24 0.2833 0.8191 13.47 0.7957 0.7669 
50H Radial = 
33.5 cm 0.5 9.52 0.2028 0.7628 11.38 0.2565 0.7237 20.70 0.6865 0.6456 
FFTF-EMC Core 0.1 2.16 0.2445 0.9611 2.62 0.3212 0.9470 4.92 0.9350 0.9189 
mockup center 0.3 6.24 0.2347 0.8911 7.47 0.2998 0.8536 13.62 0.8201 0.7784 
in ZPR-9 0.5 10.05 0.2264 0.8303 11.88 0.2816 0.7744 21.04 0.7226 0.6625 
Assembly 27 
FFTF-EMC Reflector 0.1 2.02 0.2578 1.004 2.49 0.3352 0.9792 4.83 0.9547 0.9296 
mockup edge 0.3 6.07 0.2717 1.003 7.31 0.3385 0.9377 13.51 0.8734 0.8065 
in ZPR-9 0.5 10.09 0.2837 0.9935 11.92 0.3414 0.8967 21.06 0.8033 0.7033 


Assembly 27 





Although the thermal properties of UO, have been stitutional atom causes a decrease in thermal con- 


extensively investigated, only a limited number of 
thermal-property measurements have been reported for 
mixed U—Pu oxides. Much of the early work has been 
summarized by Hetzler et al.,’° and more recent work 
has been given by Gibby,'’! Bailey et al.,’* and 
Laskiewicz et al.'* Although these experimenters used 
a variety of experimental techniques, their results are 
in reasonable agreement. 

Gibby’? determined the thermal diffusivity of a 
(Up.75PUp.25)O2.. mixture using a laser heat pulse 
technique for oxygen stoichiometries in which x varied 
from 0.00 (stoichiometric) to 0.07 [oxygen-to-metal 
ratio (O/M) = 1.93]. These data were then con- 
verted to thermal conductivity using estimated heat- 
capacity data. These data are in reasonable agreement 
with earlier steady-state measurements. The results of 
this investigation indicate that the thermal con- 
ducitivity of the mixed oxide is quite similar to that of 
pure UO, (Ref.15). The introduction of a sub- 


ductivity if there is a mass difference, an atomic size 
difference, or a bonding difference. Since tetravalent 
uranium is so similar to tetravalent plutonium in these 
respects, little change is expected. 

Gibby'! also examined the effect of changes in 
plutonium content on the thermal conductivity of 
(U,Pu)O, solid solutions and found that the con- 
ductivity decreases with plutonium content in a sys- 
tematic manner for mixed-oxide compositions contain- 
ing up to 30 wt.% plutonia. Comparisons of this study 
with studies on the effect of plutonia content on 
thermal conductivity of mixed-oxide solid solutions 
show that, in the temperature range 100 to 1200°C, 
oxygen stoichiometry has a much greater influence on 
thermal conductivity than does the plutonium content. 
Calculated values for higher plutonium concentrations 
indicate that the minimum thermal conductivity in the 
UO,—PuO, systems occur at a PuO, concentration of 
approximately 70 wt.%. 
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Laskiewicz et al.'* recently reexamined the thermal 


conductivity of (U,Puy)O,, by extending Bailey’s 
work to 2300°C and including the effect of O/M. The 
experimental technique involved heating annular 
pellets by a tungsten rod and measuring the tem- 
perature in blackbody holes using an optical 
pyrometer. The results exhibited an increase in thermal 
conductivity as the stoichiometry increased at low 
temperatures. This is in accord with the thermal- 
diffusivity measurements of Gibby'’ and Van 
Craeynest and Weilbacher,’® who found a monotonic 
increase in thermal diffusivity as the stoichiometric 
content was approached. The effect of density change 
on thermal conductivity of the mixed oxide was much 
less than in previous studies using uranium dioxide. 
The effects of varying density and plutonium com- 
position were overshadowed by the effects of varying 
stoichiometry. 

There is general agreement’’’'* on the strong effect 
of stoichiometry on the thermal conductivity. At tem- 
peratures much greater than the Debye temperature, 
the thermal conductivity is largely dependent on the 
phonon mean-free-path length. As the path length is 
reduced, phonon scattering becomes more pronounced, 
and the result is a decrease in thermal conductivity 
following a reciprocal dependence on temperature. 
When the mean free path is reduced to interatomic 
distances, further reduction in path length has no 
effect and the thermal conductivity becomes tempera- 
ture independent. This appears to be the case for 
hypostoichiometric urania—plutonia mixtures above 
1400°C. Below this temperature the thermal con- 
ductivity increases with decreasing temperature and is 
highly dependent on such impurity-scattering centers 
as oxygen vacancies. The phonon mean free path is 
reduced to interatomic distances with only about 3% 
of the oxygen sublattice replaced by vacancies, 

The thermal conductivity of a porous material 
differs from that of the 100% dense material in a 
manner dependent upon the pore fraction, con- 
ductivity, shape, and distribution. These factors are 
usually expressed in terms of a porosity correction 
factor relating the actual conductivity to that of the 
100% dense material. Belle et al.'7 reviewed the 
theories and indicated that the Maxwell—Eucken 
equation, 


Kp/K yoo = (1 — P)/(1 + BP) 


where K,p = effective 
medium 
K io = conductivity of the 100% dense material 


conductivity of the porous 
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P= volume fraction of pores 
6 = pore shape factor 


is the most accurate representation of the effect of 
porosity on thermal conductivity. Marino’® has 
examined in detail those factors which influence the 6 
parameter. For a porosity range of 0.1 to 0.25, this 
parameter was found to vary between 0.5 to 1.6. 

The work of Biancheria,’? Gibby,'’ and Laskie- 
wicz’* are in general agreement and cover the tem- 
perature range from 200 to 2500°C, which is generally 
more than adequate for present needs. Meyer and 
Buescher” recently presented a simplified mathe- 
matical approach for the utilization of available 
thermal-conductivity data for temperature calculation. 


OXYGEN POTENTIALS IN U—Pu 
OXIDE FUEL 


One of the most important parameters for the 
oxide fuel is the oxygen potential (defined as 
AGo, =RTInPo,) asa function of temperature and 
concentration of oxygen and plutonium. Some of the 
work with U—Pu oxide fuel-irradiation experiments 
suggests that the differences between stoichiometric 
urania and urania—plutonia mixtures, with O/M = 2, 
were not fully appreciated. As a result, early in the 
LMFBR program, irradiations were carried out of fuel 
having O/M = 2, and fuel specifications were set near 
this composition. 

There are two essential differences between urania 
and U—Pu oxide which involve the effects of oxygen 
concentration on fuel behavior. First, for a given O/M 
ratio and for a given plutonium concentration, the 
U—Pu oxide has a much higher oxygen equilibrium 
vapor pressure or more positive oxygen potential (i.e., 
plutonia is a more oxidizing material than urania). 
Depending on temperature and composition, the 
oxygen pressure may be many orders of magnitude 
higher than that for uranium oxide of the same oxygen 
concentration. For example, the metal-rich boundary 
of UO, at 1000°C is at an O/U ratio of 1.998. At this 
oxygen concentration the oxygen pressure over 
Up sPUo,201,998 is 10'7, 107, and 10? times that over 
UO, 99g at 1000, 2000, and 3000°C. These differences 
in pressure are greatest for O/M ratios just below 2 and 
at lower temperatures. 

The second difference between UO, and U—Pu 
oxide concerns the shape of the oxygen-pressure or 
oxygen-potential curve. For UO,,,, there is an ex- 
tremely large change in oxygen pressure with com- 
position near the stoichiometric composition (see 
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Fig. 3). At lower temperatures the change is so sharp as 
to appear discontinuous. The effect of this change is to 
keep UO, essentially at an O/U ratio of 2 even though 
the oxygen potential may vary extensively. Con- 
sequently, since all postulated mechanisms for oxygen 
redistribution depend on the oxygen-potential 
gradients, the oxygen concentration in UO, is not 
likely to change radically from the center of the fuel to 


the cladding. Nor is there likely to be much reaction of 
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Fig. 3. Calculated oxygen potentials at 2240°C for UO,.,, 
PuO, ., (fluorite single phase), and U Pu O,+x- 
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oxide with cladding, since a small amount of oxidation 
will reduce the oxygen potential of UO, to a value in 
equilibrium with the cladding. The effects of burnup 
are more complex, but, in general, UO, would appear 
to be less oxidizing for a given burnup because of the 
more reducing nature of **°U fission products as 
compared with those of 7°°Pu (see the sections on 
“Fission-Product Yields” and ‘“‘Material Distribution in 
a Thermal Gradient’). 


For U—Pu oxide the oxygen-pressure or oxygen- 
potential curve is flatter and, of course, higher than 
that for UO,. The result is expected to be a broader 
oxygen gradient from the center of the fuel to the 
cladding. Also, more cladding oxidation by U—Pu 
oxide is likely to occur for a given oxygen-potential 
change because a larger change in oxygen concentra- 
tion is required than is the case for UQ,. 


Experimental Oxygen Potentials 


Several investigators have measured oxygen poten- 
tials of U—Pu oxide, but the agreement between these 
measurements is poor. Markin and Mclver® used the 
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electromotive-force (emf) method to measure oxygen 
potentials of Up g9Puo 1102+, Uo ssPUo, 15922 x, and 
Up. wPup.39002+. between 800 and 1100°C. Oxygen 
potentials were measured for O/M ratios of 1.87 to 
2.06 and 1.95 to 2.08, respectively, for the 30 and 11 
wt.% Pu—U oxides. To calculate oxygen potentials for 
other plutonium concentrations, Markin and Mclver? 
assumed the substoichiometric U—Pu oxide could be 
treated as consisting of U*, O?, and tetravalent or 
lower plutonium cations. The hyperstoichiometric 
oxide was assumed to consist of Pu*, O7?, and 
tetravalent or higher uranium cations. The measured 
oxygen potentials of the MO, material were found to 
correlate with the average valence of plutonium for all 
three plutonium concentrations. For MQO,., the 
measured oxygen potentials were correlated with the 
average uranium valence for the 11 and 30 wt.% Pu-U 
oxides and with UO + x. 

Woodley *' equilibrated the sample with H,0/H, 
ratios and measured the oxygen potential of U—25 
wt.% Pu oxide from 950 to 1400°C. The measurements 
covered O/M ratios of 1.90 to 1.99. The coprecipitated 
oxide was calcined at 800°C as compared with 1500°C 
for the Markin and Mclver? sample. This temperature 
is certainly too low to form a solid solution; however, 
since measurements were carried out up to 1400°C, a 
solid solution may have formed in later experiments. 
Woodley’s and Markin and Mclver’s* oxygen poten- 
tials agree at 1100°C, where the data of the latter 
authors are interpolated by means of their plutonium 
valence correlation. However, because of a 20 to 30% 
difference in the oxygen partial molar enthalpies, these 
two investigations disagree at other temperatures. 

Javed and Roberts”* measured oxygen potentials of 
U—20 wt.% Pu oxide samples from 1000 to 1700°C 
with O/M ratios of 1.92 to 1.99. Their samples were 
made from sintered mechanically blended powders 
consisting of 0.6-4 UO, and 20- to 30-u PuO, particles. 
Using electron probe microanalysis, they found that 
some U—Pu inhomogeneity existed even for the ma- 
terial heated at 1700°C. The temperature dependence 
of the Javed and Roberts work is intermediate between 
that of the prior investigators, but the oxygen poten- 
tials are much more negative (10 to —20 kcal depend- 
ing on temperature and composition). In addition, 
Javed has measured oxygen potentials of UO, ,, 
finding values about 10 kcal more negative than other 
measured oxygen potentials, suggesting either a nega- 
tive bias in oxygen potentials or a positive bias in the 
measured oxygen concentration. 

Battles et al.?* used a mass spectrometer to measure 
pressures of UQO,, UO3, PuO, and PuOQ, over 
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Up gPuo.202. from 1600 to 2100°C and for O/M ratios 
of 1.92 to 1.98. Recently measured ionization cross 
sections** for UO, and UO; have been used to calculate 
oxygen pressures from the Battles et al. measured 
ion-intensity ratios of UO; to UO} together with 
thermodynamic data for the precursor**” gases. The 
calculated oxygen potentials fall between the extrapo- 
lated potentials of Markin and Mclver and those of 
Woodley and are about 20 kcal higher than those 
extrapolated from Javed and Roberts. 

Another set of mass-spectrometric measurements 
Of Up.gsPU9.1502-~ was carried out by Ohse and 
Olson.*® These authors measured ion intensities and 
ion-sensitivity factors of UO,, UO3, PuO, and PuQ,. 
Use of the Blackburn ionization cross sections for UO, 
and UO; ions, along with the Ohse and Olson ion 
intensities and the same thermodynamic data for UO, 
and UO, gas, yields oxygen potentials that are 5 to 10 
kcal more negative than those calculated from the 
Battles data. 


The experimental oxygen potentials are widely 
scattered, differing by as much as 20 kcal, and they 
also exhibit great differences in temperature 
dependence. In these measurements two variables are 
especially troublesome and deserve comment. First is 
the analysis for oxygen concentration in the U—Pu 
oxide. Most of the above-mentioned experimentalists 
used either the Lyon”? or the McNeilly and Chikalla® 
method. Battles’? used gas-fusion analysis. If the 
differences among these five sets of oxygen potentials 
are attributed solely to uncertainties in oxygen con- 
centration, the spread in O/M ratios is no greater than 
0.03. 


In a recent study of analytical methods for FFTF 
fuel by Rein, Zeigler, and Metz,” identical samples 
were sent to seven laboratories for analysis. All but one 
of the oxygen analyses were obtained by the Lyon 
method in which samples are oxidized above an O/M 
ratio of 2 and reduced in a H,—H,O0 mixture to the 
dioxide. One laboratory carried out this reaction in one 
step with a H,—-H,O mixture (the Chikalla method). 
These analyses on a total of 25 samples yielded O/M 
ratios between 1.963 and 1.994. Hence, even under 
very carefully controlled conditions, it is possible to 
obtain a variation in O/M ratio of as much as 0.03. 


The other difficult variable is the temperature 
dependence of the oxygen potential. The thermo- 
dynamic literature is filled with large errors in the 
“second law slope.” A 20 to 30% spread in AHo, is not 
unusual. Thus there is no justification for accepting 
extrapolated data over direct measurements. 
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Thermodynamic Model for U—Pu Oxide 


Blackburn®™ developed a thermodynamic model for 
calculating the oxygen potentials in a mixed urania- 
plutonia system. The model makes use of the observa- 
tions by Markin and Mclver® that the oxygen potential 
can be correlated with the average plutonium valence 
(assuming uranium is tetravalent) for MO,., and with 
the average uranium valence (assuming plutonium is 
tetravalent) for MO,,,. 

However, instead of using the empirical relations 
derived by Markin and Mclver, Blackburn relates the 
cation concentrations to phase-boundary and integral 
thermodynamic values for the oxide—fluorite phase. 
Furthermore, the model covers a much broader com- 
position and temperature range. 

The model was developed in a stepwise fashion 
starting with UO,:,, then adding plutonium, and 
finally incorporating fission products. Since UO, con- 
stitutes the bulk of the fuel, a viable model must 
predict accurate oxygen pressures over metal- and 
oxygen-rich urania. The principles involved can be 
demonstrated using only the UQ,,, system. 

The approach used is based primarily on the law of 
mass action. To develop the model, formulate ex- 
pressions for both metal- and oxygen-rich UO, and 
then combine the two sets of equations. The combined 
equations can be used to calculate oxygen pressures, at 
a given temperature, for any composition from the 
UO... boundary to UO, ,, or the equations can be used 
to calculate compositions from oxygen pressures. 

For UO,.,., the following equilibrium is assumed to 
exist between cations and anions in the lattice and 
oxygen gas: 


2U* + 207 = 2U* + 0, (7) 
For this equilibrium, one may write 


nye 


ny 2+ 





In pg = 2 In +2Inno2z + Ink (8) 


where the n’s are moles of tetravalent and divalent 
uranium and oxygen per mole of UO,. By use of a 
Gibbs—Duhem equation, one may derive an expression 
for the activity of uranium. When the boundary 
conditions are satisfied [ie., a, =1 at the UO,, 
metal-rich boundary and In (@y) <9 t1n Po,)x=0 
AGUO,,c)/RT | , Eq. 8 becomes 


ny* 





In Pg, = 2 In + In ng2- — 156,600/T + 27.2 (9) 


is | hal 
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where the temperature dependency includes both the 
enthalpy of formation of UO, and twice the tempera- 
ture dependence of In x at the phase boundary 
(UO,.x). 


For UO,+,, an equilibrium is assumed for 
2U% + 207 = 2U" +0, (10) 


By a process similar to that used for UO2 x, an 
equation is derived for the oxygen pressure in terms of 
the concentration of tetravalent and hexavalent ura- 
nium, 


nyo 





In Po, = 2 In + 2Inagz — 33,000/T + 10.4 (11) 


Ay + 


In this case the temperature dependency includes the 
difference, 2AH7U,0, — 8AHfUO,, and the tempera- 
ture dependence of 2 In y for the boundary (UQ,+,,). 
Equations 9 and 11 and two additional equations, — 


ny? + ny* + nyt = 1 (12) 
and 


ny2 + Iny + 3nys =Nor , (13) 


together with the temperature and either the oxygen 
pressure or the oxygen concentration, determine the 
system. 


The model has been used to calculate oxygen- 
pressure curves at 2705, 2390, and 2080°K for 
comparison with measured oxygen pressures. The 
curves for the atomic oxygen pressure at the three 
highest temperatures and the phase boundary are 
plotted in Fig. 4 with Tetenbaum and Hunt’s*® mea- 
surements over UO,_,.. The model shows an excellent 
fit with experimental data for UO,4,, up to ~UQ, ;. 

A similar treatment was used for PuO,.,. However, 
for plutonium, the cations consist of Pu”, Pu®*, and 
Pu**; and the concentration dependence of these 
cations On oxygen pressure was obtained from a fit to 
experimental oxygen pressures over PuO,_,. The result- 
ing equations are 


np,* 





In Po, = 4 l|n + 2 Innoz 


Apy* 


_ 109,200 


14 
7 +25.6 (14) 
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Mp,,3* 
In po, = 4 In re .+ 2 in Nor 
2 Ap y2* 
= 188,000 | 45 (15) 
i i 
Apy* + Npys + Npy = l (16) 
and 
2npy* + 1.5 Apy>* + Mpy® =No2- (17) 


For the U—Pu oxide, the eight equations (Eqs. 9 
and 11 to 17) are combined such that ng2- is the total 
oxygen for U—Pu oxide, and the sums of the uranium 
and plutonium cations are each equal to the uranium 
and plutonium concentrations in the fuel. Figure 3 
shows the individual curves for UO zs,., PuO,, 
(fluorite lattice) and the calculated curve for 
Up, gPUp,.0>+ » at 2240°K. 

Figure 5 shows calculated curves of oxygen poten- 
tial vs. temperature for Ug gPUg 05+, at O/M ratios of 
2.05, 2.01, 1.98, 1.95, and 1.91, where they are 
compared with the experimental measurements de- 
scribed above. Note that the curves are not straight 
lines, which indicates, of course, that linear extrapola- 
tions are not valid. 

The curves in Fig. 5 show that the oxygen 
potentials calculated with the model are close to an 
average for the experimental values and hence are 
believed to be a reasonable substitute for the measured 
data. This is especially true for oxygen potentials at 
temperatures and compositions outside the range of 
the measurements. The model for calculating the 
oxygen potential is invaluable in calculating such other 
fuel properties as actinide and fission-product redis- 
tribution (see the next section), O/M ratios from 
oxygen potentials based on molybdenum distribution 
(see the next section), oxygen content in fuel in 
equilibrium with the sodium—fuel reaction product 
Na3;MO, (see the section on “Attack on Fuel 
Cladding”’). 


MATERIAL DISTRIBUTION INA 
THERMAL GRADIENT 


Oxygen and Oxidized Fission Products 


When a virgin fuel pin is brought to power, the 
radial and axial temperature gradients developed lead 
to gradients in the oxygen potential, the partial 
pressures of the various gaseous-fuel-oxide species 
(UO,, UO;, PuO,, PuO), and the activities of the 
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Fig. 4 Calculated and observed atomic oxygen pressures over UO, + y. 


uranium, plutonium, and oxygen ions in the solid 
oxide. These gradients cause a tendency for the 
oxygen, uranium, and plutonium to redistribute. As 
the concentrations of fission-product elements increase 
owing to burnup, some of these elements will also 
deviate from a uniform distribution and participate in 
the transport processes in the fuel. Because of the 
gradual change in power level and overall composition 
due to burnup, a steady state is probably never 
attained in the fuel pin. Since a combination of 
solid-state diffusion and gaseous transport processes is 
involved, the concentration profiles at a given power 
level and time of irradiation will be influenced by the 
porosity of the oxide, which depends on the initial 
density of the oxide. The cladding can also react with 
oxygen or fission-product elements and therefore must 
be considered. These processes can lead to a non- 
uniform distribution, in the fuel pin, of oxygen, 
uranium, plutonium, and some of the fission-product 
elements. 

Estimates of the oxygen-potential gradient in a 
radial temperature gradient for a urania—plutonia solid 
solution have been made by Rand et al.” and Aitken 
et al.*° The estimate of Rand and Roberts™ is based on 
the assumption that a constant CO,/CO or H,0/H, 
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ratio exists in the temperature gradient. Aitken et al. 
base their estimate on experimental results obtained 
when oxide pellets were held in an axial temperature 
gradient. Both groups of workers considered qualita- 
tively the probable effect of burnup on the oxygen- 
potential gradient. A method*® for the estimation of 
the oxygen-potential gradient in the oxide phase of 
fuel pins from postirradiation studies recently has been 
developed. This new method depends on_ the 
determination, by electron microprobe analysis, of the 
ratio of oxidized to reduced fission-product 
molybdenum as a function of position in a cross 
section of an irradiated fuel pin. This method yields 
the oxygen-potential gradient in the oxide phase when 
under irradiation and therefore includes the effects of 
burnup and any oxidation of the cladding that may 
have taken place. 


Oxygen Gradient of Urania—Plutonia Solid Solu- 
tion at Zero Burnup. It is believed that, early in the 
irradiation period of a fuel pin and concurrent with the 
establishment of the fuel structure, a redistribution of 
fuel constituents occurs, resulting in O/M and Pu/U 
ratio gradients. This distribution is believed to be 
established before a substantial amount of burnup 
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Fig. 5 Calculated and measured oxygen potential vs. tempera- 
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occurs, and therefore fission-product elements are not 
involved. Rand and Roberts™ indicated that the carbon 
» hydrogen impurities in the oxide in the form of CO, 
yw HO ‘rs In this case they 
suggested that the composition of the oxide would 
tend to follow a line of constant H,O/H, or CO2/CO 
ratio. Rand and Markin? reported the results of the 
quantitative application of this concept. To apply this 
concept, use the following to compute the oxygen 
potential as a function of radius: the radial tem- 
perature profile for the oxide, an assumed value for the 
CO,/CO or H,O/H, ratio, and the free-energy data for 
the CO,/CO or H,O/H, system. Using the relation 
between oxygen potential, O/M ratio, and temperature 
for the urania—plutonia solid solution, the value of the 
O/M is computed as a function of radius. The average 
O/M of the oxide is then computed from the material 
balance. It is assumed that, when an oxide with this 
overall O/M ratio is placed in the given temperature 
gradient, a redistribution of oxygen will occur from 
which the previously calculated O/M ratio and oxygen- 
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potential gradient are obtained. The calculation is 
repeated for a series of values of the CO,/CO or 
H,O/H, ratio until the desired range of values of the 
overall O/M ratio has been covered. Very similar results 
are obtained by assuming that either the CO,/CO or 
H,O/H, ratio controls the oxygen-potential gradient. 
The effect of changes in the power level and density 
can be studied by calculation of their effect on the 
temperature profile. The results reported by Rand and 
Markin? are shown in Fig. 6. These authors offer two 
significant conclusions; (1) an oxide with an overall 
O/M ratio of nearly 2 has a negligible oxygen gradient, 
and (2) oxides with appreciable deviations from 
stoichiometry have marked oxygen gradients; the 
deviation from stoichiometry decreases from the center 
of the fuel to small values at the surface. As seen from 
Fig.6, hyperstoichiometric (or hypostoichiometric) 
fuels would be expected to become more hyper- 
stoichiometric (or hypostoichiometric) toward the 
center when placed in a temperature gradient. 

Aitken et al.*° have used the results of their 
out-of-pile axial thermal-gradient experiments to com- 
pute the radial O/M ratio gradient in an oxide when 
placed in a radial temperature profile. The basis for 
their prediction consists of a series of experiments*’ in 
which 36 to 40 urania—plutonia pellets of the desired 
composition and stoichiometry were sealed in a 
12-in.-long molybdenum tube (0.25 in. in inside diam- 
eter). Pellet density was nominally 92% of theoretical. 
The sealed molybdenum tubes were placed in a 
molybdenum-wound furnace with a N,—6% H, 
atmosphere. The axial temperature gradient ranged 
from 900 to 1600°K. After the tubes were held in the 
gradient for 100 to 2000 hr, the tubes were opened, 
the pellets were removed, and the O/M ratio for each 
pellet was determined using the Lyon method.” A 
graph of the logarithm of the deviation from 
stoichiometry, x, vs. the reciprocal of the absolute 
temperature usually indicated a region with a straight- 
line segment. The slope has been identified by these 
workers as the effective heat of transport, Ort, and 
discussed in terms of the formalism of irreversible 
thermodynamics.*”** They found empirically that this 
effective heat of transport was roughly inversely 
proportional to the square of the average deviation 
from stoichiometry, i.e., Oo = —22.5/(¥)? kcal/mole. 
This relation holds for values of x greater than 
approximately 0,02. For average values of x less than 
or equal to 0.02, a constant heat of transport of 
approximately —30 kcal/mole was observed. These 
results were used to estimate the O/M gradient in a 
mixed oxide held in a radial temperature gradient. An 
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Fig.6 Oxygen distribution in U, ,, ,Pu,., ,O0,+,. as computed 
by Rand and Markin.? 


initial estimate of the temperature profile of the fuel 
was made for an assumed power level, theoretical 
density of the mixed oxide, and fuel-surface tem- 
perature. The equation of Laskiewicz'* was used for 
final temperature profile calculations. The fuel is 
divided into about 100 concentric-shell-volume ele- 
ments. A value for the O/M of the fuel at the oxide 
surface adjacent to the cladding is assumed. If the value 
of the deviation from stoichiometry is less than 0.02, a 
value of Q* = —30 kcal/mole is used and the value of x 
at the edge of the first cylindrical segment is cal- 
culated. This calculation is repeated stepwise until a 
value of x greater than 0.02 is reached. At this point a 
new value of Q* is chosen and the stepwise computa- 
tion is continued until either the center of the fuel or a 
value of the O/M ratio of 1.90 is reached. The average 
value of the deviation from stoichiometry of the oxide 
is then calculated, and the value of Q* used is 
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compared to the value of Q* computed from the 
average value of the deviation using the inverse square 
relation given above. If a significant difference is 
found, the assigned heat of transport is changed by a 
predetermined amount and the calculation repeated. 
Since the thermal conductivity depends on the O/M 
ratio, a new temperature profile was calculated each 
time a new heat of transport was chosen. In this 
manner an O/M ratio profile consistent with the given 
O/M ratio of the fuel at the surface and power level 
was obtained. No attempt has thus far been reported to 
include the effect of the fuel structure on the 
temperature gradient. With this approach it was con- 
cluded that the average O/M ratio of the oxide for a 
given deviation from stoichiometry at the fuel surface 
increased as the fuel-surface temperature was increased. 
For example, it was predicted that, for a deviation 
from stoichiometry of 0.002 at the fuel surface, the 
average O/M ratios of the oxide are 1.950 and 1.962 
for fuel-surface temperatures of 973 and 1373°K, 
respectively. The effect of the variation in linear power 
at constant surface temperature and deviation from 
stoichiometry was investigated. The average O/M ratio 
of the fuel was found to decrease with increase in 
power level, with the greatest variation being noted 
below 12 kW/ft; also, as the average O/M ratio 
approached 2, the oxygen gradient became small. 
Those results are similar to the ones reported by Rand 
and Markin.” 

The predictions made by Rand and Markin* and 
Aitken et al.*° are for a pure urania—plutonia solid 
solution of a given average or overall O/M in a radial 
temperature gradient. To apply these predictions to a 
fuel pin, i.e., an oxide clad with stainless steel, requires 
an estimate of the amount of oxygen that reacts with 
the cladding during the time that the fuel structure is 
established but before a significant amount of burnup 
has occurred. The amount of oxygen that reacts with 
the stainless-steel cladding would be expected to 
depend on the temperature, the oxygen potential, the 
condition of the steel surface, and the time of reaction. 
The amount of oxygen that must react with the 
cladding to produce a significant change in the average 
O/M ratio of the oxide is not large. For example, the 
formation of a 0.2-mil layer of Cr,03 on the inside 
radius of the stainless-steel cladding would reduce the 
average O/M ratio of the oxide from 2.00 to 1.99. 
Aitken et al.** suggest that, at a fuel-surface 
stoichiometry of 1.998 or greater, oxidation of the 
chromium in the cladding will occur, and therefore the 
use of a fuel with an average O/M ratio of 1.96 or less 
would not result in oxidation of the cladding in the 
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early stages of irradiation. Presumably one would 
predict that, if the initial O/M ratio of the fuel were 
greater than 1.96, oxidation of the cladding would 
occur until the average O/M ratio reached 1.96 and the 
O/M ratio of the fuel surface reached 1.998. Rand and 
Roberts™ do not discuss the oxidation of the cladding 
in detail but suggest that the cladding may act as an 


oxygen getter when the fuel has an initial O/M ratio of 


2.0. 

There have been no reports of measurements of the 
O/M ratio in the oxide as a function of position in fuel 
pins that have been irradiated to “zero” burnup, i.e., 
for a period long enough to establish the fuel structure 
but before a significant quantity of fission-product 
elements has accumulated. Therefore there is no way 
to compare these predictions with experiment. Most 
observations on irradiated fuel pins have been made 
after a significant amount of fission has occurred, and 
therefore the effect of burnup on the oxygen gradient 
must be considered. 


Effect of Burnup on O/M Ratio. Ina discussion of 


the effect of burnup on the O/M ratio, it is important 
to define the O/M ratio for a fuel in which some fission 
has occurred. In this review the O/M ratio is defined as 
the ratio of the total number of moles of oxygen to the 
sum of the number of moles of uranium, plutonium, 
and fission-product elements in solid solution in the 
fluorite phase, i.e., the urania—plutonia solid solution. 


Two other definitions have been used: (1) the ratio of 


the number of moles of oxygen not combined with 
fission-product elements to the total number of moles 
of uranium and plutonium, and (2) the ratio of total 


number of moles of oxygen to the total number of 


moles of oxidized elements in the fuel. Neither of these 
two definitions is useful in discussions of the change in 
the oxygen potential of the oxide with burnup. 

The effect of burnup on the O/M ratio has been 
considered by several investigators.°? *! Common to all 
these studies is the classification of the fission-product 
elements into two main classes, those which are 
oxidized and those which are not. In the oxidized class, 
there is general agreement that Sr, Y, Zr, Nb, Ba, and 
the rare-earth elements should be included. Likewise, 
Kr, Xe, Tc, Ru, Rh, Pd, Ag, Cd, In, and Sn are 
generally included in the nonoxidized class. The major 
differences between the various studies involve the 
chemical state of Mo, Cs (and Rb), Br, I, Te, Se, and 
Sb. It would appear likely that the relatively small 
yields of Br and I form alkali halides and do not 
significantly affect the oxygen balance in the fuel. 
Similarly, the oxidation state of the low-yield elements 
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Te, Se, and Sb will have little effect on the oxygen 
balance. The chemical state of Mo and Cs(Rb) is 
important to the oxygen balance in the fuel, and the 
extent that these elements are oxidized depends on the 
oxygen potential of the system; the extent of the 
oxidation of Mo and Cs(Rb) as a function of the 
average O/M ratio of the fuel remains to be ex- 
perimentally established. 

To compute the effect of burnup on the O/M ratio 
requires that the fission yields for the fission products 
which are oxidized be available. A computer program 
has been reported* which would permit the exact 
computation of the effect of various irradiation 
parameters on the fission-product elemental com- 
position if the basic fission-yield data and branching 
ratios for the various isotopes were exactly known. 
However, the uncertainties associated with the values 
of the elemental fission yields makes such computer 
calculation of limited usefulness, and uncertainties of 
5% or greater may be expected in any computation 
that depends on fission-yield data. 

As fission-product elements accumulate in the 
oxide matrix, there will be a tendency to form the 
phase or phases most stable under the local conditions 
of temperature and oxygen potential in the fuel. For 
example, as the concentration of alkaline-earth fission 
products (Sr and Ba) and zirconium increases, the 
tendency to form a separate phase of either SrZrO, or 
BaZrO, will increase. Eventually the activity solubility 
product for the formation of these ternary oxide 
phases will be exceeded. The chemical reactions may 
be written 


SrO(soln) + ZrO,(soln) = SrZrO,(c) 
BaO(soln) + ZrO,(soln) = BaZrO3(c) 


which lead to the activity-solubility-product relations 


4510 * 4710, = &xP [AGfs:2z10, 


— AGfgr0 — AGf2,0,)/RT| 


= aw a? ed 
4g 0 * 2z10, = &XP [(AGigazr0, 


— AGfga0 — AGf2,0,)/RT] 


where the a’s are activities and the AGf°’s are standard 
free energies of formation. The enthalpies of formation 
of the ternary oxides have been reported by L’vova and 
Feodosev * and the absolute entropies by King and 
Weller: therefore the standard free energies of forma- 
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tion may be estimated and the activity solubility 
products computed. For SrZrO 3 the solubility product 
has the values 3.4 x 10° and 1.5 x 10% at 2000 and 
3000°K, respectively, whereas for BaZrO,the corre- 
sponding values are 4.8 x 10 * and 5.3 x 10°. The 
mole fraction solubility product for the oxides depends 
on the fraction of fissions that occur in **°U and ?9°Pu 
and on the burnup. For 1 at.% burnup the mole 
fraction solubility product for SrZrO, has the values 
2.8 x 10? and 7.9 x 10° if all fission occurs in 7°°U 
and 73°Pu, respectively, whereas, the corresponding 
values for BaZrO, are 2.0 x 10? and 1.3 x 10. Thus, 
at 1 at.% burnup, the mole fraction solubility product 
exceeds the activity solubility product for the two 
phases, and, unless there are large negative deviations 
from ideal solution behavior, it is expected that both 
solid phases will tend to form early in the irradiation of 
a fuel pin. The actual formation of significant amounts 
of the solid phases will require nucleation and may be 
hindered by slow solid-state diffusion. Solid phases 
tentatively identified as alkaline-earth zirconates have 
been observed in irradiated fuel pins.*° Thus it is 
assumed that all the fission-product strontium and 
barium form zirconates that form a separate phase in 
the fuel matrix. This will remove a significant fraction 
of the zirconium from the fuel oxide as well as 3 moles 
of oxygen for each mole of alkaline-earth metal. 

The O/M ratio of the fuel-oxide phase after a 
fractional burnup, B, of the heavy atoms is given by 
the equation 


_(O/M); — KB 
er 1+ B(K, — 1) 
where (O/M), = initial O/M ratio of the oxide 
K,=the moles of oxygen used to form 
SrZrO, and BaZrO, (equal to three 
times the sum of the fractional yields 
of Ba and Sr) 
K,= the sum of the fractional yields of Y, 
rare-earth elements, and the Zr not 
used to form zirconates 


The rate of change in the O/M ratio with burnup may 
be obtained by differentiation: 


d(O/M) _ (O/M);( — K,) —Ky 
dB [1 +B(K, —1)]? 





From the data given in Table 1, values of K, and K, 
can be computed. In Table 3 are values for K,, Kz, and 
d(O/M)/dB for three different cases of fast fission: 
case I, where all fission is in 7°°U; case Il, where 75% of 
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Table 3. Effect of Burnup on O/M Ratio 





(O/M) Ratio 





Fast-fission Fast-fission Fast-fission 





case | case II case III 
K, 0.474 0.427, 0.288 
i 0.684 0.655, 0.570 
d(O/M)/dB 
For (O/M); 
of 2.00 
and B of 
0.1 0.17 0.28 0.62 
For (O/M); 
of 1.95 
and B of 
0.1 0.15 0.26 0.60 


Fission cases: I, all fissions in 774 U; II, 75% of fissions in 
235U and 25% in 73° Pu; III, all fission in 73° Pu. 





fission is in 7°°U and 25% in ?°°Pu (typical of EBR-II 
irradiations of enriched urania—plutonia test pins); and 
case III, where all fission is in 7°*Pu (typical of 
irradiations of normal or depleted urania—plutonia 
fuels in demonstration LMFBRs). There is a sig- 
nificantly greater increase in O/M ratio with burnup in 
the case where all the fission is in ?°?Pu compared to 
the cases where most of the fission is in ?°°U. It should 
also be noted that in all cases the O/M ratio increases 
with burnup. 

As the O/M ratio of the oxide increases with 
burnup, part or all of the fission-product molybdenum 
being generated will be oxidized. In the hotter regions 
of the fuel, this molybdenum will dissolve in the 
fuel-oxide lattice. This will reduce the rate of increase 
in the O/M ratio with burnup. For case II the oxidation 
of about 50% of the molybdenum produced would 
maintain the O/M ratio at a constant value. For case 
III, even if all the molybdenum produced were 
oxidized, the O/M ratio would still increase, but only 
about 0.001 O/M unit per percent burnup, or about 
one-sixth of the rate shown in Table 3 when no 
molybdenum is oxidized and dissolved in the fuel 
oxide. 

It seems unlikely that any significant fraction of 
cesium will dissolve in hypostoichiometric fuel. In the 
lower temperature regions of the fuel and when the 
oxygen potential has reached values near those char- 
acteristic of stoichiometric or hyperstoichiometric 
oxide, it is possible that cesium (and rubidium) may 
react either with the oxide fuel to form a separate 
cesium uranate phase or with molybdenum to form a 
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separate cesium molybdate phase. The overall effect on 
the O/M ratio of these two possible reactions may be 
estimated. For casell the formation of cesium 
(rubidium) uranate would reduce the rate of increase in 
the O/M ratio with burnup to about 0.0003 O/M unit 
per percent burnup, whereas for case III the O/M ratio 


would increase at a rate of about 0.004 O/M unit per 
percent burnup. For case II the reaction of about 57% 


of the cesium (rubidium) with molybdenum to form 
cesium molybdate would reduce the rate of increase in 
the O/M ratio with burnup to zero; for case III the 
O/M ratio would increase at a rate of about 0.001 O/M 
unit per percent burnup if ali the cesium (rubidium) 
produced reacted to form cesium molybdate. 

Thus it appears that, for stoichiometric and hyper- 


stoichiometric mixed oxide, for case Il (typical of 


EBR-II test-fuel-pin irradiations) conditions could exist 
such that the O/M ratio would not increase with 
burnup. However, for case III (typical! of irradiations in 
a demonstration LMFBR), there do not appear to be 
conditions under which the O/M ratio of the mixed 
oxide will not increase with burnup. In this latter case, 
however, the increase in O/M ratio with burnup may be 
significantly decreased if molybenum and cesium are 
oxidized. These conclusions assume that, when fuel 
conditions lead to oxidation of molybdenum and 
cesium, the molybdenum and cesium already present in 
the fuel due to prior fission do not oxidize. This 
assumption appears reasonable in the case of molyb- 
denum, since in no case have the noble metal inclusions 


observed in irradiated fuel been found to be free of 


metallic molybdenum. The assumption is not so well 
established in the case of cesium, since it is not possible 
to ascertain the chemical form of cesium in irradiated 


fuel owing to its great chemical reactivity. Studies of 


the stability of cesium uranate and cesium molybdate 


at temperatures and oxygen potentials characteristic of 


the cooler region of the fuel may enable one to 
estimate the extent that cesium is oxidized. When a 
fuel pin is considered, the possible uptake of oxygen 
by the cladding must be taken into account. Before the 
effect of burnup on the oxygen potential gradient ina 
fuel pin can be discussed, the effect of burnup on the 
oxygen potential must be considered. It should be 
noted that, although the O/M ratio of the fuel oxide 
and its oxygen potential (at a given temperature) are 
related, the relation depends on the cation composition 
of the oxide. This relation is considered in the next 
subsection. 

Effect of Burnup on the Oxygen Potential. [he 
oxygen potential, AGg,=RT In po,, of an oxide 
phase is a function of the temperature and com- 
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position. To fix the oxygen potential at a given 
temperature for the pure-urania—plutonia solid so- 
lution requires that two composition variables be 
specified; these are usually the O/M ratio, 
[O]/({U] + [Pu]), and the Pu/M ratio, 
[Pu] /({U] + [Pu]). After some burnup has taken 
place, the oxide will contain fission-product elements 
in solid solution. For each fission-product element in 
solid solution, an additional composition variable 
would be required. However, since for a given fission 
case, i.e., a given fraction of fission in 7*°U, *°8U, and 
39Pu, the fractional composition of the mixture of 
fission-product elements that dissolve in the oxide is 
essentially constant, a single composition variable may 
be used. It must be remembered that the relation 
between this composition variable and burnup will 
depend on the details of the fission process. Since 
during burnup the relative amounts of plutonium and 
uranium in the oxide may change, the plutonium 
composition variable will also change with burnup. The 
change in Pu/U ratio with burnup depends on the 
details of the fission process. Thus, for irradiations of 
enriched-urania—20% plutonia solid solutions in 
EBR-II, the Pu/U ratio does not change very much 
with burnup, since about 75% of the fission occurs in 
235 and 25% in ??Pu. On the other hand, when 
normal or depleted urania—plutonia solid solutions are 
irradiated, almost all of the fission occurs in **°Pu and 
only a very small fraction of the fission occurs in 7*°U. 
However, capture of neutrons by **°U leads, via **’Np, 
23°Pu lost by fission is 
238 In this 


to 739Pu, and so some of the 
replaced with an equivalent decrease in 
case the change in Pu/U ratio with burnup depends on 
the ratio of fissions in **°Pu to the capture reactions in 
238U. This latter ratio depends markedly on the 
neutron-energy spectrum of the reactor, which varies 
from reactor to reactor and with position within a 
given reactor. The neutron-energy spectrum in EBR-II 
is “harder” than expected for future LMFBR power 
reactors, and hence the capture reaction in 7°°U will be 
more important in the latter reactors. The details of 
these changes were discussed earlier in the section on 
“Fission-Product Yields.” 

The change in oxygen potential with burnup is 
illustrated in Fig. 7. The curves represent, at a given 
temperature, the dependence of the oxygen potential 
on the O/M ratio for (curve A) a pure-urania—plutonia 
solid solution, (curve B) a urania—plutonia solid 
solution with the same Pu/U ratio but containing 
fission products in solution, and (curve C) the same 
fission-product mixture as curve B but with a different 
Pu/U ratio. The effect of burnup is divided into three 
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steps for the purpose of the illustration, although in 
the actual process ali three changes occur at the same 
time. In the first step the O/M ratio of the fuel oxide 
changes from the initial value, (O/M);, to the value 
associated with point 2. This increase in O/M ratio 
causes an increase in the oxygen potential of the oxide. 











Fig. 7 Change in oxygen potential with burnup (diagram- 
matic), 


Next the fission-product elements that dissolve in the 
oxide cause the oxygen potential to increase to 
point 3. Finally, the Pu/U ratio, which is changed by 
unequal fission rates and capture, causes the oxygen 
potential to move to point 4. Curve C is based on the 
assumption that more uranium than plutonium was 
lost by fission; in most real fuel pins the opposite 
would occur and curve C would lie below curve B. 

To apply this analysis of the effect of burnup on 
oxygen potential requires that data for curves A, B, 
and C be available. Limited data are available to 
construct curve A, but no experimental data are 
available for curves B and C, The model for the oxygen 
pressures Over urania—plutonia solid solutions de- 
scribed in the section on “Oxygen Potentials in U—Pu 
Oxide Fuel” may be extended to multi-cation solid 
solutions. In principle, each cation may exist in one or 
more oxidation states in the solid solution. If more 
than one oxidation state is involved, mass-action 
expressions similar to those presented in the section 
just mentioned for plutonium and uranium ions must 
be derived. To include the effect of the fission-product 
elements, the cation sum and electrical neutrality 
conditions (see Eqs. 12, 13, 16, and 17) become 


[U] + [Pu] + [Fp] =1—(1—K,)B (18) 


and 
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[O?] = [Pu2"] + 1.5[Pu®] + 2[Pu*] + [U**] 
+2[U"] +3[U%] +4 (RF) (19) 


where Z is the mean charge of the fission-product 
cations and F’, is the cation fraction of fission products 
in the oxide solid solution. The [F,] is equal to KB as 
defined previously. This scheme has been used to 
estimate the effect of burnup on the oxygen potential 
of the urania—plutonia solid solution. An example of 
calculated curves for AGo, vs. O/M ratio as a function 
of burnup at 2000°K is shown in Fig. 8. These curves 
are for a natural-urania—20% plutonia fuel in which the 
ratio of fission in ?°?Pu and capture in 7°8U is such that 
the Pu/U ratio remains constant at 0.25 during burnup. 
The oxygen potential at a given O/M ratio increases 
with burnup. This is due to the substitution of cations 
for uranium and plutonium cations which have a lower 
mean charge. To maintain electrical neutrality requires 
that the mean valence of the plutonium and uranium 
cations increase; this leads to an increase in the oxygen 
potential. 

These calculated curves include the oxidation of 
molybdenum when the oxygen potential is large 
enough. It was assumed that molybdenum forms Mo” 
ions in the urania solid solution. The inclusion of 
molybdenum oxidation has only a small effect on these 
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Fig. 8 Change in AGo, with O/M ratio for natural-urania — 
plutonia mixture at 2000°K. 


EFFECTS OF OXYGEN CONCENTRATION ON FAST REACTOR FUELS 


curves. The path taken during the burnup of an oxide 
whose initial O/M is 1.95 and which does not lose 
oxygen is also shown in Fig. 8. It is seen that each 
change in AGo, due to a 2% burnup increment is 
composed of two parts: a change in AGo, due to the 
increase in O/M and a change due to the change in the 
cationic composition of the oxide. The calculated 
oxygen potential indicates that molybdenum oxidation 
begins between 4 and 6 at.% burnup. As pointed out 
previously, when molybdenum is oxidized and is added 
to the solid solution, the rate of change in O/M ratio 
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burnup is zero, the oxygen potential still will increase 
with burnup, owing to the replacement of uranium and 
plutonium cations by rare-earth cations whose charge is 
less than that of the actinide cations. The tendency of 
the fuel oxide to interact with the cladding depends on 
the oxygen potential, which, although related to the 
O/M ratio, also depends on the cationic composition of 
the oxide. Thus the change in both the O/M ratio and 
the cationic composition of the oxide with burnup 
must be specified to fix the change in the oxygen 
potential of the fuel with burnup. 
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Fig.9 Increase in AGO, with burnup for enriched- and natural-urania—plutonia mixtures at 2000° K. 


with burnup decreases. The part of the increase in 
oxygen potential due to increase in O/M ratio is 
therefore decreased when molybdenum oxidation 
begins. The total change in oxygen potential with 
burnup is smaller after molybdenum oxidation begins. 
The oxygen potential as a function of burnup for 
fission in natural-urania—plutonia solid solutions and in 
35U-enriched-urania—plutonia solutions is shown in 
Fig.9. It is seen that the oxygen potential increases 
more rapidly with burnup in the case of the 7°°U 
natural-urania—plutonia mixture. 
molybdenum begins at about 6 at.% burnup in the case 
of the enriched-urania—plutonia mixture. As pointed 
out earlier, the oxidation of about 53% of the 
molybdenum produced is sufficient to decrease the 
rate of change in O/M ratio with burnup to zero. 
However, even though the O/M ratio change with 


Oxidation of 


Oxygen Gradient in Irradiated Fuel Pin. The 
redistribution of oxygen which occurs when a fuel pin 
is irradiated depends on both thermodynamic and 
kinetic factors. For example, the increase in the 
oxygen potential with burnup will increase the thermo- 
dynamic driving force for the oxidation of molyb- 
denum; this in turn will increase the partial pressure of 
various gaseous molybdenum oxides and other volatile 
molybdenum compounds. The increase in pressure of 
these gaseous species will increase the rate of their 
transport in the thermal gradient in the fuel. The 
increase of the oxygen potential near the cladding will 
increase the rate of oxidation of the cladding. As 
oxygen is removed from the fuel oxide, the rate of 
increase in the oxygen potential due to burnup is 
reduced. The lack of knowledge of the interchange of 
oxygen between the oxide and the cladding makes it 
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difficult to predict the change in the oxygen-potential 
gradient with burnup for a fuel pin. Aitken et al.* 
predicted that, if the cladding retains its integrity, the 
fuel will eventually reach an average stoichiometric 
composition and that at this time no measurable 
stoichiometric gradient will exist in the fuel. The 
prediction may be compared with the results of a 
method** developed for the estimation of the oxygen- 
potential gradient in irradiated fuel pins. 

The basic idea of the method is the use of the 
Mo—MoO, system as an in situ redox indicator to 
determine the oxygen potential of the oxide as a 
function of position. Thus, at each point in the fuel 
where this method may be used, the following 
equilibrium is assumed to exist when the fuel pin is 
under irradiation: 


Mo (noble-metal alloy) + O,(g) 
= MoO, (solution in oxide) 


The metallic molybdenum is present as a constituent of 
the noble-metal alloy inclusions, and the MoO, is 
present as a solute in the urania—plutonia solid 
solution. The oxygen potential is given by the 
equation: 

IMo0O, 


AGo, = AGhitoo, + RT in 
MLO 


where AGo, =RT \n Po, 


Po, = oxygen partial pressure, atm 


T = temperature, °K 


AGfM oo, = Standard free energy of formation of 


MoO, at temperature 7, kcal/mole 
4M 00, = activity of MoO, in solid solution 


Ayo = activity of Mo in noble-metal-alloy in- 
clusions 


The use of the Mo—MoO, as a redox indicator for fuel 
pins has several advantages: 

1. The standard free energy of formation of MoO, 
is about 10 kcal more positive than the oxygen 
potential of stoichiometric urania—20% plutonia, and 
hence the equilibrium MoO,/Mo activity ratio for 
hypostoichiometric systems will not be too small to 
measure. 

2. The metallic (reduced) molybdenum is a con- 
stituent of the noble-metal-alloy inclusions, and the 
MoO, is a constituent of the oxide matrix. 
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Fig. 10 Oxygen potential vs. radius (for fuel pin F2R; initial 
O/M = 2; 6.5 at.% burnup at 16 kW/ft). 


3.The Mo—MoO, system is readily frozen into 
position when the fuel pin is cooled. 

4.The yield of fission-product molybdenum is 
relatively large. 

5. The chemical reactivity of molybdenum in the 
alloy inclusions and the oxide matrix is low and hence 
sample preparation is simplified. 


In general, the noble-metal-alloy inclusions are only 
found in the columnar-grain and equiaxial-grain regions 
of the fuel pin. These regions usually comprise half or 
more of the mass of the fuel. 

A typical result for the oxygen potential as a 
function of position obtained using this method is 
shown in Fig. 10, where the oxygen potential is shown 
as a function of position before irradiation and after 
irradiation. The before-irradiation curve was computed 
for a stoichiometric (O/M = 2) mixed urania--20% 
plutonia solid solution, when placed in the 
temperature-radius profile estimated for this fuel pin at 
a power level of 16 kW/ft. The before-irradiation curve 
assumes that no redistribution of oxygen has occurred. 
The after-irradiation curve was computed from the 
electron microprobe analyses of the molybdenum 
content of the noble metal inclusions and the oxide 
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matrix and the calculated temperature-radius profile. It 
is seen that the observed oxygen potential is more 
negative than that which is calculated for 
stoichiometric fuel. It appears that a significant 
quantity of oxygen has been transported out of the hot 
region of the fuel. After 6.5% burnup, the oxygen 
potential of the fuel in the central region of the fuel 
corresponds to hypostoichiometric fuel even though 
the oxide was stoichiometric prior to irradiation. 
Similar results were obtained for three other fuel pins 
that were studied in detail using the Mo—MoO, redox 
indicator method. In all cases the oxygen potentials in 
the hot central regions of the fuel correspond to oxide 
that is hypostoichiometric. For the four fuel pins, the 
power level ranged from 11.5 to 16 kW/ft, the burnup 
from 4.7 to 10.9%, the initial O/M from 1.988 to 


2.000, and the smear densities from 82 to 94% of 


theoretical. 

The oxygen-potential data shown in Fig. 10 may be 
used to estimate the O/M ratio of the oxide using the 
relation between oxygen potential, temperature, and 
O/M developed in the section on “Oxygen Potentials in 
U—Pu Oxide Fuel.” However, since the oxygen poten- 
tials pertain to a fuel pin after 6.5% burnup, the 
presence of fission products dissolved in the oxide 
matrix must be considered. In the subsection on 
“Oxygen Gradient in Irradiated Fuel Pin,” the 
appropriate calculational procedure was indicated. The 


results of the calculation of the radial dependence of 


the O/M are shown in Fig. 11. The calculation was 
made on the assumption, first, that the oxide was 
pure-urania—20% plutonia solid solution (open circles) 
and, second, that the oxide was a urania—plutonia 
solution of fission products (open diamonds). The 
difference between the two curves is about 0.015 O/M 
units. 

When the oxygen partial pressure is computed from 
the oxygen potential and the temperature for each 
point for the four fuel pins that have been investigated 
using the Mo—MoO, redox indicator method, the 
results shown in Fig. 12 are obtained. It is seen that the 
oxygen pressure is approximately equal at the same 
radius for the four fuel pins. Since the burnups were 
different for the four pins, this result suggests that the 
oxygen pressure remains constant at a given radial 
position over the irradiation period of all fuel pins. The 
oxygen-potential profile will therefore depend on the 
temperature profile, which is determined by the power 
level and the fuel structure. After the structure has 
been established, the oxygen-potential profile is 
expected to remain constant during the irradiation if 
the power level remains constant. When this tentative 
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Fig. 11 O/M ratio vs. radius calculated from oxygen poten- 
tial. 


conclusion is combined with the analysis of the effect 
of burnup on oxygen potential, it becomes apparent 
that a significant quantity of oxygen must be trans- 
ported out of the hot central region of the fuel during 
irradiation. This oxygen probably raises the oxygen 
potential of the fuel in the cooler regions so that 
various mixed-oxide compounds (Cs,;MoQ,, Cs,UO,) 
are formed and additional oxidation of the cladding 
OCCUIS. 

The results of these studies of the oxygen-potential 
gradient indicate that the oxygen potential of the fuel 
in the hot central region of the fuel probably does not 
increase with burnup, and therefore that the O/M ratio 
of the oxide in this region will actually decrease with 
burnup. These results are not in agreement with 


. . . 2.35 
previous predictions, 


which were made without 
consideration of the possible effects of mixed-oxide 
formation (in the cooler regions) or oxidation of the 
cladding. 

The main features of the sequence of events that 
takes place during the irradiation of a fuel pin may be 
summarized. During the initial period of irradiation, 


before a significant amount of fission has taken place, 
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the differentiation of the oxide into three principal 
regions takes place, i.e., separation into columnar, 
equiaxial and unrestructured regions, and central void 
formation. Concurrent with this structural differentia- 
tion, the major part of any redistribution of uranium 
and plutonium occurs. A parallel redistribution of 
oxygen occurs during this period: oxygen is trans- 
ported out of the hot central regions to the cooler 
outer regions where the tendency for cladding oxida- 
tion is increased. Thus, during this restructuring period, 
some oxidation of the cladding may occur depending 
on the initial O/M ratio of the fuel and the temperature 
profile, which in turn depends on the power level. If 
the power level is high, the temperatures will be 
relatively high in the center of the pin, which leads to 
relatively negative values of the oxygen potential. This 
in turn leads to relatively low values of the O/M ratio 
in the central region. The difference between the initial 
O/M of the oxide and the average O/M after the 
structure is established is a measure of the quantity of 
oxygen that is available during the restructuring period 
for cladding oxidation. For each power level, i.e., each 
temperature profile, there will be a characteristic O/M 
profile and average O/M value. This average O/M value 
will tend to decrease as the power level increases. 
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Fig. 12 Oxygen pressure vs. radius. 


REACTOR TECHNOLOGY, Vol. 15, No. 4, Winter 1972—1973 


C. E. Johnson et al. 


Hence, for each power level, there will be a value of the 
initial O/M ratio below which little or no tendency for 
cladding oxidation exists and above which the 
tendency for cladding attack increases. Therefore, at 
the end of the restructuring period, the fuel structure is 
established, a change in the Pu/U ratio has taken place, 
the oxygen-potential profile characteristic of the power 
level has been set up, and some oxidation of the 
cladding may have occurred. 

During the second period of the fuel-pin irradia- 
tion, burnup occurs. According to the hypothesis 
developed above, the oxygen-potential profile of the 
oxide does not change significantly during this period 
if the temperature profile does not change sig- 
nificantly. This leads to the transport of additional 
oxygen from the hot central region of the fuel to the 
cooler outer region. The buildup of fission products in 
the fuel leads to new processes for the transport of 
oxygen down the thermal gradient, and the ease with 
which oxygen may reach the cladding is increased. The 
formation of binary oxide compounds will tend to fix 
some of the oxygen, which might otherwise attack the 
cladding. The quantity of excess oxygen generated is 
directly proportional to the burnup. The _ pro- 
portionality constant depends on the details of the 
fission process, i.e., which isotopes are fissioned and, in 
the case of natural-urania—plutonia fuels, the ratio of 
fission in **°Pu to capture in 7°8U. 


Uranium and Plutonium 


Uranium and plutonium are observed to be redis- 
tributed within the irradiated mixed-oxide fuel with 
the plutonium concentration being enhanced in the 
center of the fuel. The melting point of the mixed- 
oxide fuel decreases with increasing plutonium content 
i.e., melting point {+2870 — 470[Pu/(U + Pu)]}°C. 
The increased plutonium concentration also raises the 
center temperature for a given neutron flux. Thus a 
reduction is required in the maximum permissible 
power level over that for unredistributed fuel to 
prevent center fuel melting. Sha, Huebotter, and Lo*® 
have calculated the effects of plutonium migration on 
fuel performance. These authors conclude that a 
penalty in the maximum allowable linear heat rate, 
resulting from a 50% increase in plutonium concen- 
tration, may be as large as 3.2 kW/ft. Hence, to achieve 
optimum fuel performance, it is important to under- 
stand the effects of all variables on plutonium redis- 
tribution so that they may be readily controlled. 

The redistribution of uranium and plutonium in 
the center of irradiated fuel pins has been the subject 
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of a number of investigations. In 1967 Beisswenger et 
al.47 proposed thermal migration of plutonium to 
explain the enhanced plutonium concentration ob- 
served at the hot end (2300°C) of a urania—plutonia 
sample held in a temperature gradient for 110 hr. Rand 
and Markin? suggested that fissile-element redistribu- 
tion could occur by vapor transport primarily by UO, 
vaporization. Rand and Markin also proposed that the 
oxygen moved down the gradient would be returned to 
the center of the fuel by CO, or H,O gas impurities. 
Because of a lack of precise thermal-conductivity data, 
Rand and Markin did not attempt to predict plutonium 
redistributions. They did calculate a pseudocongruent- 
vaporizing composition of MO; , for a urania—15% 
plutonia mixture. 

Chasanov and Fischer® studied the effects of 
thermal gradients on the distribution of plutonium in a 
U—Pu mixed oxide. These authors were unable to 
establish a reasonable correlation between their mea- 
surements and thermal-diffusion theory. There appears 
to be some agreement between the observed oxygen 
and plutonium redistribution and that based on a vapor 
transport model. In addition, Chasanov and Fischer 
show that, where melting and freezing of the fuel has 
occurred, a redistribution of plutonium can be pre- 
dicted from “normal freezing theory.” 

In 1969 Aitken and Evans” proposed, as Rand and 
Markin had, that the oxygen distribution in the fuel is 
set by a constant H,O0/H, or CO,/CO ratio. These 
authors suggested that the redistribution of uranium 
and plutonium would be established by a constant 
Pu/U ratio in the vapor above the condensed state. 

Olsen, Fitts, and Lackey,°° using the approach of 
Aitken and Evans, have calculated Pu/U ratios in the 
vapor phase as a function of temperature, plutonium 
and oxygen concentrations, and equilibrium H,O/H, 
ratios. In a qualitative fashion they predict that 
uranium and _ plutonium will redistribute by 
vaporization—condensation methods. The oxygen 
concentration or activity is also shown to be important 
in the extent of redistribution. 

Lackey, Homan, and Olsen®*' have developed 
computer models for predicting radial porosity and the 
distribution of oxygen, uranium, and plutonium in 
oxide fuel. The uranium, plutonium, and porosity 
distribution are based on migration of pores by vapor 
transport. These in turn depend on an oxygen distribu- 
tion calculated with the heats of transport of Evans, 
Aitken, and Craig.** Evans et al. obtained the heats of 
transport from out-of-pile thermal gradient experi- 
ments in which 0.06 atm of hydrogen was present. The 
presence of this amount of hydrogen could have 
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resulted in oxygen redistribution by the 
Rand—Roberts™ mechanism establishing a constant 
H,O/H, ratio in the gas phase (see the subsection on 
“Oxygen and Oxidized Fission Products’). On the 
basis of an irreversible thermodynamic treatment, it 
was concluded by Aitken et al. °” that the Rand— 
Roberts mechanism could only account for part of the 
distribution of oxygen and that solid-state diffusion of 
oxygen was the principal method for oxygen move- 
ment. Aitken et al. assumed that the Rand—Markin 
oxygen-potential data for hypostoichiometric fuel 
(MO...) can be represented as 


AGo, = -—1RT Inx +K 


with n = 2. However, for Up gPug ,0, the Rand—Markin 
data yield widely varying values of » which depend on 
x and T. The Rand—Markin data yield values of n from 
1 to 12 between O/M ratios of 1.998 and 1.90 and 
temperatures of 1000 to 2000°K. Hence the above 
relation with m = 2 is not valid for Markin and Rand’s 
oxygen potentials, and Aitken’s conclusions based on 
this assumption are untenable. 

Heats of transport for oxygen calculated using 
Blackburn’s oxygen potential model (see the section on 
“Oxygen Potentials in U—Pu Oxide Fuel”) and a 
constant H,0/H, ratio agree well with those measured 
by Evans et al.*”%* Thus the measured heats of 
transport appear to be controlled by a constant 
H,0/H, ratio with Pu, = 0.06 atm. Since the fuel is not 
likely to have a significant hydrogen pressure (probably 
no more than that in equilibrium with the sodium 
coolant, ie., ~10 atm), Evans’ measured heats of 
transport are not necessarily applicable to the fuel. 
Lackey, Homan, and Olsen also conclude that hydro- 
gen pressures will be too low in an irradiated fuel for 
oxygen redistribution, but these authors did not 
consider the possibility that Evans’ measurement might 
be controlled by hydrogen. 


Lackey, Homan, and Olsen’s use of Evans’ mea- 
sured heats of transport can also be faulted on the 
grounds that the redistribution of oxygen, instead of 
being set by the heat of transport, might also occur by 
vapor transport along with the redistribution of ura- 
nium and plutonium. However, the mechanism for 
oxygen redistribution is not fully understood and 
needs further investigation. Nevertheless, with their 
model, Lackey, Homan, and Olsen have successfully 
predicted the porosity and plutonium redistribution in 
a stoichiometric sol—gel fuel. One would be interested 
in seeing the method applied to other pins. 
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Battles et al.?* measured pressures of UO3, UO,, 
PuO,, and PuO over UO,-20 wt.% PuO,. In this study 
the partial molar enthalpies of vaporization for each 
gas were also measured, and oxygen pressures were 
calculated from the ratios of UO,/UO, and PuO,/PuO. 
There was poor agreement between the two sets of 
calculated oxygen pressures. However, Blackburn and 
Danielson” have since measured the voltage-dependent 
relative ionization cross sections of uranium and the 
uranium oxide molecules. The measured ionization- 
cross-section ratio of UO, to UO; was found to be 
about 13 (at the voltages used in the mass- 
spectrometric study) instead of an assumed ratio of 
about 1. This increases the UO, pressure and the 
atomic oxygen pressure calculated from the UO,/UO, 
ratio by a factor of ~13. Since ionization-cross-section 
data are not yet available for PuO and PuO,, it is not 
possible to reliably calculate the oxygen pressure from 
the PuO,/PuO ratio. Pressures of UO3, UO,, PuO, and 
PuO, vs. O/M ratio are plotted in Fig. 13. These 
pressures have been corrected with ionization cross 
sections measured for UO, and UO; and estimated for 
PuO and PuQ,. Battles et al.?? found a pseudocon- 
gruent vaporizing composition of about 
Up.7gPug,2201.93 at ~2240°K (the temperature was not 
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Fig. 13 Partial pressures over U, ,, Pu, ,, O,-, at 2240°K. 
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held constant but was varied between 1900 and 
2240°K). This sample had an initial composition of 
UO .sPup,202 and lost relatively more uranium than 
plutonium, and this behavior is in agreement with 
assumptions for vapor transport of uranium in the fuel. 

Ohse and Olson** have also measured the pressures 
of various species over U—Pu oxide, in this case 
restricting the investigation to the U-15% Pu oxide. 
These authors calibrated their quadrupole mass spec- 
trometer to establish “sensitivity factors” from curves 
of ion intensity vs. pressure. At an electron energy of 
13 V, Ohse and Olson gave a “sensitivity factor” for 
UO,/UO; of 8.8, compared with Blackburn and Daniel- 
son’s cross-section ratio of 18 at 13 V. Ohse and Olson 
do not consider fragmentation of molecules at 13 V 
ionization voltage. However, Blackburn and Danielson 
observed fragmentation of UO 3 at voltages as low as 
10 V, and at 13 V they found the UO} fragment ions 
to be ~1.5 times the UO} parent ions. Ohse and Olson 
also found a congruent vaporizing composition of 
Uo, ssPUo,1sO1,97 at 2000°K. 

The fact that different pseudocongruent vaporizing 
compositions were found in each study does not 
indicate an inconsistency between the two studies. A 
calculation of the composition at which the O/M ratio 
is the same in the gas and solid gives Uo, 7gPU9,220;.93 at 
2240°K and Up.gsPup, 1501.9, at 2000°K in good agree- 
ment with experimental values. Contrary to assump- 
tions found in the literature, the pseudocongruent 
vaporizing composition is not fixed but shows a 
decrease in the O/M ratio with increasing plutonium 
content and with increasing temperature. 

Both mass-spectrometric studies are in qualitative 
agreement with earlier calculations by Rand and 
Markin, indicating a high UO; vapor pressure for the 
oxygen-rich fuel. However, the measured UO; pressures 
are 3 to 13 times higher (at O/M = 1.97) than those 
calculated from a long extrapolation (1100 to 
~2000°C) of the oxygen-potential data of Markin and 
Mclver.? 

Adamson and Aitken®* reexamined the earlier 
assumptions of Aitken and Evans concerning a fixed 
U/Pu_ gas ratio.** In addition, they considered a 
pseudocongruent vaporization with respect to uranium 
and plutonium—i.e., the U/Pu ratios in the vapor 
phase are assumed to be the same as the U/Pu ratios of 
the solid phase. Both criteria failed to produce U—Pu 
distributions in agreement with those observed in 
irradiated fuels. With the poor agreement in oxygen 
pressures calculated from mass-spectrometric data- 
(based on calculations made before ionization cross 
sections were available), Adamson and Aitken prefer 
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the extrapolated Markin and Mclver oxygen potentials 
and an assumed ideal solution of PuO, in UO, to 
calculate pressures over the U—Pu oxide. 

Adamson and Aitken have developed a mathe- 
matical model for calculating the redistribution of 
uranium and plutonium in a mixed-oxide fuel. This 
mode assumes: 

1.The oxygen concentration is simulated by 
experimentally observed oxygen distributions in out- 
of-pile experiments.™ 

2.In addition to vapor transport down the tem- 
perature gradient, solid diffusion up the gradient is also 
believed to occur. 

3. Most of the redistribution occurs after formation 
of the central void. 

4. Local equilibrium exists between solid and gas. 


Recently Meyer et al.** have developed a model for 
uranium and plutonium redistribution originating from 
their studies on irradiated fuels. There are some 
similarities between this model and that of Aitken. 
However, Meyer does not include solid-state diffusion 
owing to its relative unimportance. Little if any 
actinide redistribution is observed in high-density fuels 
where vapor transport is restricted. 

In contrast to Aitken’s model, Meyer et al. suggest 
that most of the redistribution occurs during the 
restructuring and formation of the central void. Meyer 
also includes a variation in temperature to account for 
the change in thermal conductivity with increasing fuel 
density. In Meyer’s treatment the assumption that the 
partial pressures of the vapor components are linear 
functions of the condensed-phase concentration is 
questionable. The relation between concentration and 
pressure for a U-Pu oxide are too complex to 
state in simple equations. The data of Fig. 13 indicate 
the pressures are approximately exponential rather 
than linear functions of the oxygen concentration. 

Bober, Sari, and Schumacher®® have measured 
redistribution of uranium and plutonium in out-of-pile 
experiments. In their investigation, cylindrical sintered 
U-—15% Pu oxide pellets were heated in longitudinal 
gradients of 1400 to 2500°K for up to 5O0hr. The 
plutonium concentration increased from its initial 
value of 15% to as high as 25% at the center. Except 
for one experiment, they attributed the increased 
plutonium to preferential loss of uranium by vaporiza- 
tion. The exception was noted to be a case of thermal 
diffusion. 

In another investigation, Bober, Sari, and 
Schumacher®”’ drilled an axial hole in the U—15% Pu 
oxide cylinder to allow the vaporized oxides to 
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condense. In qualitative agreement with thermo- 
dynamic calculations, it was observed that the first 
vapor to condense was enriched in uranium. The ratio 
of Pu/U in the gas increased with time until it reached 
~15 mole % plutonium, corresponding to the overall 
plutonium content of the pellet. In general, they 
conclude that redistribution in fuels during irradiation 
is caused by short-term vaporization—condensation. 
For fuels with very high central temperatures 
(~2700°K), Bober et al. believe that, over long periods 
of time (i.e., <35 days), redistribution of plutonium 
will occur by thermal diffusion. 


ATTACK ON FUEL CLADDING 


Chemical Driving Force 


The intergranular penetration of fission products 
into the austenitic stainless-steel cladding of EBR-II 
irradiated mixed-oxide fuel pins was first observed 
several years ago by experimenters at Argonne National 
Laboratory (ANL).** Since that time, a number of 
irradiated fuels have been examined, and many show 
evidence of this phenomenon. The correlation between 
cladding attack, fuel fabrication, and reactor operating 
parameters is not yet clearly defined. However, the 
studies that have been made suggest that the principal 
parameters initiating and controlling the cladding 
attack are the oxidizing potential, the amount of 
available oxygen, and the temperature at the inner 
cladding wall. Selected fission products (e.g., cesium 
and molybdenum) certainly contribute to corrosion, 
and their transport to the cladding wall within the fuel 
is controlled by the temperature and local oxygen 
potentials. 

Two types of cladding attack have been observed at 
the fuel—cladding interface. One mode is a general 
recession of the cladding thickness by a uniform 
oxidation (matrix attack) of the stainless steel. The 
second mode is intergranular penetration by oxygen 
and fission products along grain boundaries in the 
cladding. Examination of high-burnup irradiated pins 
indicated that the matrix attack is generally limited to 
less than 12% of the cladding thickness, whereas 
intergranular corrosion, when it occurs, can be con- 
siderably deeper. 

Results are shown in Fig. 14 from numerous 
investigations? at ANL, and the results indicate a 
general trend toward deeper cladding degradation with 
increasing cladding inner-wall temperature. No cladding 
penetration was observed below 500°C. There appears 
to be little difference in the initiation of attack of 304 
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Fig. 14 Summary of ANL results*® on depth of penetration of 
fission products in 304L and 316 stainless-steel cladding on 
EBR-II irradiated mixed-oxide fuel elements. 


and 316 stainless steels in the temperature range 500 to 
540°C, but it should be noted that the data points 
(NUMEC-C) for 316 stainless steel are from fuel 
elements that had twice the burnup of those for 304 
stainless steel. Further, the trend of this penetration- 
vs.-temperature line may be much lower than indicated 
because of the bias of the SOV-1 data. This point is 
explained in detail below. The scatter of the data 
precluded making any conclusions regarding the effect 
of burnup or of the initial fuel O/M ratio on the 
amount of cladding attack. 

Perry et al. have described experiments that also 
indicate increased cladding attack with increasing 
temperature. For these experiments no visual pene- 
tration was observed below 580°C. Tests on fuel 
stoichiometry indicated increased cladding penetration 
as the O/M of the as-fabricated fuel was increased. 
Results from both thermal and fast-flux experiments 
having similar cladding temperatures showed similar 
depths of cladding attack. 

With cladding temperatures ranging to 750°C for 
the gas-cooled fast breeder reactor (GCBR), Fitts et 
al.” have reported depths of cladding penetration vs. 
temperature for various irradiated fuels with differing 
cladding materials. These workers indicate that 316 
stainless steel shows satisfactory performance for early 
GCBR applications (temperatures to 650°C) but 
indicate that more oxidation-resistant materials are 
needed for the expected higher temperatures to be 
found in the GCBR. 

Weber and Jensen™ have recently reported similar 
results for pins irradiated in a thermal reactor. Al- 
though it can be argued that fuel-cladding reactions 
will not be identical in thermal and fast reactors, the 
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observed trends are certainly applicable to fast reactor 
fuel elements. These investigations established that the 
amount of grain-boundary penetration of stainless steel 
was a function of the initial O/M ratio of the fuel. In a 
hypostoichiometric fuel with an initial O/M ratio of 
1.94, a layer of fission products (palladium, cesium, 
and tellurium) was deposited on the inner surface of 
the cladding, but there was virtually no cladding 
penetration at any temperature. (These results® are 
confirmed by recent examination of EBR-II irradiated 
NUMEC C-19 fuel, initial O/M = 1.94, 13% burnup, no 
cladding attack.) However, for similar fuels with initial 
O/M = 2.00, there were multiple layers of reaction 
products at the cladding inner surface and intergranular 
penetration of the cladding to a depth of 4 mils at 
surface temperatures above 535°C. 

The fission products cesium, molybdenum, 
tellurium, and iodine, which populate the fuel— 
cladding yap, are significant factors in influencing the 
degree and type of cladding attack. Cesium is definitely 
an active agent, whereas molybdenum, tellurium, and 
iodine may exert independent or simultaneous in- 
fluence with cesium on the mode and rate of attack. 
The effect of impurities (e.g., F, Cl, OH) on the 
fuel-cladding attack has not been fully evaluated. 
However, recent out-of-pile studies suggest that OH” 
appears to support severe intergranular attack. 

Theoretical calculations have indicated that 
hydroxide ion can be extremely important in cladding 
attack. The pressure of CsOH(g) may be calculated for 
the reaction 


Cs(l) + % O, + 4H, = CsOH(g) 


using the relation Poy = 4s ° Pb. Pi, 
exp { —AG, [CsOH(g)]/RT}. For present FFTF speci- 
fications the hydrogen pressure in the primary coolant 
is about 1.3 x 107 atm. When an inner-cladding-wall 
temperature of 700°C and a Cr,0, cladding corrosion 
product which gives an oxygen pressure o,) of 
2.68 x 10°*4 are assumed, a CsOH pressure of 
5.6 X 10° atmis calculated. At this pressure, sig- 
nificant transport of cesium and oxygen could occur 
throughout the irradiation life of the fuel. 

Oxygen analyses were performed™ on cladding 
sections taken from a number of irradiated fuels. Each 
cladding ring was thoroughly cleaned and divided into 
sections for oxygen analysis. The analytical data are 
given in Table 4 along with pertinent information on 
fuel properties and irradiation conditions. 

The extent of intergranular attack and the large 
oxygen concentrations found in the cladding sections 
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Table 4 Oxygen Content of Irradiated Stainless-Steel Cladding 
(UO, —20 wt.% PuO, Fuels Irradiated In EBR-II) 





Smear Peak Calc. inner- 
Sample density, Average linear cladding Cladding Oxygen 
Sample location,* Cladding % theo- burnup, power, temp., attack, content,+ 
designation in, type retical at.% kW/ft "€ mils ppM 





760 
SOV-1 (plenum) 127, 304 270 
440 


SOV-1 (45G-14) 104, 304 80 5.0 20 593 9 S00 
3000 


1800 


SOV-1 (45G-16) 57, 304 80 5.0 20 545 5 2100 
1100 


1600 
1300 
870 


800 
C-11 (plenum) 15%, 316L $20 


SOV-1 (45G-18) tf, 304 80 5.0 20 530 


nN 


C-11 (66A-22) 137, 316L 88.6 9.5 14 550 1-2 1100 
C-11 (66A-11) 8%, 316L 88.6 oS 14 530 0-1 
: C-15 (66B-17) je Ae 316L 82.4 9.5 13 550 1-2 440 


i C-15 (66B-7) Lys 316L 82.4 ps 13 525 0-1 


we 
_ 
_ 
i) 
o 


SOV-6 (43B-25) 10%, 304 83.5 Zt 17 560 2 470 
F2R-1 (GE) SE 347 94.3 6.5 16.5 550 None 940 


F2R-2 (GE) Sf, 347 94.3 6.5 16.5 550 None 1100 


007 (59B-20) ie 304L 76.5 4.5 AS.7 567 0-2 


t 
l 
1 
f 021 (S9D-3) 4/, 304L 80.2 4.5 16.4 565 3 240 
f 
ASOVG-17(69E-5) 54, r 304H 80.0 & IY 13.7 535 None 670 
ASOPC-3(69B-3) 57, - 304H 82.7 BY 14.8 530 1-2 630 


ASOPC-5(69C-7) 54 304H 77.1 3.0 12.1 540 1-2 1100 
8 


- oO OC ew VW 
nN 
io.2} 
So 
o 





*Location of cladding section as measured from bottom of fuel column. 
+ Analyses of irradiated stainless-steel cladding performed by Los Alamos Scientific Laboratory and ANL by an inert-gas fusion 
technique. 
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of fuel element SOV-1 set this fuel apart from all other 
specimens examined. These effects were probably 
caused by the fuel element having an excess amount of 
oxygen at the start of irradiation. The original O/M 
ratio of a fuel is very important and could be the most 
important of all fuel-fabrication specifications. Not 
only must the O/M ratio of the fuel be rigorously 
controlled, but small amounts of moisture must be 
excluded from the fuel during fabrication; excess 
oxygen from any source will result in oxidation of the 
cladding during the early stages of irradiation. Since 
the oxygen potential of the fuel at the fuel—cladding 
interface would be significantly greater than that of 
stainless steel in equilibrium with its oxide coating, 
conditions would exist for the potential oxidation or 
“activation” of the cladding, which would make it 
more susceptible to attack during irradiation. Excess 
oxygen appears to have existed in SOV-1 in that large 
quantities of chromium oxides were found at the 
fuel—cladding interface along the entire perimeter of 
the fuel. 

An increase in power level would also result in a 
steeper thermal gradient, and, under this condition, 
one would predict an increase in the oxygen potential 
for hypostoichiometric compositions at the fuel— 
cladding interface. Increases in power level will also 
result in higher temperatures at the fuel surface and at 
the inner-cladding wall. Increases in temperatures in 
these areas would directly affect the rate of the 
corrosive reaction. These factors probably had a 
significant effect on the cladding reactions of SOV-1, 
since the peak linear power rating of over 600 W/cm 
was the highest of all the fuels examined. 

The effect during irradiation of initial fuel density 
on the oxygen potential produced at the fuel—cladding 
interface is difficult to predict. At a given power level, 
low-density fuels, which have poorer thermal con- 
ductivities, tend to have higher centerline tempera- 
tures. These higher temperatures result in more 
extensive fuel restructuring and tend to vapor-transport 
more oxygen to the fuel—cladding interface. On the 
other hand, initially high-density fuel may hinder the 
transport of oxygen to the cladding surface. This effect 
can be seen in fission-product vapor species such as 
cesium and molybdenum oxide, which will pile up at 
the equiaxed-grain unrestructured boundary region in 
fuels with over 90% initial fuel density. 

Examination of the data in Table 4, excluding 
SOV-1, indicates that variations of the parameters in 
the ranges just discussed do not appear to have 
adversely affected the fuel lifetimes. These fuels had 
various starting densities and were irradiated at dif- 
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ferent power levels, but the original O/M ratios were 
similar. The most severe intergranular attack (for 
element 021) was 3 mils. These fuels also have burnups 
ranging from 2.7 to 9.5 at.%. Burnup is expected to 
increase the overall oxidation level of the fuel; how- 
ever, the oxygen available for reaction with the 
cladding depends on many factors. These factors have 
been fully discussed in the previous sections of this 
article. All of the fast reactor fuels examined were fully 
enriched and do not show the effects of extra oxygen 
that would be generated for fuels of all plutonium 
fission. If all the oxygen generated by burnup resulted 
only in an increase in the oxygen potential of the oxide 
matrix, conditions would be readily established in 
which oxidation of the cladding would occur at low 
burnups. This obviously does not happen since burnups 
of up to ~10 at.% have been reached without 
significant attack. 


Mechanisms 


In the austenitic stainless steels, such as those with 
18 wt.% Cr—8 wt.% Ni, resistance to corrosion is 
usually afforded by a protective-complex film of 
Cr,03, which may contain various amounts of Si, Mn, 
Mo, Fe, and Ni ions, depending on the alloy com- 
position and the corrosive environment. When the 
stainless steel is heated for short periods (5 to 30 min) 
at about 550 to 650°C, a process referred to as 
“sensitization” occurs. In this process the austenitic 
matrix is locally depleted of chromium by the pre- 
cipitation of Cr,,C, in the grain boundaries. These 
locally depleted regions are then more prone to 
oxidation, which promotes a severe intergranular 
attack. The cladding operating temperatures of 
LMFBR fuels are usually in the sensitization region, 
and the protective film of Cr,0;3 is probably rendered 
ineffective by thermal cycling and by reactions with 
such fission products as cesium. However, the actual 
mechanism of cladding attack produces some ob- 
served® variations in irradiated fuels. In many fuels, 
iron and chromium are transported out of the grain 
boundaries with nickel enhancement in the grain 
boundaries. In some fuels, iron and nickel are trans- 
ported out with chromium (probably Cr,,C,) enhance- 
ment in the grain boundaries. The attacked grain 
boundaries always contain cesium and rubidium with 
varying amounts of molybdenum, ranging from un- 
detectable quantities to more than the cesium content. 
Electron-microprobe images show oxygen concen- 
trating in the intergranular areas. Minor amounts of 
iodine are usually deep in the grain boundaries asso- 
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ciated with the cesium, but tellurium, although it is 
transported to the inner-cladding wall and is easily 
detected there, shows no tendency to penetrate deep 
into the intergranular areas. 

The intergranular-attack mechanism, we believe, is 
best characterized as stress corrosion or stress-assisted 
corrosion with metal physics, microstructure, and 
electrochemical features combining to produce accel- 
erated penetration of the stainless-steel cladding. The 
attack can produce either an intergranular or a trans- 
granular path in stainless steel. 

There is wide acceptance®™ that the phenomenon of 
stress corrosion is basically an electrochemical one 
whose rate is determined by the speed with which 
surface features are influenced by moving dislocations 
and surface adsorption processes. 

The irradiated LMFBR fuel pin contains a very 
complex gas—liquid—solid environment that can 
change with operating power density and fuel oxygen 
potential. This in turn affects the fission-product mix 
at the cladding. Molybdenum is especially sensitive to 
these parameters. A feature of stress corrosion is the 
pronounced specificity of the environment, and only 
certain metal—electrolyte combinations produce the 
intergranular attack. 

Thus a galvanic cell that could produce anodic 
dissolution at the crack tip and in the grain boundaries 
is as follows: 


Cr + 4,07 = ¥%Cr,0; + 3e 


Fe + O* = FeO + 2e 


The cathodic reaction would involve the reduction of 


some Oxygen carriers, e.g., 
Cs,MoO, + 2e = Cs,0 + MoO, + O? 
The overall cell reaction would be 
Cr + ¥%Cs,MoO, = 4 Cr,03 + %Cs,0 + 4,Mo0, 


For the local galvanic cell to be effective, a thin 
film of electrolyte would have to be condensed in the 
reaction zone area (grain boundary) to provide ready 
transport of ionic species. The electrolyte has not yet 
been identified, but if it exists it would probably be an 
oxide mixture of fission products and cladding com- 
ponents. Relatively pure cesium oxide would not be 
thermodynamically stable at the oxygen potentials 
expected; however, the formation of such compounds 
as Cs,MoO, would stabilize the cesium and oxygen 


activities in this environment. In a study of the 
Cs,0—MoO, system, Spitsyn and Kuleshan®’ found a 
temperature minimum of 450°C for an equimolar 
mixture of Cs;MoO, and MoQ3. This temperature is 
well below that expected near the fuel—cladding 
interface. Intergranular attack has not been reported to 
occur below 500°C; a possible explanation might be 
that the electrolyte on the cladding surface was solid 
below 500°C and therefore that an electrochemical 
reacticn could not take place at a reasonable rate. 


Postirradiation Mechanical-Property Tests 


The intergranular penetration of stainless steel by 
fission-product cladding is of great importance; also of 
concern is its effect on mechanical properties of the 
cladding. Complete understanding has yet to be 
achieved in that cladding tested after irradiation has 
not exhibited intergranular fission-product penetration. 

Pacific Northwest Laboratory (PNL) has con- 
ducted® several cladding burst tests on UO,—25 wt.% 
PuO, fuel-element cladding (304 stainless steel) after 
irradiation in EBR-II. These elements were irradiated in 
EBR-II to ~4 at.% burnup at 10.5 kW/ft peak linear 
power. Testing was done by internally pressurizing the 
specimens at a rate of 20 psi/sec after thermal 
equilibrium had been established at 900°F. The results 
of these tests are summarized in Table S. 

The irradiated burst specimens are seen to have 
higher hoop strengths and lower ductilities than the 
unirradiated controls tested at the same temperature. 
Such results are consistent with the tensile properties 
of austenitic stainless-steel structural material tested 
after irradiation in EBR-II. The stress-rupture life of 
the irradiated cladding specimens decreased with 
applied hoop stress at a test temperature of 1200°F, 
again an expected result. General Electric Company 
(GE) also has conducted biaxial burst tests on cladding 
sections from irradiated specimens with results some- 
what similar to those of PNL. Again, metallographic 
examination showed no intergranular penetration of 
the specimens. Thus to date there is a paucity of 
mechanical-property-test data on fuel cladding that is 
known to have experienced intergranular penetration 
by fission products. 


Out-of-Pile Corrosion Tests 


Out-of-pile experiments have been run by a number 
of experimenters to study the effects of fission 
products (Cs, Mo, Te, and 1) and other potential 
impurities (OH and Cl) on austenitic stainless steel. 
The object of these studies is to establish the reactants 
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Table 5 PNL-2 Biaxial Test Results 








Fluence, 
Specimen Av. cladding- 107? neutrons/cm? Test Max. hoop Rupture Strain 
identification ID temp., °F (E > 0.1 MeV) temp., °F stress, psi life (AD/D) , % 

Control (unirrad.) 1000 54,000 Burst 15—20 
PNL 2-3:K 935 pa) 1000 70,600 Burst 5.0 
Control (unirrad.) 1200 45,000 Burst 15—20 
PNL 2-6:J 930 23 1200 57,100 Burst 6.4 
PNL 2-6:K 930 1.9 1200 25,760 108 hr 3.0 
PNL 2-6:E 870 2.0 1200 19,500 870 hr N.A. 





and reaction conditions to simulate the type of 
corrosive attack observed in irradiated fuels. 

Stainless steels have been tested for compatibility 
with metallic cesium. Results” tend to show that 
stainless steels are relatively unaffected by pure cesium 
liquid and vapor to temperatures as high as 870°C. 
Another investigator” indicated that surface wetting of 
stainless steel by liquid sodium could lead to embrittle- 
ment and lowering of the tensile properties. Liquid 
cesium wetting of 302 stainless-steel tensile specimens 
decreased the room-temperature yield strength from 
73,000 to 6400 psi and lowered the ultimate tensile 
strength from 105,000 to 81,000 psi. Elongations were 
somewhat reduced in the cesium-wetted tests. The 
tensile properties appeared to recover when tests were 
conducted after removal of the cesium from other 
wetted specimens. Several other iron- and nickel-based 
materials behaved similarly, but the results were not 
statistically significant. 

Rubin and coworkers at GE Vallecitos” have been 
conducting a series of out-of-pile tests on the migration 
of oxygen and cesium within a temperature gradient. 
UO, pellets of O/M= 1.97, an amount of cesium (in 
the form of Cs,MoO, or Cs,O) corresponding to ~2 
at.% burnup, and a 316 stainless-steel-cladding material 
were placed inside either a molybdenum or 316 
stainless-steel capsule. A 550 to 1500°C (1022 to 
2732°F) temperature gradient was established over the 
length of the 6-in. capsule. Oxygen migration in the 
presence of cesium was found to be much greater than 
in similar tests without cesium. A hot-end temperature 
of ~1500°C previously had been found necessary to 
effect oxygen migration in the absence of cesium. 

When cesium was added as Cs,MoO, with excess 
uranium (to decompose to UO,+Cs,O+Mo at 
900°C), molybdenum was also found to migrate to the 
cool end of the capsule. This observation agrees with 
GE and ANL microprobe results, which show both 
cesium and molybdenum to concentrate at the fuel— 
cladding interface. In tests using Cs,MoO, as the 
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cesium additive, GE found that intergranular attack 
occurs over the temperature range 650 to 750°C. 

Out-of-pile tests conducted at ANL by Battles and 
Johnson™ have sought to establish the reactant and 
reaction conditions that will simulate the type of 
corrosion attack observed in irradiated fuel elements. 
This experimental program has followed two lines of 
attack. Differential thermal analysis has been used to 
examine compositions in the Cs,0—Cs,MoO,—MoO, 
system to provide information on temperatures and 
compositions of potential electrolyte systems. A 
second approach has utilized capsule tests to determine 
the specific conditions under which different types of 
cladding corrosion occur. Examination of all containers 
was carried out using metallographic and microprobe 
techniques. 

The types of corrosive attack observed were classi- 
fied into three distinct categories: 

1. Shallow intergranular attack. 

2. Uniform oxidation. 

3. Intergranular attack. 


Shallow intergranular attack was the predominant 
attack mechanism for capsule tests in which 
stoichiometric Cs,0 was the principal sample material 
present. The attack was probably initiated along grain 
boundaries as the cesium oxide melted, but, as reduc- 
tion of cesium oxide occurred, attack by the liquid 
metal appears to have centered on the triple point of 
grain boundaries. The oxygen pressure was estimated 
to be very low, ~107° atm, in these experiments. 
Uniform oxidation was found to occur in experiments 
containing Cs,0,, Cs,03, and MoO. In these experi- 
ments a uniform layer of oxide product (Fe, Cr, 
oxides) adhered to the interior wall of each capsule. 
Results suggested that the liquid electrolyte which 
could have been present early was later consumed 
by interaction with the stainless steel or interaction 
of sample materials [e.g., Cs,O(1) + MoO;(s) = 
Cs,MoO,(s)]|. Intergranular attack was noted in experi- 
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ments having Cs,0, and Cs,0; present. It is very 
evident that liquid electrolyte was present throughout 
all experiments, and, with the high oxygen pressure 
(~10% atm), accelerated attack took place. Ex- 
periments were also run in which the additives CsCl 
and CsOH were added to the electrolyte systems; CsCl 
was found to have little influence on the extent of 
reaction even at the highest oxygen pressure. Addition 
of CsOH to electrolyte systems of high oxygen pressure 
resulted in severe intergranular attack. The combina- 
tion of CsOH and CsCl proved to be the electrolyte 
system that gave the most extensive attack with 
rupture of the capsule occurring. 

Maiya and coworkers’ conducted laboratory 
studies on the kinetics of the reaction of cesium oxide 
with 304, 304L, 304H, 316, and 318 stainless steels. In 
all cases, severe intergranular attack of the stainless 
steel was observed. It is believed that these laboratory 
studies were run at too high an oxidizing potential and 
are not representative of what is found in irradiated 
fuel systems. 

Recently, G6tzmann and Hofmann” have reported 
on a cesium-accelerated oxidation of stainless-steel 
cladding. They subjected a piece of stainless steel to 
the following environments at 800°C for 1000 hr: 

1. UO, og + stainless steel. 

2. UO, og + Cs stainless steel. 

3. UO, o9 + Cs stainless steel. 


No attack was noted in experiments 1 and 3, whereas 
severe intergranular attack was noted for experiment 2. 
This is a very clear influence of the cesium in a 
high-oxygen environment and how it accelerates attack 
along grain boundaries. Again, as with the experiments 
of Battles and Johnson, it is almost certain that a liquid 
cesium oxide electrolyte was formed in experiment 2 
which facilitated the intergranular attack. 

Intergranular penetration of stainless-steel cladding 
during irradiation is as yet not well enough understood 
to make an adequate prediction of depth of 
penetration as a function of such parameters as 
cladding temperature, oxidizing potential (O/M ratio), 
and linear power rating. There are, however, certain 
trends observed in the out-of-pile tests. Intergranular 
penetration of stainless steel by fission product cesium 
has not been observed at temperatures below ~500°C, 
and the degree of attack increases with temperature. 
The intergranular attack appears to be the result of an 
oxidation reaction with cesium, oxygen, and some- 
times molybdenum participating. These conditions also 
enable the formation of an electrolyte, and when the 
appropriate, but as yet undefined, oxygen potential is 
established, corrosive attack can occur. 
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SODIUM—FUEL REACTION 


The chemical interaction of the sodium coolant 
with the U—Pu oxide fuel in the event of any breach in 
the stainless-steel cladding can cause fuel swelling and 
cladding rupture. The severity of the problem may 
determine how long failed fuel elements can be left in 
the reactor core and also whether or not failed fuel 
elements can be stored in the reactor sodium coolant. 


Reaction Product 


Although no reaction occurs between stoichio- 
metric UO, and oxygen-free sodium at normal reactor 
temperatures, the presence of excess oxygen, i.e., as 
UO,4, or oxygen dissolved in the sodium, leads to the 
formation of a voluminous compound of Na,UO, 
(density = 5.6 g/cm*). A number of other products 
(i.e., Na,0, NaUO,, Na,UO,) have been proposed as 
the reaction product, but only Na,UO, can be in 
equilibrium with sodium and UQ). 

Pepper, Stubbles, and Tottle”’ were the first to 
identify Na3;UO, by chemical and mass-balance analysis 
as the phase which can exist in equilibrium with 
sodium and UO,. These authors also established the 
structure and lattice parameter of Na3UQ, as face- 
centered cubic (fcc) with a = 4.79 to 4.80 A. Bartram 
and Fryxell” found the structure of this compound to 
be more complex than that proposed by Pepper, 
Stubbles, and Tottle.” In addition to the strong lines 
found by Pepper, Stubbles, and Tottles, many addi- 
tional lines were found by Bartram and Fryxell, who 
indexed them as a cubic cell with the space group 
P4,3,. On the basis of the structure alone, Bartram and 
Fryxell proposed the compound to have the formula 
Na,,U;0,,. However, because of the presence of other 
phases, they were unable to obtain confirmation for 
this composition by chemical analysis. 

O’Hare et al.”* obtained X-ray data for a reaction 
product which agree with those of Bartram and 
Fryxell, but chemical analysis of the compound shows 
it to have the composition Na3,UQ,. Furthermore, 
O’Hare et al. found that this compound must have a 
fairly narrow composition range because slight varia- 
tions in overall concentration always produced at least 
one other phase. Blackburn et al.*' established the 
phase relations among UO,, Na, Na;UO,4, Na,UO,, 
Na,U,0,, and Na,UO, for the Na—U—O system which 
showed, in agreement with Pepper, Stubbles, and 
Tottle, that Na,UO, is the phase in equilibrium with 
UO,. The phase relations are shown in Fig. 15. 

Martin and Mrazek”*® have also equilibrated 
Up, gPtp,.202., with sodium to find by X-ray analysis that 
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the reaction product for the mixed oxide is 
Na3U;.)PuyO,(Na3;MO,). The compound is _iso- 
morphous with Na,UO, and has a Pu/U ratio which is 
nearly the same as the U—Pu oxide from which it is 
formed. X-ray analysis of Na;MO, showed it to have a 
face-centered cubic structure with a lattice parameter 


of 4.77 A. 
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Fig. 15 Isothermal section of the Na—U—O phase diagram at 
about 800°C. 


Equilibrium Conditions 


The reaction of sodium and fuel results in three 
phases, sodium, fuel, and Na;MO,, in local equilibrium. 
The important variable is the oxygen pressure in 
equilibrium with these three phases. If both the oxygen 
pressure produced by the fuel and that produced by 
oxygen as dissolved in sodium are less than the pressure 
in equilibrium with the three phases, no reaction will 
occur. Oxygen pressures greater than the equilibrium 
value will cause the reaction to proceed until the 
oxygen pressure is reduced to the equilibrium value. 
The equilibrium oxygen pressure is assumed to be 
nearly the same for reaction of sodium with either UO, 
or (U—Pu)O,. Equilibration reactions of sodium and 
fuel support this assumption. 

Recently, measurements have been made to acquire 
thermodynamic data for Na;UO,. These data include 
enthalpy of formation by two methods: mass- 
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spectrometric measurements of sodium and oxygen 
pressures over the Na—U-O system® and solution 
calorimetric experiments on Na,UO, (Ref. 79). Heat- 
capacity data for Na,;UO, were obtained from 5 to 
350°K with an adiabatic calorimeter.® In addition, the 
enthalpy of Na,;UO, was measured between 500 and 
1200°K by drop calorimetry. The thermodynamic 
properties of Na,UO, at 298.15°K are summarized in 
Table 6. (Note the excellent agreement for the 
enthalpy of formation of Naz;UO, by two independent 
methods.) 


Table 6 Thermodynamic Properties of 
Na3 UO, at 298°K 





Cp 41.35 cal mole/°K 

Ss 47.37 cal mole/°K 

AH} —477.7 kcal/mole (solution calorimetry) 
AHf —477.8 kcal/mole! (mass spectrometry) 





These data, together with auxiliary data, are used to 
calculate the equilibrium oxygen pressure over 
Na—UO,—Na,UO, for the reaction 


Na,UO, = 0, + 3Na(I) + UO,(s) (20) 


The calculated oxygen pressures are those in 
equilibrium with stoichiometric UO, confirming the 
observation that only UO,,, reacts with sodium. To 
calculate oxygen pressures for Na3;MO, and fuel, 
assume that the difference in the free energy of 
formation between Na3MO, and (U,Pu)O, is the same 
as that between Na,UO, and UOQ,. A correction is also 
made for the oxygen deficiency in the fuel. That is, the 
oxygen pressure is calculated for the reaction 


ZT x 
Na3U.)PuyO, = pe’ ae 0, 3 3Na(I) 


+ Uy Pu, O2.. (21) 


The oxygen potential for Eq.21 may be calculated 
from 


2 : * 
AGo, = Fey| 46,1%a,U09) - ( a 


+ F)AGAUO,) -= AGPu0,) — sa6y,| (22 


where the AG;’s are the free energies of formation and 
AG’s are the oxygen and sodium potentials (partial 
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molar free energies). The sodium potential in the case 
of sodium—fuel contact is zero. If there is a tempera- 
ture gap, the sodium potential will be negative. 

The calculated equilibrium oxygen content of the 
fuel is uncertain because experimental measurements 
of oxygen-potential data for the fuel are in poor 
agreement (see the section on “Oxygen Potentials in 
U—Pu Oxide Fuel’). 

A recently developed oxygen-potential thermo- 
dynamic model for the fuel (see the section on 
“Oxygen Potentials in U-—Pu Oxide Fuel’) gives 
oxygen potentials that are close to an average of those 
measured experimentally. Calculations of the oxygen 
content of uranium oxide and mixed oxide (U—Pu) 
fuels based on these oxygen potentials, and with 
thermodynamic data for Na;UO, considered to be a 
close approximation to those for Na;MOu, give results 
in excellent agreement with experimental equilibration 
measurements for mixed fuel oxides containing 10 and 
20 wt.% plutonium (see Table 7 and below). An 
estimated uncertainty in the oxygen potential of +3 
kcal results in O/M uncertainties from +0.004 (10 wt.% 
plutonium at 1200°K) to +0.02 (30 wt.% plutonium at 
600°K). 

At a given temperature the equilibrium oxygen 
content of the sodium for fuel—sodium--Na,;MO, 
equilibrium is practically constant for either UO, or 
U—Pu oxide fuel. However, as shown in Table 7, the 
equilibrium oxygen content of the sodium is strongly 
dependent on the temperature at the sodium—fuel 
interface. These calculations are based on solubility 
data for Na,O in sodium,® Henry’s law, and oxygen 
pressures from Eq.21. The calculated oxygen con- 
centration of sodium is completely independent of the 
fuel oxygen potentials so that uncertainties in the fuel 
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oxygen potential will not affect the sodium oxygen 
concentration. 

The data in Table 6 refer to the conditions 
prevailing at the sodium—fuel interface. That is, for the 
temperature at the interface and for a fuel of given 
plutonium content, the O/M ratio of an oxygen-rich 
fuel will be decreased to the level indicated. Other 
calculations, based on densities of the fuel and the 
product, indicate a volume increase of about 1.0% for 
each 0.01 O/M decrease in the fuel oxygen concen- 
tration as the fuel reacts to form Na,;MO,. One of the 
problems still to be resolved is the distribution of 
oxygen expected in the rest of the fuel away from the 
sodium—fuel interface and toward higher temperatures. 
Further work is required to establish the oxygen 
distribution in fuel under normal operating conditions. 
Work at ANL, involving measurements using irradiated 
fuels, has not yet included fuel with an O/M ratio as 
low as that expected after reaction with sodium (see 
the section on “Attack on Fuel Cladding’). At GE, 
oxygen distribution in out-of-pile thermal-gradient 
experiments is being measured with and without 
sodium, *"%6 


Isothermal Study 


Sodium has been reacted with fuel under iso- 
thermal conditions” *®"*’ and with the fuel in a 
thermal gradient.** An example of the reaction be- 
tween sodium and Up,sPuo,.0, pellets at 1173°K is 
shown in Fig. 16. The darker material is Na;MQg, 
which forms a surface layer and surrounds the lighter 
oxide fuel grains. The dark circles are fuel porosity. 
The growth pattern appears to consist of the formation 
of a surface layer and also grain-boundary penetration. 


Table 7 Calculated and Experimental Isothermal Oxygen Contents of 
Fuel and Sodium for Fuel—Sodium Reaction 











Temp., Temp., Calc. O/M for O/M for U, , Pu, ., 0, O/M for U, ,, Pu,., O, O/M for UO, ppM O, in Na 
°F °K U,,., Fez.29; Calc. Meas.* Calc. Meas.* Calc. Meas.* Calc Meas.+ 
621 600 1:93, = 002 1.95, 2 0.055 1.98, + 0.007 2.00 0.01 
801 700 £393. 5.396,, 1.98, 2.00 0.20 
981 800 1.94, 1.96, (1.95, )t 1.98, 2.00 1.6 

1161 900 1.94, 1.96, 1.98, (1.98)¢ 2.00 2.00 8.0 

1341 1000 1.94, 1.96, 1.98, (1.98)¢ 2.00 2.00 30 0.2 
1521 1100 1.94, b26- j Pe oF 1.98, 1.98 2.00 2.00 95 0.3 
1701 1200 1.95, + 0.01 1.96, + 0.007 Eo, 1.98, + 0.004 1.98 2.00 2.00 170 0.5 





*Mastin.’?*?* 
+D. L. Smith (private communication). 


{May be kinetically controlled. lor an operating fuel pin: ag = f (radius) because of large AT. 
2 
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Fig. 16 Na,MO, layer formed on surface of mixed-oxide pellet 
exposed to sodium at 900°C for 2 days (as polished; magnifica- 
tion, 300 X before a 60% reduction for printing). 


The initial grain-boundary penetration proceeds by 
formation of a very thin intergranular layer (as 
indicated by the thin lines). This is followed by 
thickening of the Na,UO, layer, leading to cracking 
and further intergranular attack. In extreme cases the 
oxide crumbles into a powder. 

An objective of the isothermal study was to 
establish the equilibrium O/M value of the mixed oxide 
as a function of temperature and of the Pu/(U—Pu) 
ratio of the mixed oxide. A lattice parameter method 
was used. The reaction products were treated with 
ethyl alcohol to remove the free sodium and were then 
examined by X-ray diffraction. The lattice parameter 
of the mixed oxide was used to establish the O/M 
composition of the mixed-oxide phase. The measured 
O/M results are given in Table 7. Data from the 
thermal-gradient experiments®® are in general agree- 
ment with those obtained in the isothermal experi- 
ments. That is, U—20 wt.% Pu oxide is reduced to an 
O/M ratio of slightly less than 1.96 for U—20 wt.% Pu 
oxide. 
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Aitken has discussed out-of-pile and in-pile ex- 
periments on the sodium—fuel reaction.®* He assumes 
that the sodium oxygen concentration threshold for 
reaction is 10 to 100 ppM. This corresponds to 
calculated concentrations for the higher fuel—sodium 
temperatures in Table 7. The measured oxygen con- 
centrations are much lower than those calculated. 
Because of the assumed high equilibrium oxygen 
concentrations in the sodium, he concludes that the 
sodium—fuel reaction will be limited or even reversed 
by sodium removal of oxygen. However, since the 
oxygen content specifications for sodium are between 
1 and SppM, only those reactions in which the 
sodium—fuel interface is in excess of 600°C can involve 
oxygen removal by sodium (see Table 7). The 
measured equilibrium sodium concentrations are below 
specifications for sodium so that no oxygen removal 
could occur. The removal rate, if it occurs, must be 
faster than the rate of reaction of sodium with fuel to 
keep the reaction from reaching the maximum possible 
destruction of the fuel pin. 

To establish whether or not the reaction can be 
limited by sodium carrying away oxygen, one needs 
data for the oxygen-removal rate and for the sodium— 
fuel reaction rate. No data exist for the former. The 
reaction rate of sodium with fuel may also be 
important in establishing how long failed elements may 
be left in the reactor. 

For reaction rates of sodium with fuel, the data are 
still incomplete. Martin®’ has shown that the reaction 
rate is temperature dependent. For fine powder the 
isothermal reaction reaches completion in <2 days at 
1173°K, <7 days at 1073°K, and >8 days at 873°K. 

Some in-pile reaction rates of defected pins with 
sodium appear to vary as the fourth power of the 
time.8® The mechanism for the reaction is proposed to 
be one in which the rate is controlled by sodium 
diffusion from a point defect through the product 
layer. The product is assumed to form a hemisphere 
with its center at the point defect. However, the 
defected pins showing a one-fourth power time de- 
pendence had slit defects rather than point defects, and 
even an initial point defect will form a larger opening 
in time. Furthermore, when an equation was being 
derived for the rate of increase of fuel-pin diameter 
with time, an incorrect expression was used for an 
assumed sodium-diffusion-controlled reaction. The 
increase in mass of product with time is assumed to be 
inversely proportional to the radius of the product 
hemisphere growing from the pinhole cladding defect 
(dm/dt =k/r). This leads to cladding swelling, which 
varies as t4. However, the assumed mass increase with 





EFFECTS OF OXYGEN CONCENTRATION ON FAST REACTOR FUELS 


time applies to slab growth, not spherical growth. For 
the spherical case, as in other geometries, the rate of 
increasing thickness of the product layer varies in- 
versely with the total thickness (dx/dt =k/x). This 
leads to parabolic growth or at dependence for both 
the growth of product and fuel-pin diameter. Since the 
observed increase is not parabolic, it must not be 
sodium-diffusion controlled. In fact, because of the 
surface attack, intergranular attack, and fuel cracking, 
it is unlikely that any simple mechanism can be 
proposed for this apparently complex reaction. 

If the rate-limiting mechanism is diffusion of 
sodium through an impervious product layer, one 
would expect the product to form a layer of near- 
uniform thickness on the surface of the fuel. Since 
intergranular attack is observed, reaction control by 
sodium diffusion appears unlikely. On the other hand, 
if diffusion of sodium to the fuel—product interface is 
faster than oxygen diffusion within the fuel, then 
intergranular attack could be explained as follows. 

Presumably sodium can move rapidly to the fuel— 
product interface, i.e., it is limited only by the rate of 
reaction at the interface. The movement of sodium 
may be by diffusion or vapor transport depending on 
the nature of the Na,;MO, film. The slow step is 
assumed to be oxygen diffusion from the fuel to the 
fuel—product interface. Since oxygen is believed to 
diffuse more rapidly along grain boundaries, more 
reaction occurs along intergranular paths than normal 
to fuel surfaces or grain surfaces. Under these con- 
ditions the length of sodium-diffusion paths can show 
considerable variation (compare the surface-layer thick- 
ness to the circuitous length of the intergranular 
material in Fig. 5). Presumably, sodium diffuses rapidly 
through these paths, whereas oxygen moves slowly to 
the grain boundaries. These arguments refer to the case 
where the oxygen for the reaction is furnished by the 
fuel. Further work is in progress to establish the 
mechanism and rate of reaction between sodium and 
fuel. 


Control Methods 


The most reliable method for controlling the 
sodium—fuel reaction is to limit the oxygen both in the 
fuel and in the sodium to a level where the reaction 
product will not be stable. Present oxygen concen- 
tration specifications for the sodium appear to meet 
this criterion on the basis of calculated concentrations, 
except at the lower temperatures (i.e., <773°K) where 
the oxygen in the sodium could produce some 
Na3;MO,. The measured equilibrium oxygen concen- 
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trations are also below specifications so that oxygen 
from the sodium could cause a reaction. 

The oxygen levels necessary in the fuel are not as 
clear-cut because of uncertainties concerning the 
oxygen redistribution in the fuel. The uncertainties are 
emphasized by the in-pile reactions of two defected 
pins encapsulated with sodium (B9A), as discussed by 
Murata et al.” One pin contained fuel with an initial 
O/M ratio of 1.95. The first pin swelled 10%, and the 
second pin swelled 21%. The measured swelling was at 
least four times that calculated by Aitken et al. Part of 
the discrepancy between calculated and measured 
volume expansion is assumed to be the result of 
“nonuniform swelling.” Presumably, nonuniform swell- 
ing means the formation of the product in a narrow 
sector such that it forces the pin open without filling 
all the void space that is created. So far no experi- 
mental evidence exists for this phenomenon. 
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TECHNOLOGICAL CONSIDERATIONS IN 

EMERGENCY INSTRUMENTATION 

PREPAREDNESS. 

Phase Il-A, Emergency Radiological and 

Meteorological Instrumentation Criteria 

for Reactors 

Authors: J. M. Selby, B. V. Andersen, L. A. Carter, J.P. 
Corley, C. E, Elderkin, R. L. Kathren, L. C. Schwendiman, 
C.M. Unruh, E.C. Watson, and R.K. Woodruff, 
Pacific Northwest Laboratory, Richland, Wash. 

Abstract: Criteria for several instrument and instrumentation 
systems are presented. These include: 


— 


. Instrumentation for monitoring liquid and gaseous radio- 
logical effluents, both within the facility and in the 
environs. 

2. Instrumentation for monitoring ambient radiation. 

. Portable radiological monitoring instrumentation. 

4. Meteorological instrumentation, both for synoptic and 

predictive purposes. 


Ww 


The criteria are for instrumentation for determining the 
nature and extent of a radiological incident immediately 
after its occurrence. Emergency instrumentation criteria are 
examined for the contemporary reactor, in consonance with 
what is needed, reasonable, and achievable with the state of 
the art, while at the same time identifying areas in which 
additional study is indicated. Before criteria for reactor 
emergency instrumentation could be developed, it was 
necessary to establish the maximum fission-product inven- 
tory that could be involved in an accident and also to make 
certain assumptions about the fraction of the inventory 
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which might escape from the core. Once the fission products 
were characterized, instrumentation requirements for quickly 
assessing the release and predicting the trajectory and 
airborne concentration downwind were defined. On the basis 
of the quantity and characterization of airborne fission 
preducts which could escape, guides to instrument selection 
are given along with criteria for performance and placement, 
and some techniques are presented for estimating the 
magnitude of the release. Instrumentation for meteorologi- 
cal-parameter measurement is considered in detail because of 
the important role such measurements will play in determin- 
ing the direction of flow and dispersion in the atmosphere, 
and in the initial response by monitoring personnel. 115 ref- 
erences. (auth) 

Available as BNWL-1635. Date of report, May 1972. 


DESIGN STUDIES OF A MOLTEN-SALT REACTOR 

DEMONSTRATION PLANT 

Authors: E. S. Bettis, L. G. Alexander, and H. L. Watts, 
Oak Ridge National Laboratory, Oak Ridge, Tenn. 

Abstract: The 350-MW(e) Molten-Salt Reactor Demonstration 
Reactor is based on much of the technology that was 
demonstrated by the Molten-Salt Reactor Experiment. The 
cylindrical vessel (26 ft in diameter by 26 ft high) houses a 
matrix of graphite slabs forming salt passages having a 
volume fraction in the core of 10%. The flow of fuel salt is so 
distributed that the temperature rise along any path is the 
same—from 1050 to 1250°F. In the primary exchanger, 
heat is transferred to barren carrier salt (in at 900°F and out 
at 1100°F). In the secondary exchanger, heat is transferred 
to a stream of Hitec salt (in at 800°F and out at 1000°F). 
The Hitec oxidizes tritium to tritiated water, which is 
removed and disposed of. The Hitec generates steam at 
900°F, 2400 psi in a boiler, superheater, and reheater. 
Electricity is produced at an overall efficiency of 36.6%. 
Soluble fission products are removed by discarding the 
carrier salt every 8 years after recovery of the uranium by 
fluorination. Volatile fission products are removed by sparg- 
ing the fuel salt with helium bubbles in the reactor primary 
system. The fuel-cycle cost was estimated to be 0.7 mill/kWh 
for inventory, 0.3 mill/kWh for replacement, and 0.1 mill/ 
kWh for processing, giving a total of 1.1 mills/kWh. (auth). 

Available as ORNL-TM-3832. Date of report, June 1972. 
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340 SELECTED RECENT OR RECENTLY RELEASED PUBLICATIONS 


NUCLEAR FUEL CYCLES: A SUMMARY REVIEW nuclear fuel cycle are discussed. Topics discussed include: 
Author: R. G. Wymer, Oak Ridge National Laboratory, source materials, uranium enrichment, reactor requirements, 
Oak Ridge, Tenn. plutonium recycle, fuel reprocessing, fuel refabrication, and 
Abstract: The major aspects of the nuclear fuel cycle are safety aspects. (D.C.C.) 
described. Future developments and requirements for the Available as TID-26071. Date of report, 1971. 





SYMPOSIUM ON ENVIRONMENTAL BEHAVIOR OF RADIONUCLIDES 
RELEASED IN THE NUCLEAR INDUSTRY 


Aix-en-Provence, France, May 14—18, 1973 


It is expected that the symposium will emphasize the presentation of new data and information related to 
the behavior and fate of radionuclides released to terrestrial and freshwater ecosystems. Of particular 
interest will be papers considering the environment as a whole, including the complex interactions of all 
components and sectors of the environment, which lead to an estimation of the impact of radionuclide 
releases on man and other sensitive elements of the environment. 

The symposium will be sponsored by the International Atomic Energy Agency, the Nuclear Energy 
Agency of the OECD, and the World Health Organization. 

For details, persons in the United States should write Mr. Jack Kane, Office of Information Services, 
U.S. Atomic Energy Commission, Washington, D.C. 20545; persons outside the United States should 
write Miss C. de Mol van Otterloo, Administrative Secretary, Division of External Relations, International 
Atomic Energy Agency, Karntner Ring 11, P. O. Box 590, A-1011, Vienna, Austria. 
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Recent USAEC publications of special interest 
RADIOACTIVE WASTE PROCESSING AND DISPOSAL [TID-3311(Suppl. 3) ] 


The 2249 references cited in this bibliograph’ supplement those in TID-3311 and TID-3311(Suppls. 1 and 
2). The references were selected from Nuclear Science Abstracts, Vol. 21 to issue 4 of Vol. 26 (1967 to 
February 1972), by a computer search technique and are arranged by volume and abstract number. 
Corporate and personal author, report number, and subject indexes are included. All literature cited is 
publicly available. 


TITLE LISTING OF CIVILIAN POWER REACTOR DOCKET LITERATURE IN 
NUCLEAR SCIENCE ABSTRACTS, Volumes 21—25 (1967-1971) (TID-3324) 





This bibliography cumulates, in docket number order, 4938 technical documents (announced in NSA) 
submitted to the USAEC in support of U.S. nuclear power reactor license applications. Docket numbers 
are cross referenced to VSA abstract numbers. A subject index is included. 


Available from TID-3311(Suppl. 3) 
Paper, 8", x 11, 592 pp. $6.00 
National Technical Information Service Microfiche 0.95 
U. S. Department of Commerce 
Springfield, Virginia 22151 TID-3324 


Paper, 8'/, x 11, 358 pp. $6.00 
Microfiche 0.95 
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AEC Critical Reuiew Series 


BOILING CRISIS 
AND CRITICAL HEAT FLUX 


L. S. Tong 


Westinghouse Electric Corporation 
August 1972 


The ratio of the predicted critical heat flux (CHF) to 
the reactor heat flux is used as an index of the power 
capability of water reactors. Numerous boiling heat- 
transfer and two-phase-flow studies over the past 15 
years have put emphasis on exploring the boiling- 
crisis mechanisms and on perfecting the CHF predic- 
tions. Unfortunately, the problem is more compli- 
cated than was thought. Thus far no overall analytical 
solution has been obtained, and the design equations 
for CHF predictions are empirical. In this review 
various types of boiling crisis are categorized accord- 
ing to four major flow patterns: pool boiling, bubbly 
flow, annular flow, and normal flow across a cylinder. 
The physical mechanisms of these boiling crises are 
described. Experimental observations are summarized 
and used for verifying various analyses. The missing 
links in the analyses are identified, and future work is 
recommended. 


92 pages, 7 by 9%, inches, paperbound 
Library of Congress Catalog Card Number: 72-600190 
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